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Abstract 

This manual reviews the derivation of the similitude rela- 
tionships believed to be important to ice accretion and 
examines ice-accretion data to evaluate their importance. 
Both size scaling and test-condition scaling methods 
employing the resulting similarity parameters are de- 
scribed, and experimental icing tests performed to evalu- 
ate scaling methods are reviewed with results. The mate- 
rial included applies primarily to unprotected, unswept 
geometries, but some discussion of how to approach 
other situations is included as well. The studies given 
here and scaling methods considered are applicable only 
to Appendix-C icing conditions. Nearly all of the ex- 
perimental results presented have been obtained in sea- 
level tunnels. Recommendations are given regarding 
which scaling methods to use for both size scaling and 
test-condition scaling, and icing test results are described 
to support those recommendations. Facility limitations 
and size-scaling restrictions are discussed. Finally, ap- 
pendices summarize the air, water and ice properties 
used in NASA scaling studies, give expressions for each 
of the similarity parameters used and provide sample 
calculations for the size-scaling and test-condition scal- 
ing methods advocated. 

1. Nomenclature 

cii constants in curve fit for water vapor pressure 

over ice, where i ranges from 0 to 6 
A empirical constant in expression for Nu, 

dimensionless 

A c accumulation parameter, dimensionless 

AOA angle of attack, ° 

b relative heat factor, dimensionless 

bi constants in curve fit for water vapor pressure 

over water, where i ranges from 0 to 6 
B empirical power dependency of Nu on Re, 

dimensionless 
c airfoil chord, in 

c p a constant-pressure specific heat of air, Btu/lbm R 

c pws specific heat of water on model surface, 

Btu/lbm R 

c pjs specific heat of ice on model surface, Btu/lbm R 
c Ka constant-volume specific heat of air, Btu/lbm R 

Ca capillary number, dimensionless 

C D drag coefficient of drop, dimensionless 

Cf skin-friction coefficient, dimensionless 


d cylinder diameter or twice the airfoil leading- 

edge radius, in 

D v diffusivity of water vapor in air, ft 2 /s 

el, e2, e3 powers in supplemental similarity 

parameter,/ 5 , equation (3.65), dimensionless 
E exponent in latent heat of evaporation expres- 

sion, dimensionless 

h c convective heat transfer coefficient, Btu/hr fi 2 R 

hfiim water-film thickness, in 

he gas-phase mass transfer coefficient, lbm/hr fi 2 

k thermal conductivity, Btu/hr ft R 

ki constant in equation (3.15) - (3.17) 

K inertia parameter (eq. (3.5)), dimensionless 
K Bragg’s trajectory scaling parameter, dimen- 

sionless 

K n modified inertia parameter, dimensionless 

l span- wise length of control volume, in 

L undefined length proportional to model 

diameter or chord, in 

L WC liquid water content of cloud, lbm/ft 3 
m mass flux of water per unit time, lbm/ft 2 s 

m e mass flux of water evaporated, lbm/ft 2 s 

m, mass flux of ice sublimated, lbm/ft 2 s 
M Mach number, dimensionless 

n freezing fraction, dimensionless 

no freezing fraction at stagnation, dimensionless 

n a analytical value of freezing fraction, dimen- 
sionless 

n e freezing fraction found from leading-edge 

thickness of experimental ice shapes, 
dimensionless 

Nu Nusselt number based on cylinder diameter 
or twice the airfoil leading-edge radius, 
dimensionless 

Nu c Nusselt number based on airfoil chord, 

dimensionless 
p pressure, psi 

Pni vapor pressure of water over ice, psia 

p mv vapor pressure of water over liquid water, psia 

p w vapor pressure of water in the atmosphere, psia 
P supplemental similarity parameter, equation 

(3.65), dimensionless 

Pr Prandtl number of air, dimensionless 

q c surface heat loss due to convection, Btu/hr ft 2 

q con d surface heat loss due to conduction, Btu/hr ft 2 

q e surface heat loss from evaporation, Btu/hr fi 2 

qf surface heat gain from release of latent heat of 
fusion, Btu/hr ft 2 
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< 7 , surface heat gain from ice cooling to surface 

temperature, Btu/hr ft 2 

qk surface heat gain from kinetic energy of water 
drops, Btu/hr ft 2 

q r surface heat loss from radiation, Btu/hr ft 2 

q,b surface heat loss due to water run-back, 

Btu/hr ft 2 

q s surface heat loss from sublimation, Btu/hr ft 2 
q w surface heat loss to raise temperature of imping- 
ing water to freezing, Btu/hr ft 2 
r recovery factor, dimensionless 

r ie leading-edge radius, in 

r A evaporative-to-convective heat transfer term, 

dimensionless 

R a gas constant for air, lbf ft/lbm R 
Re Reynolds number based on cylinder diameter 

or twice the airfoil leading-edge radius, 
dimensionless 

Re c Reynolds number based on airfoil chord and 
free-stream temperature, dimensionless 
Refiim Reynolds number based on airfoil chord and 
fdm temperature, dimensionless 
Re re i Reynolds number of drop relative to airstream, 
dimensionless 

Res Reynolds number of drop, dimensionless 
s surface distance from stagnation, in 

Sc Schmidt number of air, dimensionless 
t temperature, °F 

tbi temperature of air in boundary layer, °F 

tfiim film temperature, °F 

T absolute temperature, R 

u air velocity vector, dimensionless 

U a local air velocity, mph 

Ug local drop velocity, mph 

U w water-film velocity, mph 

V free-stream velocity (air speed), mph 

We c Weber number based on model length dimen- 

sion, dimensionless 

Wes Weber number based on water drop diameter, 
dimensionless 

Wei, Weber number based on water-film thickness, 

dimensionless 

We, t Weber number based on water-film thickness 

and water-film edge velocity, using Feo correla- 
tions, dimensionless 

We L Weber number based on length L, dimen- 
sionless 

We t Weber number based on water-film thickness, 

using Kind analysis, dimensionless 
We w Weber number based on water-film thickness 
and water- film edge velocity, using Feo analy- 
sis, dimensionless 

We wK Weber number based on water-film thickness, 
using alternate Kind analysis, dimensionless 


x horizontal coordinate, in 

X drop position vector, dimensionless 

y vertical coordinate, in 

p catch efficiency, dimensionless 

Pn catch efficiency at stagnation, dimensionless 

y ratio of specific heats for air, 1 .4 

8 drop median volume diameter (MVD), pm 

A ice thickness, in 

A C p change in pressure coefficient on airfoil, 

dimensionless 

s emissivity of surface, dimensionless 

6 air energy transfer parameter, R 

0 dimensionless time 

k exponent on Reynolds number for {C^yRel 24) 

dependency (eq. 3.15), dimensionless 
X drop range in absence of gravity, ft 

XI Xstokes drop range parameter, dimensionless 
Xstokes drop range in absence of gravity if Stokes' law 
applies, ft 

Afi latent heat of freezing, Btu/lbm 

A s latent heat of sublimation, Btu/lbm 

A v latent heat of vaporization, Btu/lbm 

jU viscosity, lbm/ft s 

p density, lbm/ft 3 

cr Stefan-Boltzmann constant, .1714xl0' 8 Btu/hr 

ft 2 R 4 

o wa Surface tension of water against air, I b t 7 ft 

r icing time, min 

$ drop energy transfer parameter, R 

8, chord-wise dimension of control volume, in 

General Subscripts: 

0 stagnation value 

a air 

c based on airfoil chord 

/ at the freezing point of water 

1 ice 

R reference conditions 

s at the surface 

S scale conditions 

st static 

tot total 

w water 
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2. Introduction 

Icing poses a hazard both to aircraft on the ground pre- 
paring for flight and after take off. Ice can accumulate 
on upper surfaces of aircraft on the ground due to the 
freezing of precipitation. The ice layer is sometimes 
invisible but always introduces additional risks at take- 
off, because it typically increases drag and decreases lift. 
This type of icing is protected against by spraying de- 
icing fluids, which depress the freezing temperature, on 
the upper surfaces of wings and tail planes prior to take 
off. The scaling discussed in this manual applies only to 
in-flight icing encounters, and ground icing will not be 
discussed further. 

In-flight icing results when aircraft fly through super- 
cooled clouds. Supercooling is the state in which water 
exists as a liquid at a temperature less than 32°F. Sus- 
pended water drops in clouds are frequently supercooled 
when ambient temperatures are below freezing. Super- 
cooled water freezes when it impacts surfaces, and the 
resulting ice accretions can have a dangerous effect on 
aircraft. Aircraft flying through supercooled clouds are 
susceptible to ice formations on the leading edges of 
engine inlets, tail planes, wings or rotors and instrumen- 
tation probes. Engine inlet icing can distort the airflow, 
and when shed ice is ingested compressor blades can be 
damaged and combustors may experience flame-out. On 
aerodynamic surfaces icing generally decreases lift and 
increases drag. Iced probes produce false instrument 
readings leading to possible pilot confusion and some- 
times fatal mistakes. Large transport aircraft are pro- 
tected against ice by bleeding hot air from the engine 
compressor to maintain critical surfaces warm enough to 
vaporize water. Some surfaces may not need to be pro- 
tected, and smaller aircraft often use intermittent impulse 
methods to remove small amounts of ice repeatedly. 

Aircraft and component manufacturers must thoroughly 
test new products to determine the effect of icing on their 
performance. This testing is performed both during the 
design process and for certification purposes. Flight- 
testing is necessary but is expensive and can only be 
done when atmospheric icing conditions exist. Further- 
more, it can be very time consuming to find in nature the 
extremes in the cloud drop size and liquid-water content 
envelope required for certification testing. Icing wind 
tunnels can simulate natural icing with water-spray and 
refrigeration systems and provide control of cloud condi- 
tions, temperature and airspeed to permit safe, conven- 
ient and relatively inexpensive testing. Some measure- 
ment of lift and drag changes can be made in the icing 
tunnel, and ice shapes are usually recorded. Knowing 
the ice shape which results from testing at given simu- 
lated flight and cloud conditions, researchers can obtain 


more precise aerodynamic-penalty information either in 
flight or in an aerodynamic wind tunnel. This can be 
done by attaching a reproduction of the ice shape fabri- 
cated from suitable materials (e.g., wood, plastic) to the 
leading edge of the airfoil. 

A scaling method is a procedure to determine the scaled 
test conditions to produce the same result as exposing the 
reference model to the desired cloud conditions. Be- 
cause of test-section blockage limitations, many compo- 
nents cannot be tested full size in an icing wind tunnel. 
Furthermore, facilities that simulate natural icing can 
provide only limited ranges of air speed, cloud drop size, 
and liquid-water content. When the reference (full-size) 
model is too large for a given facility, model-size scaling 
is applied, and when the desired test conditions are out- 
side the facility operating capability, test-condition scal- 
ing is required. 

The objectives of such scaling are usually to insure that 
the amount of ice relative to the model size, the shape of 
the ice accreted, dimensionless ice accretion limits, and 
the aerodynamic penalties due to the ice are the same as 
would have been obtained with the desired, or reference, 
model size or test conditions. It has always been as- 
sumed that for adequate scaling, characteristic features, 
such as glaze horns, need to be simulated in size, loca- 
tion, angle and shape, but how closely these characteris- 
tics need to match between scale and reference accretions 
has never been shown. For testing of ice-protection sys- 
tems, the main objectives may be to match scale and 
reference non-dimensional impingement limits or water 
loading. In any case, it is important to establish reliable 
methods to scale model size or test conditions. In this 
manual, the emphasis will be on ice-accretion scaling 
rather than scaling for ice-protection systems. 

Accurate scaling requires that similitude of geometry, 
flowfield, drop trajectory, drop catch and heat transfer 
are satisfied. The similitude analyses result in equations 
that can be solved for the scale icing conditions that will 
simulate the reference icing encounter. However, these 
relationships may sometimes lead to conflicting values 
for one or more of the test conditions; thus, the similitude 
equations cannot always be rigorously applied, and sim- 
plifications are usually necessary to permit their use in 
practical scaling situations. To address this problem a 
number of scaling methods have been derived from the 
basic set of similitude equations with different simplifi- 
cations and assumptions. These methods provide guid- 
ance for the researcher in establishing proper scaling 
conditions. 

For the testing of airfoils, when the desired model size 
cannot be accommodated in a test facility, an alternative 
approach to attempting to apply simplifications to the 


NASA/CR— 2004-212875 


3 



similitude equations is to use a hybrid scaling method 1 ’ 2 ’ 3 . 
In this approach, a full-size leading-edge portion of a 
wing section is attached to a flapped section of reduced 
size. The flap is adjusted such that the flowfield around 
the leading-edge region of the scale model is the same as 
that of the full-size airfoil. In reference 1, this flowfield 
match was accomplished by trial and error. References 2 
and 3 presented a method in which the design of the test 
article behind the full-size leading edge was determined 
using computer flowfield analysis tools. In the hybrid 
method, scale test temperature, airspeed, liquid-water 
content, median volume drop diameter and icing time are 
the desired full-size values. However, multiple models 
need to be fabricated for different flow conditions and 
angles of attack. The hybrid method will not be ad- 
dressed in this manual. 

Efforts to establish icing scaling methods in which either 
the entire model is scaled in size or test conditions are 
scaled began in the 1950’s and continue to the present. 
In general, to test scaling methods, an ice shape is re- 
corded for a reference condition that represents the full- 
size or desired test, the scaling equations are applied to 
find the appropriate scaled conditions from the reference 
values, and the scaled ice shape is recorded. The two ice 
shapes are then compared. If size has been scaled, the 
comparison can be facilitated if the ice-shape coordinates 
are normalized with respect to some characteristic di- 
mension of the model; for example, the chord of an air- 
foil. Whether the two shapes agree is frequently depend- 
ent on subjective judgment, and the quality of agreement 
when the match isn’t perfect has always been difficult to 
define. 

In this manual, the similitude relationships are reviewed, 
a number of scaling methods are described and com- 
pared, and experimental icing tests performed to evaluate 
scaling methods are reviewed with results. The primary 
emphasis in the studies discussed has been on insuring 
that scale and reference ice accretions have the same 
non-dimensional quantity and shape. The analyses and 
test results reported are for unprotected, unswept geome- 
tries, but some discussion of a preliminary study of inter- 
cycle-ice scaling is included as well. The icing condi- 
tions for which the tests and analyses of this manual were 
made fall entirely within the FAA Part 25 Appendix-C 
envelope, 4 which includes drop sizes from 15 to 50 pm 
median volume diameter (MVD). Supercooled large 
drop conditions (SLD), for which drop sizes are larger 
than 50 pm MVD will not be discussed because studies 
of the physics required to understand phenomena related 
to SLD icing are just beginning. For some SLD situa- 
tions (low airspeed and drop sizes not far above the Ap- 
pendix-C envelope) the scaling methods described will 
be adequate. In general, however, the scaling described 


here should be considered valid only for Appendix-C 
conditions. Although some of the scaling methods dis- 
cussed deal with altitude scaling, nearly all of the ex- 
perimental results presented have been obtained in sea- 
level tunnels. Recommendations are given regarding 
which scaling methods might be most effective for both 
size scaling and test-condition scaling. Restrictions on 
scale size ratio are discussed with respect to practical 
considerations, including facility limitations. Finally, 
recommendations for additional study are outlined. 


3. Similitude Analysis 

To achieve similarity in ice accretion, it is necessary for 
the scale test to simulate the geometry, the flowfield, the 
drop trajectories, the total water catch, the heat transfer 
and, probably, the surface phenomena of the desired ic- 
ing encounter. The similitude analyses performed for 
scaling studies through the 1970’s was made with the 
objective of deriving a series of closed-form equations 
that could be solved for the scale test conditions. Fewer 
simplifications need to be made if equations can be pro- 
grammed for computer solution, and most of the results 
to be given here used the latter approach. 

3.1. Geometric Similarity 

The alternative to the hybrid scaling mentioned in the 
Introduction is to make the reference and scaled models 
geometrically similar over the entire model. In this man- 
ual, the discussion will assume geometric scaling. As ice 
grows on the model, the shape must continue to be simi- 
lar for flowfield similarity to be maintained. This basic 
requisite will be assumed satisfied in subsequent deriva- 
tions from other similitude requirements. 

3.2. Flowfield Similarity 

Similarity of the flowfield would suggest that the Rey- 
nolds and Mach numbers for the scaled test need to be 
matched to their respective hill-size, or reference, values. 
The usual definitions for these parameters will be used: 


Re„ 


Vd Pa 

Pa 


(3.1) 


M 


a 


V 

lyRX, 


(3.2) 


In equation (3.1) the length scale d represents either the 
diameter if the model is a cylinder or twice the leading- 
edge radius if the model is an airfoil. For many analyses, 
the chord of the airfoil has been chosen as the character- 
istic length; however, the main ice accretion occurs in the 
region near the leading edge. Therefore, in studying 
scaling similarity it seems reasonable to define similarity 
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parameters in terms of a length scale representative of the 
leading-edge region. The convention of using d to repre- 
sent either cylinder diameter or twice the airfoil leading- 
edge radius will be used in the definition of all similarity 
parameters in this manual. 

Icing occurs in the atmosphere only within the limited 
range of absolute temperature, T, of 420 to 492°R (233 to 
273 K). Thus, because y and R a are constant properties 
of air, equating the Mach number for the scaled and ref- 
erence cases requires that the velocities for the two cases 
be very nearly the same. However, to match the Rey- 
nolds numbers when the scale model is, for example, half 
the reference size, it is necessary for the scale velocity to 
be approximately double the reference. Clearly, then, 
Mach and Reynolds numbers cannot be simultaneously 
matched when the scale model is other than full size. 

Most scaling analyses avoid this dilemma simply by dis- 
regarding both of these parameters. The rationale for this 
approach is that, for most icing conditions, the Mach 
number is relatively low and compressibility effects can 
be neglected. The Reynolds number has been ignored by 
arguing that ice accretions occur mostly at the leading 
edge of surfaces where the boundary layer is initially 
thin, and viscous effects are small. Farther downstream, 
as ice accretes it often causes the boundary layer to tran- 
sition to turbulent flow, whose characteristics are nearly 
independent of Reynolds number. However, when the 
effects of surface-water phenomena are considered (see 
section 3.6), the Reynolds number may need to be in- 
cluded in the analysis. 

In addition to Mach and Reynolds number considera- 
tions, simulation of the flowfield requires that the veloc- 
ity, pressure and temperature distributions over the scale 
model must simulate those of the reference case. How- 
ever, scaling methods proposed to date do not explicitly 
attempt to satisfy these requirements. In spite of this, 
some success with scaling methods has been achieved, 
probably because significant ice accretions are limited to 
the leading-edge region. 

3.3. Drop Trajectory Similarity 

The mass of water reaching each part of the surface of 
each model must be similar. The drop trajectory deter- 
mines if an individual drop will impinge on the surface 
of the model and, if so, where on the body. Thus, simi- 
larity of drop trajectory between scale and reference 
cases must be satisfied. An analysis of the drop trajecto- 
ries will be made in this section, and expressions will be 
derived which relate the scaled and reference values of 
some of the operating and test parameters of importance. 
The quantity of water reaching the surface will be dis- 
cussed in section 3.4. 


Several simplifications are required to derive workable 
equations. It is assumed that the trajectories in both the 
scaled and the reference clouds can be adequately repre- 
sented by that of one typical drop diameter for each. 
Bragg, et al 5 showed that for the typical distributions of 
drop size and relative velocity of Appendix-C icing ex- 
periments, the drop motion in the flowfield and eventual 
impingement on a surface can be accurately represented 
by using the drop MVD. Note that this simplification 
may not be valid for SLD drop sizes, for which the drop 
distribution may have an effect on the accreted ice shape. 
The small size of cloud water drops involved in Appen- 
dix-C icing minimizes the effect of gravity relative to 
other forces; therefore, gravity will be neglected. The 
non-dimensional drop momentum equation describing 
the motion of a single drop relative to the airstream is 
then given by (see, for example, Bragg 6 and Ruff 7 ) 


d 2 X _ C D Re rd ( dX\ 
d & 2 24 K \ dGJ 


(3.3) 


where the relative drop Reynolds number, Re re i, is based 
on the local relative velocity between the air and drop: 

S Pa 


dX 
u 


Pa 


In equation (3.3), K is the non-dimensional inertia pa- 
rameter defined by Langmuir and Blodgett 8 as 


K = 


pj 2 v 

18 d/i a 


(3.5) 


In Langmuir and Blodgett’s paper the cylinder radius 
was used in place of the dimension d. Thus, for cylin- 
ders, the radius should be used in equation (3.5). Bragg 6 
and Ruff 7 used a length c which they defined as a char- 
acteristic dimension. Airfoil chord is commonly used in 
evaluating the inertia parameter for aircraft icing applica- 
tions. For scaling, however, the use of chord means that 
for a general trajectory analysis the same inertia parame- 
ter could result for airfoils of the same chord but different 
forms, whether thick or thin. By setting d to twice the 
leading-edge radius, this problem is avoided. Maximum 
airfoil thickness might also be considered for the charac- 
teristic length, but the use of twice the leading-edge ra- 
dius leads to collection efficiencies consistent with those 
from the LEWICE ice-accretion code, 9 as will be shown 
below. 


Because equation (3.3) is non-dimensional, if the refer- 
ence and scale flows approaching the models are similar, 
the trajectories are the same for both reference and scaled 
cases if 
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and 


C p Re rei 1 c D R e r A 

24 J s 24 J R 


(3.6) 


*s = r r (3.7) 


C D Re r J 24 is the ratio of the actual drag on a sphere, C D , 
to the Stokes-law drag, 24 IRe re i. It is a function only of 
Reynolds number. Thus, if the Reynolds number and 
inertia parameter, K, for the scaled case were matched to 
the values for the reference case, drop trajectory similar- 
ity would be satisfied. Unfortunately, in practical situa- 
tions it is not always possible to satisfy both parameters 
simultaneously. The modified inertia parameter, K 0 , 
introduced by Langmuir and Blodgett 8 , offers a solution 
to this problem. It combines the effects of the inertia 
parameter and the Reynolds Number and thereby pro- 
vides a single similarity parameter to satisfy. The Lang- 
muir and Blodgett analysis will be presented next, then in 
section 3.3.2 a simplified form will be given. 

3.3.1. Langmuir and Blodgett Trajectory Analysis 

Langmuir and Blodgett’s expression for modified inertia 
parameter is 

K 0 =- + -?— I AT-— |, forW>- (3.8) 

8 8j 8 

In equation (3.8) A Slokes is the range of the drop if it were 
released in still air at an initial velocity of V with the drag 
given by Stokes' law, and A is its range with the actual 
drag. The dimensionless range parameter, A! A Stokes was 
defined by Langmuir and Blodgett as the average drag 
ratio over the Reynolds number range, 0 to Re s , 
experienced during the drop trajectory 


2 _ 1 R ? d Re rd 

Aaoh. Re s o C D Re re ,/ 24 


(3.9) 


The drop Reynolds number, Re s , is given by: 


R e s 


VZpg 

/L 


(3.10) 


It can be seen that the range parameter is a function only 
of drop Reynolds number. 

By using the parameter Kg, equations (3.6) and (3.7) can 
be combined and replaced by equation (3.11): 


K os =K or (3.11) 

Langmuir and Blodgett noted that for K < 1/8 drop im- 
pingement would not occur. This fact is reflected in 


equation (3.8), which is only valid if K is greater than 
1/8. Figure 3.3.1 shows plots of inertia parameter as a 
function of V/d for drops of 4 - 10 pm. The K = 1/8 
limit is noted on the figure. Several runs of the LEWICE 
ice-accretion code 9 were made with drops of 4, 5, 6 and 7 
pm. Although the code’s drop trajectory routine has not 
been validated for drops smaller than 10 pm, its results 
were reasonably consistent with the Langmuir and 
Blodgett impingement lower bound. At a V/d of 7.7 
(m/s)/cm LEWICE predicts no accretion for 5-pm drops, 
accretion only at the leading edge (J5 n = 1 .4%) for 6-pm 
drops, and accretion with a leading-edge collection effi- 
ciency of 2.7% for 7-pm drops. The Langmuir and 
Blodgett limit would suggest no accretion would occur 
below a drop size of 7 pm. At a V/d of 19.3 (m/s)/cm 
LEWICE predicts no accretion for 4-pm drops, but ac- 
cretion with a leading-edge collection efficiency of 2.7% 
for 5-pm drops. The Langmuir and Blodgett limit also 
indicates no accretion with a 4-pm drop, but that accre- 
tion should occur with a 5-pm drop. 



V/d, (m/s)/cm 

Figure 3.3.1. — Drop Impingement for Small Values of 
Inertia Parameter. For K < 1/8 drops will not impinge on 
the model. Solid symbols are the result of LEWICE 2.0 9 
calculations showing some impingement. Open symbols 
are the result of LEWICE 2.0 calculations for which no 
impingement occurred. All LEWICE calculations were 
for a 1.83-m chord NACA 0012 airfoil ( d = 5.78 cm). 
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The range parameter, AJ&stokes, in equation (3.8) was tabu- 
lated by Langmuir and Blodgett as a function of the drop 
Reynolds number, Reg. A fit to the Langmuir and 
Blodgett tabulation is 


X r0.8388 + 0.001483Re/ 

K tokes ^+0.1847 ^Re s 


(3.12) 


Equation (3.12) gives range parameters within ±0.2% of 
Langmuir and Blodgett’s values for 14 < Res < 600. The 
error increases to 0.4% for a Reg of 800. For a static 
temperature of 0°F, Reg for icing encounters range from 
18 to 760, for airspeeds from 100 to 400 mph, drop me- 
dian-volume diameters from 10 to 50 pm and air pres- 
sures of 6.8 to 14.7 psia (corresponding to altitudes from 
20,000 ft to sea level, respectively). 



The local catch efficiency, or collection efficiency, of the 
model is defined by considering the projection of a 
stream tube from the undisturbed flow upstream of the 
model onto the model surface at the location of interest. 
The stream tube contains water drops with the distribu- 
tion of sizes and the liquid water content of interest. The 
fraction of the original water content that actually im- 
pacts the model is the local catch efficiency. The smaller 
drops in the distribution will tend to follow the air 
streamlines around the model, while the larger drops will 
turn less easily, and are more likely to strike the surface. 
The catch efficiency, fi varies over the model surface 
with a value pn at the stagnation line. Langmuir and 
Blodgett 8 published tables of the stagnation-point catch 
efficiency as a function of the inertia parameters, K and 
Kn for cylinders. They showed that fairly accurate values 
of p n could be calculated from equation (3.13) if the iner- 
tia parameter, K, is less than or equal to 7.5: 

L 40 Kf 

Po = V 84 (3- 13 ) 

i+L 4 r- 8 J 

where K n is the modified inertia parameter from equation 
(3.8). A sample set of conditions which would produce a 
value of K 7.5 at 0°C is d= 1.0 in (2.5 cm), V = 330 
mph (147.5 m/s) and S= 20 pm. For an NACA 0012 
airfoil, d= 1 in corresponds with a chord of 31.6 in (80.3 
cm). Larger model sizes, lower velocities or smaller 
drops would result in smaller values ofK. We will see in 
the next paragraph, however, that equation (3.13) gives 
values of p 0 that are consistent with LEWICE calcula- 
tions over a wide range of conditions. 

Figure 3.3.2 compares the collection efficiency calcu- 
lated from equation (3.13) with values found from the 
LEWICE ice accretion code. The LEWICE code deter- 


Modified Inertia Parameter, Kq 

Figure 3.3.2. — Stagnation Collection Efficiency for 
NACA 0012 Airfoils at 0° Angle of Attack. Static Tem- 
perature, 10°F; Static Pressure, 14 psia; Airspeed, 120 to 
300 mph; Water Drop Median Volume Diameter, 10 to 
50 pm; Liquid-Water Content, 1 g/m’. Open Symbols, 
7-in Chord; Shaded Symbols, 21-in Chord; Solid Sym- 
bols, 31. 5-in Chord. Data Represented by Symbols are 
from LEWICE 9 Predictions. 

mines the local catch efficiency by integrating the equa- 
tions of motion of the drops assuming they are rigid 
spheres acted upon by drag and gravity forces. This 
process gives the trajectories, thus determining the drop 
impingement locations on the model surface. By com- 
paring the release locations of several adjacent drops in 
the undisturbed flow upstream with the trajectory inter- 
sections on the model surface, the local collection effi- 
ciency can be determined. The code was run with 
NACA 0012 models with chords of 7.0 - 31.5 in (17.8 - 
80.0 cm), velocities from 120 - 300 mph (53.6 - 134.1 
m/s) and mono-dispersed drop sizes of 10 - 50 pm. Fig- 
ure 3.3.2 shows that equation (3.13) gives leading-edge 
collection efficiencies that closely match those found by 
the more detailed procedure of the LEWICE code. The 
equation (3.13) value for p n ranged from 4.7% lower 
than the LEWICE value at K 0 = 0.4 to 0.23% higher at K n 
= 34. The conditions considered in figure 3.3.2 provide a 
range of inertia parameter, K, of 0.72 to 202, suggesting 
that Langmuir and Blodgett’s upper limit of K = 7.5 for 
equation (3.13) is very conservative. 

Figure 3.3.3 compares LEWICE determinations of lead- 
ing-edge collection efficiency at angles of attack of 0 and 
10° with those from equation (3.13). For the NACA 
0012 airfoil used for these computations Po varies little 
from 0 to 10° AOA. If the scale model is mounted at the 
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o ... „ „ . -v.u -vj.iw w.uu vj.ikj 

Figure 3.3.3. — LEWICE Predictions of Effect of Angle 

of Attack on Collection Efficiency for NACA 0012 Air- Distance From Stagnation Line, s/c 


foils of Various Chord. Static Temperture, -30 to 25°F; Figure 3.3.4._LEWICE-Predicted Reference and Scale 

Static Pressure, 14 psia; Airspeed, 100 and 400 mph; Collection Efficiencies. NACA 0012 Airfoils at 0° 

Water Drop Median Volume Diameter, 10 to 100 pm; A0A . Reference Conditions: c R , 21 in (53.3 cm); V R , 

Liquid-Water Content, 1 g/m . Open Symbols, AOA, 0°; 150 .0 mph (67 m/s); S R , 30.0 pm. Scale Conditions: c. s -, 

Solid Symbols, AOA, 10 . 10 .5 in (26.7 cm); V s , 262.5 mph (1 17 m/s); 8s, 15.6 pm. 


identical AOA as the reference, matching of K 0 will still 
produce the correct drop trajectories for the scale test 
although there may be a small error in fio if equation 
(3.13) is used. 

Langmuir and Blodgett also gave a fi n expression for 
cylinders for K > 7.5. However, because equation (3.13) 
has been validated in figure 3.3.2, there is no reason to 
use any other expression for wing icing situations. 

Equation (3.13) shows that to match fi n between scale 
and reference values, it is only necessary to match K 0 . 
Because K n is easy to calculate, it has been universally 
adopted in scaling methods as the similarity parameter to 
use to insure similarity of drop trajectories. 

The initial capture efficiency is assumed to be that at the 
stagnation line of a clean surface. It can be argued for 
scaling, however, that aft of the stagnation line, because 
scale and reference models are geometrically similar, the 
collection efficiencies must vary in the same way for 
both models. This agreement of scale and reference col- 
lection efficiencies is shown in figure 3.3.4, which gives 
LEWICE predictions for two NACA 0012 airfoils of 
different sizes. The conditions for the smaller model 
were scaled from those of the larger such that the stagna- 
tion K n (and therefore fin) for the two matched. The two 
curves are indistinguishable over the range for which 


accretion occurs. Therefore, to help identify the curve 
for the reference case the area under that plot has been 
shaded. Other reference conditions and their resultant 
scale conditions have also been tested with LEWICE. 
The scale and reference beta curves matched consis- 
tently. Consequently, it is only necessary to match the 
stagnation K n to properly scale drop trajectories over the 
entire clean airfoil. It is assumed that as ice accretes, 
because the geometry changes in the same way for both 
models, the time-varying collection efficiency will con- 
tinue to match everywhere. 

Derivation of the equations for trajectory similitude as- 
sumed that a single drop diameter, the median volume 
diameter, adequately represents the cloud and its trajec- 
tory. At any instant during icing, the size of a drop im- 
pacting a particular location on the model surface cannot 
be predicted. Thus, in glaze ice conditions in particular, 
features of the ice shape tend to be subject to a certain 
randomness, with the history of drop impact at each sur- 
face location having some small effect on the final accre- 
tion there. The resulting surface variations mean that 
icing encounters repeated with the same model and iden- 
tical test conditions cannot be expected to produce iden- 
tical ice shapes. 

Noticeable macroscopic-scale differences in ice shapes 
produced with repeated test conditions are usually minor, 
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however. Although Olsen and Walker 10 reported a strong 
effect of intermediate ice shape on final shape, others 
have not observed as much influence. Anderson and 
Shin 11 found that 10-min glaze ice shapes were not af- 
fected by initial A-min ice accretions at different test 
conditions. In these tests, the final conditions determined 
the shape. Hansman and Turnock 12 found that ice shapes 
were independent of the initial surface material. Appar- 
ently, for glaze ice, initial roughnesses features, whether 
inherent in the uniced surface or produced by a short- 
time ice accretion, are masked by subsequent ice deposits 
and eventually have little or no noticeable effect on the 
final shape. This is fortunate, since there is no easy way 
to include the random effect of drop size in the scaling 
equations. The cloud drop distribution and individual 
drop impingement have been found to be adequately 
represented by the median volume diameter, MVD. Fur- 
thermore, characterization of the initial uniced surface is 
not necessary when dealing with substantial ice accre- 
tions. 

3.3.2. Simplified Trajectory Analysis 

For making quick calculations or estimates, a closed- 
form expression for scale drop size can be developed by 
applying two simplifications to Langmuir and Blodgett’s 
analysis. The first is to define the modified inertia pa- 
rameter as: 


K 0 =K — 

/l Stokes 


(3.14) 


The second simplification is to develop an approximate 
expression for C- D Re rd t24 that can be easily integrated. 
Figure 3.3.5 is a plot of CD^e re /24 as a function of Re rd . 
Although in general the data cannot be fit by a simple 
power law, over limited ranges of relative Reynolds 
number the relationship 


C pRe rel 

24 


= k x Re 


K 

rel 


(3.15) 



Reynolds Number, Re yel 

Figure 3.3.5. — Ratio of Actual Drag, C D , to Stokes-Law 
Drag, 24/Re re i. Drag values from Langmuir and 
Blodgett. 8 


are substituted into (3.14), the modified inertia parameter 
in terms of the test conditions and physical properties is 
obtained: 


K 0 = 


r: 


18^(1 -k)h 1 ; 


g2—K-y\ — Krj-rK 


d p K 


(3.17) 


The density of water (/>„,) and the viscosity of air vary 
little over the relatively narrow temperature range (-40°F 
to 32°F) at which icing occurs. Thus, these parameters 
have been included in the brackets along with the con- 
stants in equation (3.17). When the scale and reference 
values of K 0 are equated, the term in brackets cancels, 
and an expression is obtained which approximately satis- 
fies drop trajectory similarity and from which the scale 
drop size can easily be determined: 


can be used to represent the data. k } and k are constants. 
Figure 3.3.5 shows that the k = 0.38 power approxi- 
mately describes the C D Re rd /24-vs-Re re i data over the 
range of relative Reynolds number from 5 to 100, and k 
= 0.45 is satisfactory for the range 10 < Re re i < 300. 
When equation (3.15) is substituted into the expression 
for the range parameter, equation (3.9), and the integra- 
tion performed, the range parameter becomes: 


1 


k Stokes 



1 

0 -*) 


Res 


(3.16) 


When the definition of Re s from equation (3.10) is in- 
cluded in equation (3.16) and equations (3.5) and (3.16) 


Ss 

Sr 


f ds ) 

i2-/c 

( \ 

Ps 

2-k j 

f 

\2—k I 

( T \ 

l s 

\d r ) 


[pr) 


Ur) 


k T R; 


(3.18) 


This expression has been widely used with different val- 
ues for k. Olsen 13 at NASA used k= 0.38, Charpin and 
coworkers 14 ’ 15 at ONERA and Sibley and Smith 16 at 
Lockheed used k= 0.39 and Dodson at Boeing 17 used k 
= 0.45. Equation (3.18) can be used to obtain the scale 
drop size that gives an approximate match to scale and 
reference K n . Alternately, equation (3.11) can be solved 
using equations (3.8), (3.5) and (3.12). Depending on 
the software used, this latter approach may require an 
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iterative solution for which equation (3.18) can be used 
to provide an initial guess of scale drop size. 

As a practical matter, for typical sizes and conditions for 
wing icing, K n found from equation (3.14) will differ 
from that using equation (3.8) by less than 5%. The scale 
8’s found from equation (3.18) are almost always within 
1 pm of those using equations (3.11) and the analysis of 
section 3.3.1. 

Bragg 6 presented some alternative approaches. The first 
was to replace the simplification of equation (3.15) with 
an integratable expression that represents the drag curve 
over the entire Re r( , t range. The result was a closed-form 
expression for K n : 


K 0 =\&K\ 


Re, 


-y/6 Re 


r Re s ' /3A 


(3.19) 


Although equation (3.19) may represent greater accuracy 
over a wider range of Res, than the method that led to 
equation (3.18), finding scale drop MVD by matching 
the scale and reference K 0 from equation (3.19) is not as 
convenient as using equation (3.18), and we have already 
seen that equation (3.18) produces results with accept- 
able accuracy. 


Bragg also proposed an alternative scaling parameter to 
K 0 . It was derived by applying the drag approximation of 
equation (3.15) directly to the drop trajectory equation 
(3.3) instead of using it in the averaging process that de- 
fined the range parameter, equation (3.9). Bragg defined 
his trajectory scaling parameter as 


K = — (3-20) 

Re/ 

Because Re re i (equation (3.4)) and Res (equation (3.10)) 
differ only by the dimensionless relative velocity be- 
tween the air and drop, equation (3.3) was now rewritten 
with the single similarity parameter of equation (3.20). 
Thus, to satisfy drop trajectory similarity it is only neces- 
sary to match scale and reference values of K . Match- 
ing scale and reference values of K from equation 
(3.20) also results in equation (3.18); thus, K differs 
from the approximate K n of equation (3.17) by only a 
constant. Bragg also found that a value of /cot' 0.35 best 
represented icing conditions, consistent with most of the 
studies cited following equation (3.18). 

3.3.3. Effect of MVD on Ice Shape 

It is necessary to establish the importance of each simi- 
larity parameter and test parameter. In this section we 
will review the effect of drop size on ice shape. In 1998 
Chen ls performed icing tests to compare shapes pro- 
duced with drop MVD’s of 20, 36 and 55 pm. The tests 


were made in the IRT using a 61-cm-chord GLC 305 
airfoil model with constant t st , V, and accretion time. 
LWC was adjusted to maintain constant of 0.3 and 0.5. 
With fixed velocity and model size fin varied with MVD. 
The quantity of ice collected is determined by the prod- 
uct PoA„ and this value was maintained within about 
10%. Results for drop sizes of 55 and 20 pm are given 
in figure 3.3.6. Although j3 n changed from about 92 to 
74% for this change in 8. \ the reduction in drop size and 
collection efficiency appeared to have no effect on the 
main ice shape. Undoubtedly, the impingement limit 
would have changed with drop size (that is, with J3 0 ), but 
this characteristic was not measured. 

Additional ice-shape comparisons from tests using 
NACA 0012 airfoils at 4°AOA were reported in Ander- 
son 19 and are shown in figure 3.3.7. The ice shapes were 
recorded in the NASA Glenn Icing Research Tunnel in 
1997 tests by Harold E. Addy and 1998 tests by Colin S. 
Bidwell and Judith K. Van Zante. The digitized data are 
available from a 2-CD set accompanying Wright. 20 As 
with the GLC 305 ice shapes of figure 3.3.6, the shapes 
compared in figure 3.3.7 reveal no significant effect of 
drop MVD on the main ice shape. 

The implications of these results for scaling are signifi- 
cant. Within the Appendix-C envelope, the main ice 
shape appears to be virtually independent of drop size, 
and, thus, Kg. Therefore, if calculated conditions for a 
scale model test involve drop sizes outside the available 
range for a particular test facility a different 8 can be 
substituted providing the key parameters j5gA c and n n are 
maintained at the reference values. The accumulation 
parameter, A c , and freezing fraction, n, will be introduced 
and discussed in the next section. When scaled tests 
need to provide information about impingement limits, 
however, the scale K n must match the reference. If nec- 
essary, impingement limits can be found from tests with 
scale velocity and model size chosen such that K n 
matches the reference. The test conditions for impinge- 
ment-limit tests would likely be different from those 
required to satisfy proper scaling for the main ice shape. 

3.4. Water Catch Similarity 

The quantity of ice accreted at a specified location on the 
model depends on (1) the amount of water in the cloud, 
(2) the portion of that water that reaches the given sur- 
face location (dependent on the drop trajectories as dis- 
cussed in the previous section) and (3) the fraction of 
water having reached the surface that freezes in that loca- 
tion. The concept of the freezing fraction, n , is due to 
Messinger. 21 For some conditions it is possible that there 
is no local freezing of water; in that case, n = 0. On the 
other hand, if all the water that reaches the surface of the 
model freezes on impact, n = 1. The freezing fraction is 
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x/c 

A/ID = 55 pm. Run 624.40 

A/ID = 20 pm, Run 624.43 

(a) n = 0.28 


x/c 

“ Ml '/.) = 55 pm, Run 924.40 

Ml’D = 20 pm. Run 924.43 

(b) n = 0.52 
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(a) 624.40 

24 

14 

200 

55 

1.16 

6.1 

91.6 

2.53 

2.32 

0.28 

0.84 

16.1 

20.3 

10.9 

6.76 

2.49 

2.02 

624.43 

24 

14 

198 

20 

1.31 

6.1 

73.9 

2.83 

2.09 

0.30 

0.76 

15.9 

20.2 

10.8 

2.44 

2.44 

1.98 

(b) 924.40 

24 

2 

200 

55 

1.16 

6.1 

91.7 

2.53 

2.32 

0.51 

0.84 

28.6 

36.8 

11.4 

6.77 

2.56 

2.02 

924.43 

24 

2 

200 

20 

1.30 

6.1 

74.0 

2.83 

2.09 

0.54 

0.76 

28.5 

36.7 

11.4 

2.46 

2.54 

2.00 

Figure 3.3.6. 

— Effect of Drop MVD 

on Ice Shape. GLC 305 Airfoil at 0°AOA. 

Ice-Shape Data from Tests by Chen ls 


another similarity parameter and will be discussed further 
in the following section. The rate of growth of ice with 
time when n = 1 is given by 


dA m 
dr p i 


(3.21) 


In equation (3.21) m is the total mass of water striking a 
unit area of surface per unit time: 


m = LWCV P (3.22) 


Here L WC is the liquid-water content of the cloud and P 
is the local catch efficiency. Although /? can be expected 
to vary with time as the ice shape changes, assuming 
both P and p, are constant as ice is accreted, equation 
(3.21) can be integrated over the accretion time to give 
the ice thickness, A: 


LWCV Pt 
Pi 


(3.23) 


One goal of properly scaling the ice accretion is that the 
local thickness of ice normalized by the characteristic 
model dimension will be the same for the scaled and 
reference cases; thus. 



(3.24) 


Then, 


LWCV Pt'' 
Pi d ) s 


LWCV Pt'' 

Pi d Jr 


The accumulation parameter, A c , is defined as 


(3.25) 


A. = 


LWCVt 
Pi d 


(3.26) 


From the previous section we know that to insure trajec- 
tory similarity it’s necessary to match scale and reference 
values of Langmuir and Blodgett’s modified inertia pa- 
rameter, K 0 . And from equation (3.13) if the K 0 match, 
so do the P„. Figure 3.3.4 indicated that if p n is matched, 
so too are all local catch efficiencies. Therefore, equa- 
tion (3.25) is the same as requiring that equation (3.11) 
be satisfied along with 


A c ,s ~ A c .r (3.27) 

to insure properly scaled ice thickness everywhere on the 
model when n = I everywhere. If it is not possible to 
find scaled conditions that permit a match of K n (and 
therefore /?«), it is recommended that the product of p (l A c 
be matched. Thus, 
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8= 20 pm, Run 303 

8= 15 pm, Run 314 

(a) April 1997 Tests by Addy. 


8= 40 pm. Run 413 

8= 25 pm, Rim 414 

(b) January 1998 Tests by Bidwell and Van Zante. 


Run 

c, 

in 

A, 

°F 

v, 

mph 

3 

pm 

LWC, 

g/m 3 

T, 

min 

Pn, 

% 

A c 


b 

A 

°F 

a 

°F 

Rc a , 

10 4 

We s 

10 3 

We„ 

10 3 

We u 

10 6 

(a) 303 

21 

13 
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20 

0.55 
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75 
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13.0 

3.25 
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2.74 

314 

21 

13 
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15 

0.60 

6.4 

66 
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0.52 
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16.7 

19.5 

13.0 

2.44 

3.40 

2.74 

(b) 413 

21 

12 
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40 

0.50 

7.7 

88 

1.54 

0.51 

0.38 

17.7 

20.9 

13.0 

6.50 

3.40 

2.74 
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21 

13 
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25 

0.55 

7.0 

80 

1.54 

0.48 

0.38 

16.7 

19.5 

13.0 
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Figure 3.3.7. — Effect of Drop MVD on lee Shape. NACA 0012 Airfoil at 4°AOA. Ice-Shape Data from Tests by 
Addy, April 1997 and Bidwell and Van Zante, January 1998, reported in data CD’s accompanying Wright. 20 


Po,sA,s=Po,rA,r (3-28) 

should be applied in this case. The use of equation (3.28) 
is not as desirable in general as satisfying equations 
(3.11) and (3.27) separately, and should only be used 
when no other alternatives are available. Only when the 
two collection efficiencies, /3 0 ,s and J3 0 ,r, are within 10%, 
can equation (3.28) be expected to give acceptable 
agreement of the quantities of the scale and reference ice 
accretions. 

For glaze ice, the freezing fraction is less than unity and 
varies with position on the model. The local accumula- 
tion depends not only on what impinges but is the result 
of a mass balance that includes evaporation and surface 
water flowing into and away from the location of inter- 
est. Thus, for glaze it is not simple to quantify the in- 
coming liquid water except at the stagnation line, where 
there is no incoming surface water flow. The general 
non-dimensional ice thickness at the stagnation line de- 
pends on three things: the potential ice accumulation, A c , 
the collection efficiency, fto, and the fraction of water 
reaching the surface that freezes there, no- Thus, 

A j = n 0 A,A, (3.29) 

d 


The freezing fraction, no, represents the heat balance at 
the ice surface; as such it is an additional similarity pa- 
rameter that needs to be matched between reference and 
scale situations. It will be discussed in the next section. 
With both n 0 and fin matched, when the scale non- 
dimensional ice thickness at stagnation is equated with 
the reference, equation (3.27) again results. With the 
additional assumption that p, is the same for scaled and 
reference cases, equations (3.26) and (3.27) can be rear- 
ranged to solve for the ratio of scale-to-reference accre- 
tion time: 


fj_ = 


(v \ 

v s 

' LWC s ' 

t r 

1 ^3 

V 

\ y R J 

LWC 

V R J 


For airfoils, of course, the ratio ( c^Cr ) can be used here 
in place of ( ds/d R ) because these ratios have the same 
value. 

3.5. Energy Balance Similarity 

If all the water impinging on a surface were to freeze on 
impact, equations (3.11) (or (3.18)) and (3.30) would be 
sufficient to satisfy scaling. This situation applies to 
rime ice formation. However, the fraction of water that 
freezes, n, is not always unity, and energy-balance simi- 
larity must be satisfied for glaze ice accretions. 
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The rate at which water freezes on a surface depends on 
the local heat balance. Tribus, et al 22 studied heat trans- 
fer for a heated cylinder in icing conditions. Their analy- 
sis is useful for thermal deicing or anti-icing systems. 
Messinger 21 later performed the classical energy-balance 
analysis for freezing at an unheated surface. Messinger's 
work was further developed by Ruff. 7 The similitude 
relations to be presented here are for an unheated surface 
and are based primarily on the Messinger and Ruff 
analyses. The energy terms for this case, expressed per 
unit area, are listed here. The analyses will be done at 
the stagnation line. 

(1) Heat lost from the surface by convection through the 
boundary layer: 

q c ^ ( T - t/j/) 

(2) Heat lost from the surface due to evaporation of wa- 
ter: 

q e =m e A v 


(11) Heat gained by the surface by water flowing into 
the control volume from a neighboring location. 

These terms will be discussed individually in greater 
detail in the following subsections. 

3.5.1. Convective Heat Transfer 

In term (1) the film coefficient, h c , varies with position 
on the model. For scaling, this term in the heat balance 
will be evaluated only at the stagnation line. It is as- 
sumed that if scale and reference conditions are matched 
there, they will also match everywhere. In addition, the 
heat balance on the clean airfoil (before ice begins to 
accrete) is used; if the clean heat balance matches for two 
models and similarity is maintained in all the physical 
phenomena that affect the shape of the accreting ice, the 
heat balance, although changing with time, will be the 
same for the two models as ice develops. 

Experimental convective heat transfer data is typically 
correlated in the form 2 ’ 


(3) Heat lost from the surface due to sublimation of ice: 


q s = m s A s 

(4) Heat lost from the surface due to radiation: 
q, = a&iT s 4 -Tj) 


(5) Heat lost from the surface to raise the temperature of 
the impinging liquid to the freezing point: 


q» 


= fhc p , ws ( t f 


(6) Heat lost from the surface by water flowing out of 
the control volume to a neighboring location: 

q,b = [(! - n) rh - m e ] c pws (t s -t f ) 


(7) Heat lost from the surface due to conduction through 
the ice into the model: 


q cond 


= ki — 
1 


ts-tst-r 


V 


2c 




(8) Heat gained by the surface due to release of latent 
heat of fusion from the freezing water: 


q f = rim A, 

(9) Heat gained by the surface due to release of sensible 
heat from the ice as it cools from the freezing tempera- 
ture to the surface temperature: 

q t = mnc P ,,s (tf ~ k ) 

(10) Heat gained by the surface from the kinetic energy 
of the water drops striking the surface: 

q k =rii(v 2 / 2 ) 


Nu a =ARe B (3.31) 

where A and B are constants, 

h d 

Nu a = -r— (3.32) 

K 

and Re a is the form defined in equation (3.1). Kreith 24 

gives the following relationship for Nu a for the stagna- 
tion line of a cylinder: 

Nu a =l.l4Pr° A Re° a 5 (3.33) 

where 

c u 

Pr a = (3.34) 

k a 

and Re a was given by equation (3.1). The air properties 
for Nu a , Pr a , and Re a in equation (3.33) are evaluated at 
the film temperature, which is taken as the average of the 
free-stream and surface temperatures, 

^ Aim = + ^<) (3.35) 

For glaze ice the surface temperature is 0°C, and over a 
free-stream temperature range of —40 to 0°C, the range of 
interest to icing, the Pr a evaluated at film temperatures 
varies only from 0.708 to 0.705; thus, 1.14 Pr a A has a 
nearly constant value of 0.992. Then 

Nu a = 0.992 Re 0 / (3.36) 
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Equation (3.33) was used to find the convective heat 
transfer coefficient by Ruff 7 for both cylinders and air- 
foils. Heat transfer measurements were made along the 
surface of airfoils by Gelder and Lewis" and Poinsatte 26 
and around a cylinder by Van Fossen, et. al. 27 At the 
stagnation line of a clean, smooth NACA 0012 airfoil at 
0° AOA in the NASA Glenn Icing Research Tunnel 
(IRT), Poinsatte found 

Nu ac =6.mRe° a A J 2 (3.37) 

where Nu ac and Re a>c are based on airfoil chord rather 
than twice the leading-edge radius: 

(3-38) 

k a 

and 

Re a , c = ^ iL (3.39) 

Pa 

For consistency with the above expressions for Nu a based 
on twice the nose radius, equation (3.37) can be con- 
verted using (for an NACA 0012) d = 0.0316c. Then 
Nu a ,c = Mr fl /0.0316 and Re ac = Re a /0.0316, and equation 
(3.37) is equivalent to 

Nu a = 1.100Re°' 472 (3.40) 

Poinsatte’s IRT data at 4° AOA correlated with 

Nu a = 0.903Re° 498 (3.41) 

The coefficients in Poinsatte’s study had an uncertainty 
of about ±5%. 

Poinsatte also conducted flight tests for which the turbu- 
lence level can be expected to be less than that in the 
IRT. At the leading edge of a clean, smooth NACA 
0012 airfoil, the Nusselt number was 

Nu a = 2.61 8Re° 388 (3.42) 

Using a turbulence level of 0.5%, the approximate level 
in the NASA Glenn IRT, and for a smooth cylinder with 
no ice accretion, Van Fossen, et al reported 

Nu a = 1 ,244Re°' 477 (3.43) 

Thus, the Poinsatte 0012 IRT data correlations gave co- 
efficients and exponents within 10% of those quoted by 
Kreith for a cylinder. Van Fossen’s clean, smooth cylin- 
der tests for IRT turbulence levels produced a coefficient 
about 25% higher than Kreith’ s and an exponent on Re a 
within 5%. For scale testing in an icing tunnel with tur- 


bulence levels similar to those of the IRT, it is recom- 
mended that Kreith’ s expression for Nu a be used because 
it is consistent with these recent experimental data for 
both cylinder and airfoil models. 

Charpin and coworkers 14 ' 15 used a Nusselt number of the 
form 

Nu a = A Re™ (3.44) 

Exponents on Re a near 0.8 were obtained by both Poin- 
satte and Van Fossen, et al only in regions where separa- 
tion would probably have occurred or for very high tur- 
bulence levels. For example, for an NACA 0012 with 
roughened surface in the IRT Poinsatte reported an ex- 
ponent of 0.809 at a location well aft of the leading edge. 
And Van Fossen, et al found that for a cylinder with 
simulated smooth 15-min horn-glaze ice and with 0.5% 
turbulence an exponent of 0.798 occurred just aft of the 
horns. For the same test, near the leading edge the Nu a 
dependency was Re°' 62 . With roughness added to the 
ice, the exponent at the measuring station nearest the 
leading edge was 0.68. However, scaling calculations 
are based on parameters evaluated for a clean airfoil at 
the leading edge, so these exponents are not relevant in 
justifying equation (3.44). 

The effect of using Nu a = Ref * on the value of the freez- 
ing fraction will be discussed in section 3.5.14. 

The boundary-layer temperature, t u = t st + rV 2 l2c Pia , 
where r is the recovery factor. The recovery factor has a 
theoretical value of 0.85 for laminar airflow and 0.87 - 
0.88 for turbulent airflow over a flat plate at temperatures 
of interest to icing. 28 Experimental measurements for a 
flat plate give values of 0.88 - 0.90 for both flow re- 
gimes. 28 For a cylinder in crossflow at subsonic speeds, 
the recovery factor based on free-stream conditions de- 
creases from 1 at the stagnation point to about 0.6 at 80° 
(ref. 28). Messinger 21 used r = 0.875 as an average value 
to cover most situations. For the purposes of this man- 
ual, the analysis will be limited to the stagnation line, 
where r= 1, and, therefore, the convection term becomes 
h c (t s - -V 2 /2cp >a ). 

3.5.2. Evaporation 

The mass of water which evaporates, m e , in term (2) is 

7 [Pww Pw ) /n 

m e =h G - (3.45) 

Ps, 

Here the driving force for evaporation is the difference 
between the partial pressure of vapor at the surface, p mv , 
and that in the atmosphere, p w . The vapor pressure over 
water at the surface is assumed the saturation pressure at 
the surface temperature. 
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A compressible form of equation (3.45) was developed 
by Ruff. It is 


= h. 


P m 
T 

St 

1 P,o. 

0.622 T , 


P,o, Pn 

T tot Pst 


(3.46) 


A curve fit for the saturation pressure of vapor over wa- 
ter was given by Pruppacher and Klett 29 for the range of - 
50°C to 50°C and can be found in Appendix A of this 
manual. 

The gas-phase convective mass transfer coefficient, h G , is 
analogous to the convective heat-transfer coefficient. 
The theoretical relationship between h G and h, can be 
expressed in terms of the Schmidt and Prandtl numbers 
as (see, for example, ref. 24): 


useful curve fit for the saturation pressure of vapor over 
ice (see Appendix A). However, sublimation of ice is 
small compared with evaporation of liquid. Therefore, 
this term will not be included in the energy balance 
equation. 

3.5.4. Radiation 

When hot-air or electrothermal anti- or deicing systems 
are not used, the surface temperature and ambient tem- 
perature are not sufficiently different for term (4) to be 
important. Therefore, for this analysis, radiation will be 
neglected. 

3.5.5. Sensible Heat of Water 

In term (5) the impinging mass of water is found from 
equation (3.22) for a general location. Because other 
terms of the energy balance have been formulated for the 
stagnation line the following form will be used 




/ \ 0.67 

' Pr ' 


p,a \ SC a J 


(3.47) 


In equation (3.47) the Prandtl and Schmidt numbers, Pr a 
and Sc a , are those for air, with properties evaluated at the 
film temperature, tf,i m , from equation (3.35). The Prandtl 
number was given by equation (3.34) and the Schmidt 
number is 


Sc a =~^~ (3.48) 

P D 

r a v 

An expression to calculate the diffusivity of water vapor 
in air was given by Pruppacher and Klett and is presented 
in Appendix A. 

Measurements of h G do not agree precisely with the val- 
ues from equation (3.47), but that equation is useful for 
making estimates. At 0°C equation (3.47) gives h(Jh c = 
4.81 lbm R/Btu. Messinger used an experimental value 
of 2.90 lbm R/Btu for this ratio. Results shown in this 
manual were based on the value from equation (3.47). 

3.5.3. Sublimation 

Similar to the mass of ice that evaporates in term (2), the 
mass of ice which sublimes to vapor, m s , in heat-balance 
term (3 ) can be written 

m s = h G y p ™- p ™> (3.49) 

Pst 

where h G is the convective mass transfer coefficient, 
equation (3.47). The driving force for sublimation is the 
difference between the partial pressure of vapor at the 
surface and that in the atmosphere. At the surface, the 
partial pressure of vapor is the saturation pressure of 
water vapor over ice. Pruppacher and Klett provide a 


m = L WC V ji g (3.50) 

The impinging water is assumed to be at the ambient 
static temperature; therefore, the temperature has to be 
increased from the static temperature to the freezing 
temperature before freezing can occur. 

3.5.6. Runback 

Term (6) evaluates the heat carried away from the con- 
trol volume by water run-back. For rime ice n = 1, and 
the term has a value of 0. On the other hand, for glaze 
ice t s = tf, and the term again disappears. Thus, this term 
can be ignored except for the case of thermal ice protec- 
tion for which t s * tf. Scaling for thermal ice protection 
will not be considered in this manual. 

3.5.7. Conduction 

Ruff 7 included an expression for conduction into the 
model surface. Conduction effects can be neglected for 
icing times greater than a few seconds because the ice 
quickly produces an insulating layer over the model sur- 
face. Even for short times, however, conduction simili- 
tude can be satisfied implicitly by requiring that the scale 
and reference models be identical in material and struc- 
ture. Therefore, term (7) will be omitted. 

3.5.8. Latent Heat of Fusion 

The mass flux entering the control volume at stagnation 
is given by equation (3.50). The portion of this incom- 
ing mass that freezes is found by multiplying by the 
stagnation freezing fraction, n n . Then term (8) becomes 

q f =L WC V P 0 n 0 A f (3.51) 
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3.5.9. Sensible Heat of Ice 

Term (9) is valid only for rime ice. For glaze, the surface 
temperature, t„ is equal to the freezing temperature, tf ; 
therefore, this term disappears. 

3.5.10. Drop Kinetic Energy 

Term (10) needs no further explanation. 

3.5.11. Water Inflow 

The contribution due to term (11) cannot be assessed 
easily except at stagnation. There, water enters the con- 
trol volume only through impingement, and this term can 
be ignored in the heat balance. 


Two other parameters often used for convenience are (ft 
and 8 which have dimensions of temperature and relate 
to the drop energy transfer and air energy transfer, re- 
spectively: 

9 i = t f -t s ,-^— (3.56) 

2c 

p,ws 

and 


8 = 


‘ s ‘st 


V 


2c 


P,“ J 


r Pw^_Ptot_Pw_ ^ 
T st T tot Pst 
1 Plot Pwv. 
0.622 T tot T st J 


A. (3-57) 


3.5.12. Energy Balance Equation 

When the surface temperature reaches steady state, the 
net heat transfer is zero. Thus, the energy-balance 
equation at an unheated stagnation surface for glaze ice, 
neglecting terms (2), (4), (6), (7), (9) and (1 1) is 

<lc +c le + Q w ~ Qf + ‘Ik (3.52) 


or 


t. -t„-~ 


V 


2 c 


+ h. 


Pto, P, 


T 

± St 

T, 0 , Pst 

1 

P tot Pww 


u. 


(3.53) 


v 0.622 T tot 

, \ . A . V 1 

+ mc P ,™ Vf ~ t s,)= mn o A f + m -j 


This equation includes compressibility effects for the 
evaporation term. The incompressible form of the en- 
ergy equation is 


V 


2c 


+ hr 


PP J 


Pww Pw 
Pst 


A. 


(3.54) 


+™ c p,ws [tf ~ t st) = ™ n o A f + m -y 


Several similarity parameters can be incorporated to 
simplify the expression of equation (3.53). Tribus 22 in- 
troduced the dimensionless parameter b, known as the 
relative heat factor: 


b = 


me 


p,ws 


LWCVPqC 


p,ws 


(3.55) 


It relates the total heat capacity of the impinging water to 
the ability of the surface to convect heat. The convective 
heat transfer coefficient, h, can be found from equations 
(3.32) and (3.33). 


Equation (3.57) includes compressibility effects. It was 
used by Ruff and will be used in this manual. The in- 
compressible form, based on equation (3.45) has also 
been widely used. It is 


8 = 


‘ S ‘ St 


V 2 


2c 


pp J 


Pst 


A. (3-58) 


If equation (3.47) is substituted into equation (3.58) 
along with Messinger’s experimental value of 0.693 
gK/joule for ( l/c p , a )(Pr/Sc ) 0 ' 667 the form of 8 used by 
Olsen and by Charpin and coworkers results. For aircraft 
icing conditions, compressibility is rarely an issue. 
Therefore, it is not surprising that the values of 8 from 
equations (3.57) and (3.58) are nearly the same. 

The energy balance, equation (3.53), can be rewritten in 
terms of the parameters b, <j> and 8. Then, the energy 
balance for glaze ice becomes 


n 0 = 


f c p,ws 






J 



(3.59) 


This formulation of the freezing fraction includes simpli- 
fications and evaluations that apply strictly to the stagna- 
tion line. It is often the practice in scaling publications to 
omit the subscript 0 from n, but for this manual, the sub- 
script will be retained to emphasize that the analysis is 
for the stagnation line. By using as the heat-balance 
parameter to match for scale studies, the implicit assump- 
tion is that if no,s = Ho.r then n will also match every- 
where. Furthermore, equation (3.59) applies to a clean 
airfoil. As ice accretes, the stagnation freezing fraction 
can be expected to change because the geometry is con- 
tinually changing. Thus, for scaling it is necessary to 
assume additionally that no varies with time in such a 
way that if n os = nn.R at the initiation of spray, freezing 
fractions will continue to match during the course of ice 
accretion. This is not an unreasonable assumption be- 
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cause the goal of scaling is to maintain similar ice shapes 
for scale and reference. 

It was shown earlier that if the stagnation collection effi- 
ciency, /3o, matched for scale and reference, LEWICE 
predicted collection efficiencies that matched every- 
where on the airfoil. Similarly, Oleskiw, Kind and 
McCullough 30 presented LEWICE predictions showing 
that when the stagnation n n matched, reference and scale 
local freezing fractions also matched from the stagnation 
line back to about s/c = ±0.02. This is the region of 
greatest interest to icing because most of the main ice 
shape is accreted in this range. However, it’s important 
to keep in mind that LEWICE predictions of freezing 
fraction have not been validated experimentally the way 
collection efficiencies have, and therefore there is less 
confidence in the LEWICE values of/7. 

Experimental evidence to validate the formulation of the 
expression for stagnation freezing fraction will be pre- 
sented below in section 3.5.2. 


For rime ice, because water freezes on impact, the energy 
balance does not affect the ice shape. Consequently, no 
energy balance needs to be considered when scaling 
rime. However, for glaze ice, the energy balance similar- 
ity is typically satisfied by matching scale and reference 
values of one or more of the parameters no, (j), <9 and b. If 
any three of these are satisfied, it can be seen from equa- 
tion (3.59) that the fourth is, as well. 

Kind 31 ’ 2 used an additional parameter that related the 
evaporative and convective heat transfer rates. He 
matched scale and reference values of a convective-to- 
latent heat ratio that will be called r A \ 


le 




( 

1 + 

l 9c) 


KAp-h,) 


LWCVA 


(3.60) 


/ 


Typically, q e « q c . Then, substituting the definition of b 
from equation (3.55), equation (3.60) becomes 


(be ^ 

ro'~p,ws 


A 


7 J 



constant 



3.5.13. Effect of Temperature unALWC 

Studies with both cylinders and NACA 0012 airfoils 33 ' 34 
have shown how ice shapes change with temperature and 
LWC. Some of the 0012 results will be given here. The 
effect of incremental changes in temperature on ice shape 
is shown in figure 3.5.1. The tests were performed in the 
NASA Glenn IRT using 21 -in- (53.3-cm-) chord NACA 
0012 airfoil models. In each portion of the figure, two 


ice shapes are compared for which the temperature was 
different with all other test conditions held constant. For 
all tests K n was the same to give a fixed value of /C of 
0.85, A c was 1.88 and b was 0.58. As temperature de- 
creased, the parameters n 0 , ^and <9 increased. 

The ice shape for the warmer temperature of each pair is 
shown shaded while the colder of the two is represented 
with a solid line. With each incremental decrease in 
temperature (and consequent increase in n 0 ) the horns 
moved forward and the leading-edge ice thickness in- 
creased (consistent with equation (3.29)). The total 
quantity of ice produced appeared to be the same for 
each pair of profiles. The same trends were observed 
with a 36-in-chord GLC 305 model when temperature 
decreased from 14 to 1°F. Again, the horns moved for- 
ward and the leading-edge thickness increased when the 
temperature decreased. 18 

With constant test conditions c, V, 8, LWC, and r and 
constant similarity parameters K n ( f } 0 ), A c and b for the 
range of temperatures tested in figure 3.5.1, it must be 
reasoned that only the varying parameters no, (j) or 9 can 
be responsible for the changes observed in ice shape. 
Further evidence will be shown next to isolate which of 
these three has the greatest impact on shape. 

Figure 3.5.2 presents results of another series of tests; 
this time, the test conditions c, t st , V, and 8 were main- 
tained constant while LWC varied. The icing time was 
adjusted to keep A c the same from test to test. For these 
experiments, the similarity parameters K 0 (Pd), A c , (j) and 
9 were unchanged as LWC was decreased. n 0 increased 
with decreasing L WC, and b also decreased. Again, the 
model was a 21-in- (53.3-cm-) chord NACA 0012 tested 
in the NASA Glenn IRT. The shaded ice shape in each 
portion of the figure was recorded for the higher L WC. 
As the LWC was decreased from 1.4 to 1.0 g/m 3 (fig. 
3.5.2 (a)) the horns on the main ice shape moved forward 
and the leading-edge thickness tended to increase, similar 
to the effect of decreasing the static temperature shown 
in figure 3.5.1. A further reduction in LWC from 1.0 to 
0.8 g/m 3 (fig. 3.5.2 (b)) produced a small additional nar- 
rowing of the horn angle. 

Finally, figure 3.5.3 gives results of tests in which both 
temperature and L WC were varied in such a way that the 
freezing fraction was maintained constant. The LWC 
was again reduced from 1.4 to 1.0 g/m 3 (fig. 3.5.3 (a)) 
and from 1.0 to 0.8 g/m 3 , (fig. 3.5.3 (b)). The fixed simi- 
larity parameters were K n (Po), A c and n 0 , while b, (j) and 
9 varied from test to test. Ice shapes now were un- 
changed as L WC varied. 

For the range of temperature and LWC tested, these re- 
sults led to the conclusion that neither temperature nor 
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Figure 3.5.1. — Effect of Static Temperature on Ice Shape. NACA 0012 Airfoil at 0°AOA. 




NASA/CR— 2004-2 12875 


18 



0.08 
0.04 

^ 0.00 
-0.04 
-0.08 

-0.04 0.00 0.04 0.08 0.12 0.16 -0.04 0.00 0.04 0.08 0.12 0.16 

x/c x/c 

LWC = 1.4 g/m 3 , n 0 , 0.40 LWC = 1.0 g/m 3 , n 0 , 0.52 

LWC = 1.0 g/m 3 , n 0 , 0.52 LWC = 0.8 g/m 3 , n 0 _ 0.59 

(a) LWC’s of 1.4 and 1 g/m 3 . (b) LWC s of 1.0 and 0.8g/m 3 . 


Run 

c, 

hh 

V, 

S, 

LWC, 

T, 

Po, 

A 

n n 

h 

A 

e. 

Rc a , 

We s 

We a 

We„ 

in 

°F 

mph 

pm 

g/m 3 

min 

% 

A-c 

u 

°F 

°F 

10 4 

10 3 

10 3 

10 6 

(a) 10-25-95/11 

21 

11 

150 

30 

1.40 

5.2 

80.3 

1.90 

0.40 

0.77 

20.4 

28.2 

8.90 

2.08 

1.52 

1.17 

10-25-95/9 

21 

10 

150 

30 

1.00 

7.3 

80.4 

1.90 

0.52 

0.55 

21.2 

29.4 

8.91 

2.08 

1.51 

1.17 

(b) 10-25-95/9 

21 

10 

150 

30 

1.00 

7.3 

80.4 

1.90 

0.52 

0.55 

21.2 

29.4 

8.91 

2.08 

1.51 

1.17 

10-25-95/10 

21 

11 

150 

30 

0.80 

9.1 

80.3 

1.90 

0.59 

0.44 

20.4 

28.2 

8.87 

2.08 

1.51 

1.16 


Figure 3.5.2. — Effect of Liquid-Water Content on Ice Shape. Constant Static Temperature. NACA 0012 Airfoil at 
0° AOA. 




L WC has an effect on ice shape independent of the freez- 
ing fraction. The results of the three test series shown in 
figures 3.5.1, 3.5.2 and 3.5.3 demonstrate not only the 
very strong effect of freezing fraction on ice shape, but 
also the apparent unimportance of the other energy 
parameters by themselves, b, (j> and 0 are related to n 0 
through equation (3.59) but appear to have little or no 
independent effect on ice shape. These results also show 
that when tests must be performed to simulate an L WC 
that cannot be obtained in a test facility, it is only neces- 
sary to adjust temperature in such a way that freezing 
fraction matches the desired value. This approach to 
scaling test conditions is known as the Olsen method. 1 ' 
The Olsen method will be discussed in section 5.2. 1.2. 

This ability to match shape by matching c, V, 5 and n n 
fails if the scale total temperature is too high. Obviously, 
ice will not accrete if the scale total and static tempera- 
tures are above freezing. If the static temperature is 
below freezing while the total is above, ice will accrete 
over some part of the model, but the quantity of ice at the 
leading edge will be less than would result at the same 
freezing fraction but with a total temperature below 
freezing. An example ' with a static temperature of 20°F 
and a total temperature of 32°F is shown in figure 3.5.4. 
Flere, instead of well-defined horns, ridges of ice formed 


adjacent to the stagnation region while there was little 
accumulation at the leading edge. Inspection showed this 
accretion to be wet. Other investigators 36 ' 37 have reported 
that slushy ice forms when total temperatures are near 
freezing. 

The temperatures reported in figure 3.5.4 were averaged 
over the run time. For this test, the variation in tempera- 
ture over time was approximately ±0.5°F. Therefore, at 
times, the stagnation temperature would have been below 
freezing, and at others, above. Away from the stagnation 
line, the recovery temperature would have been some- 
what below the stagnation (total) value. Consequently, 
more freezing was possible aft of stagnation, accounting 
for the ice ridges adjacent to the leading-edge region. 
Another case, for a total temperature of 33°F, was also 
reported in reference 35. The ice shape was very similar 
in appearance to that in figure 3.5.4. 

Unpublished ice shapes recorded by the author in the 
IRT showed that even for temperatures as much as 
4 - 5°F below freezing accretions look different from 
those formed at lower temperatures at the same n 0 . 
Figure 3.5.5 gives an example. Flere shapes obtained at 
different total temperatures but with constant n n are com- 
pared. In figure 3.5.5 (a) ice shapes are shown for total 
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Figure 3.5.3. — Effect of Liquid-Water Content. Temperature Adjusted for Constant n f) . NACA 0012 Airfoil 
at 0°AOA. 


temperatures of 24.9 and 27.6°F. The same model was 
used for both tests and the velocity and drop MVD were 
unchanged. To maintain the freezing fraction at 0.30, the 
LWC was decreased from 1.22 to 0.80 g/m 3 when the 



x/c 

Figure 3.5.4. — Ice Shape for Total Temperature of 
32°F. 35 NACA 0012 at 0° AOA Tested in the IRT, 3-23- 
98 Run 2; c, 10.5 in; t s „ 20°F; V, 264 mph; 8, 21 pm; 
LWC, 0.49 g/m 3 ; r, 4.3 min; p 0 , 85.8; A c , 1.96; n n , 0.22; 
b, 0.29; <j>, 9.1°F; 6, 6.5°F; Re, 7.11xl0 4 ; We, 4.44xl0 3 ; 
We L , 1.80xl0 6 . 


temperature was increased. The two ice-shapes agreed 
within the typical repeatability of shapes for low freezing 
fractions such as this. There were no major differences 
in the horn size, horn location or leading-edge thickness 
between the two. 

In figure 3.5.5 (b) the 27.6°F test shown in 3.5.5 (a) is 
compared with one at 28.8°F. Again, the freezing frac- 
tion was maintained the same for the two tests by lower- 
ing the L WC when the temperature was raised. All other 
test conditions were the same. In contrast to the minimal 
effects of a 2.7°F total temperature change in 3.5.5 (a), 
the increase in total temperature of just over 1°F in 3.5.5 
(b) produced a dramatic change in the shape. For the 
28.8°F shape the leading-edge thickness has decreased 
and ice ridges have appeared adjacent to the stagnation 
region. Horn-like structures were located significantly 
farther aft than those seen at temperatures from 24.9 to 
27.6°F. 

These observations suggest that, with constant model 
size, velocity and drop size, matched freezing fractions 
will insure matched ice shapes only for total tempera- 
tures below 27°F in the IRT. Test results to define this 
temperature limit have not been published for models 
other than the NACA 0012, but there is no reason to sup- 
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Figure 3.5.5 — Effect of Total Temperature near Freezing on Ice Shape. NACA 0012 Airfoils at 0°AOA. 


pose the effects shown in figures 3.5.4 and 3.5.5 would 
be different for other airfoil profiles. The results shown 
are for 10.5-in-chord models, and it is not known how 
the temperature limit might vary with model size. Other 
facilities would be expected to show similar practical 
limitations, as well, although the value of the upper tem- 
perature limit may be different from that in the IRT. 


3.5.14. Validation of the Freezing Fraction 
Expression 


When experimental ice shapes are available, there are 
two ways to determine the stagnation freezing fraction, 
n n . First, an experimental stagnation freezing fraction, n e , 
can be found using measurements from the ice shapes. 
From equation (3.29), 


« = ■ 


d AJ 0 


(3.62) 


where A 0 is the leading-edge thickness of the ice, and j3 0 
and A c are calculated from equations (3.13) and (3.26), 
respectively, using the conditions at which the ice was 
accreted. Second, an analytical stagnation value can be 
determined from equation (3.59), also using the condi- 
tions of the test, 


n 


a 


f c ^ 

p,ws 

\ A fJ 



(3.63) 


No experimental results are needed to calculate n a . The 
freezing fractions found from equations (3.62) and (3.63) 
should agree. 

Anderson and Tsao 38 compared values of n e and n a for a 
number of experimental ice shapes obtained in the IRT 
using NACA 0012 models with chords of 10.5 - 31.5 in 
(26.7 - 80 cm), velocities of 106 - 257 mph (48 - 115 
m/s) and drop MVD’s of 21.5 - 46 pm. The experimen- 
tally-derived freezing fractions are plotted against the 
analytical values in figure 3.5.6. Open symbols were 
measured from ice tracings taken at mid-span (CL) and 
the solid from the tracings 1 in (2.5 cm) above mid-span 
(CL+2.5 cm). The solid line represents perfect agree- 
ment of the two freezing fractions and the shaded band 
indicates the limits for ±10% agreement. A linear fit to 
the data is shown as a dashed line. 

There was significant uncertainty in the ice thickness 
values found from tracings at low freezing fractions. The 
data scatter in figure 3.5.6 at the lowest freezing fraction 
tested was a consequence of this uncertainty. The linear 
fit to the data fell within the ±10% agreement band, al- 
though on average n e values tended to be about 5% lower 
than the n a . The value of n e found from equation (3.62) 
is inversely proportional to the value of L WC through the 
accumulation parameter (eq. (3.26)). Because the uncer- 
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Freezing Fraction from Analysis, n a 

Figure 3.5.6. — Experimental (eq. (3.62)) and Analytical 
(eq. (3.63)) Stagnation Freezing Fractions of Anderson 
and Tsao. 38 

tainties for the LWC and A c were both estimated to be 
about ±10%, this agreement is quite good. 

While n a = 1 suggests fully rime ice, n e is susceptible to 
the uncertainties in both the LWC calibration and the 
tracing technique and can therefore be greater or less 
than unity. n e values for experimental rime shapes in this 
sUidy were within ±10% of 1. 

Bilanin' 9 performed a similar sUidy using a limited set of 
ice-shape data from Ruff. Bilanin found agreement of n e 
and n a at rime, but the two deviated significantly for 
the lowest freezing fractions. Fie concluded that the 
Messinger analysis for freezing fraction had serious defi- 
ciencies, probably because surface-water effects, particu- 
larly splashing, were not included in the heat balance. 
Flowever, Bilanin did not consider the possibility that his 
results might be explained by experimental complica- 
tions, such as possible problems with the facility LWC 
calibration and the uncertainty in the leading-edge thick- 
ness for low freezing fractions. 

Unlike Bilanin’s results, the experimentally-determined 
freezing fractions for the Anderson and Tsao study in 
figure 3.5.6 showed no systematic deviation from the 
analytical values at low freezing fractions. This consis- 
tently good agreement between n e and n„ over a range of 
conditions including both fully rime and fairly warm 
glaze argued against the contention that the Messinger 
freezing fraction formulation neglects any important ef- 
fects, such as surface-water phenomena. On the con- 
trary, for Appendix C icing conditions, the experimental 


results appeared to validate the Messinger freezing frac- 
tion formulation expressed in equation (3.59) (or eq. 

(3.63) ). 

Besides the consistency between n e and n a over the range 
of freezing fractions evaluated, Anderson and Tsao 
showed that when conditions were such that equation 

(3.63) gave n a = 1, fully rime ice resulted, as evidenced 
both by the shapes and the white appearance. The maxi- 
mum temperaUire for the given model size, velocity, drop 
size, and LWC at which rime should occur was deter- 
mined using equation (3.63). Figure 3.5.7 compares ice 
shapes obtained at this temperature, which should have 
produced fully rime, with ones at slightly higher tem- 
peraUires. Ice shapes at these higher temperaUires would 
be expected to deviate from Hilly rime appearance if the 
formulation for n a is correct. 

The rime shape (10-12-02 run 1 1) used for comparison in 
both portions of the figure is shown with a shaded cross 
section. It was accreted at a static temperaUire of-13°F. 
Figure 3.5.7 (a) compares this shape with one produced 
at — 7°F, giving an analytical stagnation freezing fraction 
of 0.91. The outline of this shape is plotted with a solid 
line and no shading. The shape at the lower freezing 
fraction differed from the rime in that a narrow strip of 
ice at the leading edge changed from the characteristic 
white appearance of rime to a slightly transparent form. 
At the same time, the smooth, convex shape of rime was 
replaced with a slight valley right at the leading edge. 
The remainder of the ice aft of this narrow strip main- 
tained the rime appearance. A further increase in tem- 
peraUire to -5°F reduced the analytical stagnation freez- 
ing fraction to 0.86 (fig. 3.5.7 (b)). This change in n n 
caused the leading-edge depression to deepen slightly 
without changing the appearance aft of the leading edge. 

These results show that the use of the Messinger analysis 
leads to a consistency between the analytical freezing 
fraction and experimental ice shapes for rime and nearly- 
rime ice. 

Anderson and Tsao also looked at the effects of calculat- 
ing the analytical freezing fraction with 2 minor modifi- 
cations. The first was to use the model chord instead of 
twice the leading-edge radius, and the second substituted 
the average Nu a around a cylinder 24 

Nu a = 0.0239Re° 805 (3.64) 

for the stagnation value of equation (3.33). Each of these 
changes led to somewhat poorer agreement between n e 
and n a , but more significantly gave analytical freezing 
fractions of 0.82 to 0.89 for the conditions yielding ice 
with the appearance of fully rime ice. 
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Figure 3.5.7. — Sensitivity of Ice Shape to Freezing Fraction at and near n a = 1. NACA 0012 Airfoil at 0°AOA. 


These comparisons illustrate how it is possible to apply 
simple checks of experimental consistency to expose 
problems with particular analytical methods for calculat- 
ing freezing fraction. At the same time, because fully 
rime is easy to recognize, equation (3.62) provides a 
simple way to validate facility L WC calibrations. Condi- 
tions to give n a = 1 can be calculated using equation 
(3.63), then validated by the appearance of the ice actu- 
ally produced. With J3 n determined from equation (3.13) 
and the stagnation ice thickness measured from the ex- 
perimental accretion, when equation (3.26) is substituted 
for A c , equation (3.62) can be rearranged to solve for the 
actual IJVC. 

3.6. Similarity of Surface-Water Dynamics 

Scaling methods developed from the 1950’s through the 
1970’s relied on matching scale and reference values of 
some or all of the parameters K 0 , A c , and n 0 . Given a set 
of reference conditions and choosing scale model size, 
ONERA and Ruff (see section 5.3) used these parameters 
plus <!> and //from equations (3.56) and (3.57) to provide 
five equations to solve for scale temperature, pressure, 
drop MVD, cloud LWC and spray duration. Some of 
these earlier methods permitted the user the arbitrary 
choice of scale velocity. In this section, it will be shown 


that velocity has an effect on ice shape independent of its 
effect on the similarity parameters identified so far. This 
effect of velocity appears to be associated with surface- 
water phenomena, as will be seen. 

In 1988 Bilanin’ 9 presented a Buckingham - n analysis in 
which he concluded that surface-water phenomena had to 
be included in icing scaling methods. Olsen and 
Walker 10 and Hansman, et al 12 ’ 40 ’ 41 studied surface effects 
and surface water during ice accretion, presenting addi- 
tional evidence that these were important phenomena to 
consider in ice accretion. From the close-up photographs 
of these research studies, it was observed that for glaze 
ice accretion unfrozen water on the ice surface tended to 
coalesce to form beads. These beads sometimes were 
swept downstream and sometimes froze in place. 
Bilanin' 9,42 also argued that drop splashing on impact 
might affect the shape of the ice accreted. 

Hansman and Turnock 12 found that when a surfactant 
was added to the icing spray water, the ice shape appear- 
ance and shape changed significantly, with the glaze 
horns moving toward the leading edge. Clearly, then, 
surface tension, and by implication, surface phenomena, 
have a significant role in the physics of ice accretion. 
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3.6.1. Forms of the Weber Number 

Similar observations were made by Bilanin and Ander- 
son. 42 Figure 3.6.1 includes published and unpublished 
data from their studies and compares ice shapes with and 
without surfactant addition to the spray for velocities of 
105, 150 and 209 mph. The shaded ice shapes were ob- 
tained with the NASA Glenn Icing Research Tunnel 
(IRT) spray bar system de-mineralized water supply with 
no surfactant added, and the shape represented by the 
solid line resulted from surfactant addition to the water. 
Other than the addition of surfactant the same test condi- 
tions were used for each pair of tests, as shown in the 
accompanying table. Surfactant addition reduced the 
surface tension to roughly half that of water. Because no, 
Po and A c were the same for each pair of tests, the lead- 
ing-edge ice thickness was also nearly the same. How- 
ever, the included angle between the horns decreased 
dramatically when surfactant was added. Horn angle 
also decreased when velocity increased as can be seen by 
comparing figure 3.6.1(a) with (b) and (c). All the re- 


sults shown for different velocities were from tests made 
at approximately the same freezing fraction. Thus, both 
surface tension and velocity have an effect on ice shape 
independent of freezing fraction. Clearly, then, scale 
velocity cannot be chosen arbitrarily, and a similarity 
parameter dependent on the ratio V eI lo w/a e2 must be in- 
cluded in scaling methodology, where the powers el and 
e2 need to be determined. 

Chen’s 18 results for GLC 305 airfoils, shown in figure 

3.3.6, as well as NACA 0012 results, shown in figure 

3.3.7, show that drop size and collection efficiency ap- 
pear to have little effect on the main ice shape for at least 
much of the range of Appendix-C conditions. Therefore, 
the parameter being sought cannot be dependent on 8. 

The effect of model size was also evaluated by Chen. 
Figure 3.6.2 illustrates that with p n , A c , no, b, </> and 9 
constant a reduction in chord from 91 to 30 cm moved 
the glaze horns rearward. This is the same effect shown 
in figure 3.6.1 for decreasing velocity. 
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Figure 3.6.2. — Effect of Model Size on Ice Shape. GLC 305 Airfoil at 0° AOA. Ice-Shape Data from September 
1998 IRT Tests by Chen. 18 


The previous section showed that the temperature and 
L WC do not have effects on the ice shape independent of 
the freezing fraction. Thus, the effect of these test condi- 
tions are adequately accounted for by n 0 , and the sup- 
plementary parameter does not need to include them, and 
Its general form must then be 

v e, c e 3 

P = constant — — — (3.65) 

a wl a 

This form suggests a Weber number based on chord: 

V 2 co 

We c = ^ (3.66) 

® 'wl a 

Studies by Bartlett 43, 44 and Oleskiw, et al 45 found no 
measurable effect of pressure on ice shape. These 
observations rule out the dependence on air density in 
equation (3.66), making water density a better choice. 
Furthermore, the length may not be chord itself but rather 
some physical characteristic related to the accreting ice 
that is proportional to chord; for example, the water-fdm 
thickness. Because this length is not yet identified, L will 
be used to represent it, with the understanding that 


Lee c 


(3.67) 


or, since the leading-edge radius is proportional to chord, 
Lozd (3.68) 

and the likely form of the Weber number is then 

We, = V Lf> " (3.69) 

<*■„/ a 

The trends apparent in figures 3.6.1 and 3.6.2 show that 
an increase in We L has the same effect on glaze ice shape 
as an increase in n n \ that is, horns move back with an 
increase in the included horn angle as either n 0 or We L 
decrease. Weber numbers with other lengths have been 
considered in the past, such as one based on drop size 
and water properties, We s : 


We s 



(3.70) 


Before data like those of figures 3.3.6 and 3.3.7 were 
available, this Weber number seemed to be the logical 
choice, based on the assumption that drop splashing 


NASA/CR— 2004-2 12875 


25 



plays a role in establishing shapes for glaze ice, and that 
drop size affects splashing characterization. Anderson 
and Ruff 35 , Anderson 46 and Kind' 2 have all used the We- 
ber number of equation (3.70). Because in most studies, 
S was scaled along with chord, matching Weg would lead 
to scale velocities not much different from those obtained 
from matching Wep, thus, positive scaling results from 
these studies do not give us clues about the correct form 
of Weber number. 

From equations (3.68) and (3.69) the scale velocity found 
from We Li s = We LR is 


K 



(3.71) 


As a practical matter for tests with nearly constant air 
density, this result is little different from that obtained 
using constant We c (eq. (3.66)), but has a stronger basis 
in experimental evidence. Equation (3.69) provides one 
potential additional similarity parameter to supplement 
Ruff s basic scaling method. 


Another approach presented by Feo 50 involved simulta- 
neously matching the water-film thickness and a Weber 
number that used the water-film velocity and water-film 
thickness: 


UfafUmPw 
We k = = 

a w!a 


(3.75) 


where the water-film velocity was found by equating 
shear stresses at the air-water interface 


U,„ = 


^ film 


\2 


V 2 Re„ 


(3.76) 


When equation (3.76) is incorporated into equation 
(3.75) and scale and reference values of hjuJd are 
equated, the scale free-stream velocity for this Feo 
method is 


v s =v R 




\ d s J 


(3.77) 


The Weber number based on the water film thickness, 
hflim, was proposed by Feo. 47 


... ^filmPw 
We h = 

a w/a 


(3.72) 


Feo and Urdiales 48 had measured the water film thick- 
ness for heavy rain and developed the following correla- 
tion: 


LWC V2 Re~ V4 Re l d 5l4 Wef 14 (3.73) 
d 

Later studies by Feo 47 were made with L FTC’s in the 
range of 2.6 to 14.5 g/m'. Although still somewhat 
higher than desired to represent typical icing conditions, 
these L JVC’s were lower than those used in the heavy 
rain work. From this effort, the correlation 


^ film 

~d~ 


' LJVC ^ 


Re. 


- 1/4 


(3.74) 


was found. Evaluations of the water-film-thickness We- 
ber number in scaling tests using hfn m from equation 

(3.74) have been made 49 and the results compared well 
with those using the Weber number of equation (3.70). 
While these results are encouraging, the conclusions of 
section 3.5.13 that LWC does not have an effect on ice 
shape independent of freezing fraction are not consistent 
with a Weber-number length having the form of equation 

(3.74) . It’s possible that L may prove to be a water-film 
thickness of some form, however. 


While this approach provides a simple expression to find 
scale velocity, there is no fundamental reason that both 
film thickness and Weber number should simultaneously 
match. Another approach is to equate only the scale and 
reference We k from equation (3.75) without requiring 
film thickness to match. With water-film velocity again 
taken from equation (3.76) and water-film thickness from 
equation (3.74), when the We k are matched, the scale 
velocity becomes: 


r* = v„ 


LWC n 


LWC. 


sj 


5 

W 


\ d S J 


Pa,R 

{Pc,sJ 


1 


Pa,S 

Pa.R 


(3.78) 


• 30 5 1 52 53 

Kind and coauthors ’ ’ ’ have also argued for the use 
of a Weber number based on water-film thickness as a 
similarity parameter for scaling. They introduced two 
alternatives, We w and We„ given by 


We... = 


p a V 2 h 


film, A 


(3.79) 


' wl a 


and 


We,= 


Pa V 2 h 




(3.80) 


a w/a 

where hfi im/1 and hf,i mB are two forms of analytically de- 
rived water-film thickness (t wa and t wh , respectively, in 
Kind’s nomenclature). Each Weber number is similar to 
equation (3.72), but uses air density in place of water 
density. Expressions for Kind’s water film thickness 
were derived by first applying continuity, then develop- 
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ing two different mean water-film velocity terms. The 
first, for hfu m A, was found by assuming the water film 
velocity was the result of aerodynamic shear stress alone 
(Couette flow). The resulting water-film thickness was 


h film, A 


6 LWC' 

1/2 

f \ 
Pw 

\ Pw ) 


V Pa y 


1/2 


Re. 


- 1/2 


(3.81) 


The scale velocity found by matching scale and reference 
Weber numbers using equations (3.79) and (3.81) was 
then 


s-il/3 


Fc =Vn 


lwc r 

JWCsj 


\ d s ) 


Pa,R 

Pa,S 


(3.82) 


As noted for other expressions discussed in this section, 
existing data supports the view that L WC does not have 
an effect on ice shape independent of the freezing frac- 
tion; therefore, equation (3.82) does not seem to be fully 
consistent with experimental observations. However, 
because the effect of LWC on V s in equation (3.82) is 
weak (L WC 1/3 ), this expression cannot be disregarded. 

The mean water-film velocity used by Kind to develop 
an expression for h/u mB came from the additive effect of 
both shear-, as for hfu m A, and pressure -gradient- 
(Poiseuille-flow-) driven contributions. The resulting 
expression for water-film thickness is not a closed form, 
but must be solved by iteration: 




A 


fi w LWCV p(\-n)-smO 
~p/ r ' 2 


c f 

^F + 0.42AC„ 
2 p 


(3.83) 


Yl film, B 


Kind evaluated this film thickness at an angle of 45° 
from the stagnation, and A C p was the change in pressure 
coefficient from stagnation to 45°. 


In addition to the various Weber numbers described here, 
there are several other similarity parameters related to 
water film. These will be discussed next. It is possible 
to find a scale velocity by matching scale and reference 
values of any of the water-film similarity parameters. 


3.6.2. Other Water-Film Parameters and Additional 
Methods to Determine Scale Velocity 

Feo 54,55 suggested that the non-dimensional water film 
thickness itself might be a critical dimension controlling 
surface water dynamics. Although no experimental data 
exist for the actual film thickness in glaze icing 
encounters, the correlation of equation (3.74) was 
obtained for cloud conditions not too much different 
from icing clouds. This form of if considered as a 
similarity parameter, suggests that L WC should have an 


effect on ice shape independent of its effect through n 0 . 
As noted previously, this conclusion is not consistent 
with the observed effects of LWC shown in section 
3.5.13, where ice shapes were seen to be invariant with 
changing L WC as long as no was fixed. Nevertheless, the 
correlation of equation (3.73) was used in two recent 
studies 35,46 as a possible description of hpjd. Scale 
velocities calculated by matching scale and reference 
values of hfi m from equation (3.73) were lower than the 
reference and gave glaze ice shapes that did not match 
the reference. The general conclusion of references 35 
and 46 from looking at a range of scale-to-reference 
velocity ratios was that scale velocities needed to be 
higher than the reference for good simulation of glaze ice 
shapes when the scale model size was less than the 
reference. 


When scale and reference values of water-film thickness 
were equated to find V s using equation (3.74), the scale 
velocity values were typically close to those determined 
with the average- F method discussed next and provided 
generally good matches of scale and reference ice 
shapes. 49 

Equation (3.73) suggests that the water-film thickness 
depends on both the Re a and We s . Anderson and Ruff 3 " 
and Anderson 46 used a scale velocity that was the aver- 
age of those found by matching Wes and by matching 
Re a . This average velocity is 


V s 


Vr d R P a ,R Pa,S 
2 yd s p aS p aR 



(3.84) 


Excellent scaling results were typically obtained in refer- 
ences 35 and 46 using this average -velocity method of 
finding scale velocity. 

Finally, the capillary number, Ca, is the ratio of viscous 
forces to surface-tension forces. This parameter has been 
used by Kind, et al' 1 ' 2 to describe the phenomena in the 
surface water on glaze ice. The capillary number is 

Ca = (3.85) 

a w/a 

For practical scaling, the scale test will use water with 
viscosity and surface tension that match the reference 
values; thus, matching of the capillary number is equiva- 
lent to matching the scale and reference velocity. As 
noted previously, studies of the effect of velocity on sub- 
scale ice shapes 35,46 showed the best match of reference 
ice shapes occurred when the scale velocity was higher 
than the reference. Thus, the capillary number by itself 
is not likely to be a useful similarity parameter. How- 
ever, Kind noted that when the We w (eq. (3.79)) is 
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matched, the capillary number is coincidentally matched 
as well. 

Section 5.3.4 compares the scale velocity that results 
from using each of the parameters presented in sections 
3.6.1 and 3.6.2 when scaling model size. Section 5.3.7 
gives recommendations for how best to determine the 
scale test conditions, including velocity, when model size 
is scaled. 

3.7. Other Phenomena 

The physical models for icing on which the similitude 
parameters have been based have not been fully verified 
and probably do not include all phenomena that may 
have some effect on the ice accretion process. Close-up 
studies of water impact and freezing on surfaces were 
made by Olsen and Walker. 10 They believed that shed- 
ding of water from the surface is a significant part of the 
ice-accretion process. Such shedding is most likely to 
occur from the tips of the horns. Whether such shed wa- 
ter is subsequently re-entrained and deposited aft of the 
main ice shape is not known at this time. However, wa- 
ter shed from the horns would not impact the accretion 
process at the leading edge, and therefore does not affect 
the analytical expression for freezing fraction at stagna- 
tion or the n e found from the measured ice thickness 
there. In addition to water shedding, ice shedding is a 
major consideration for situations involving high speeds, 
rotorcraf! and deicing systems. Olsen and Newton 1 ' and 
Ruff 7 discuss the similitude requirements for ice shed- 
ding. 

Both of these shedding phenomena take place away from 
stagnation where the scaling analyses are focused. Thus, 
the present approach to scaling cannot easily incorporate 
such effects. Furthermore, until high-speed, close-up 
studies can better identify the processes taking place, 
shedding events will not be included in scaling methods. 
For typical Appendix C conditions, reasonable success 
with scaling has been demonstrated without including 
shedding. 

3.8. Summary of Potential Scaling Parameters 

For convenience, all the similarity parameters discussed 
above that appear to have the most influence on ice ac- 
cretion and are therefore most likely to form the basis 
of a scaling method will be listed here again. All are 
dimensionless except for the water-energy and air-energy 
transfer terms ((/) and 6) which have dimensions of tem- 
perature. 

• Reynolds number, Re a (eq. (3.1)) 

• Modified inertia parameter, K 0 (eq. (3.8)) 

• Stagnation collection efficiency, /3 n (eq. (3.13)) 

• Accumulation parameter, A c (eq. (3.26)) 


• Relative heat factor, b (eq. (3.55)) 

• Water-energy transfer parameter, $ (eq. (3.56)) 

• Air-energy transfer parameter, <9(eq. (3.57)) 

• Freezing fraction, n 0 (eq. (3.59)) 

• Weber number, for example: We c We L , Wes, 
We,„ We w or We, (eqs. (3.66) (3.69), (3.70), 
(3.72), (3.79) or (3.80)) 

• Water film thickness, hjiiJd, (eq. (3.74)) 

Note that in place of the modified inertia parameter, K 0 , 
the collection efficiency at the leading edge, J3 0 , can be 
used. Pn is a more sensitive parameter than K„. These 
two parameters are directly related by equation (3.13). 

4. Experimental Evaluation of Scaling Methods: 
Studies at NASA 

Many of the test results cited in this manual were ob- 
tained in the NASA Glenn Icing Research Tunnel (IRT). 
This facility and the experimental procedures used in 
testing scaling methods will be described here. In addi- 
tion, facility effects on scaling results will be discussed. 

4.1. NASA Glenn Icing Research Tunnel 

The NASA Glenn (formerly NASA Lewis) Icing Re- 
search Tunnel (IRT) has been updated and improved 
several times in the last 10 years. Several of these modi- 
fications occurred over the period of time the data pre- 
sented in this manual were taken. Consequently, only a 
general description of the tunnel will be given here. For 
more information, the reader is directed to the relevant 
IRT publications 56,57 ’ 58 ' 59 and references to individual 
scaling studies. 

The IRT has a test section width of 9 ft and a height of 
6 ft. The current refrigeration system and heat exchanger 
allows accurate control of the test-section total tempera- 
ture from approximately -20 to +40°F. 

A water spray system with ten spray bars simulates the 
conditions in a natural icing cloud. The test-section 
cloud drop MVD, 8, and liquid-water content, LWC, 
depend on spray-bar air and water pressures. The rela- 
tionships among these pressures, the tunnel airspeed and 
the cloud properties are established periodically by a 
series of tunnel calibration tests. 60,61,62 For Appendix C 
icing the cloud has been calibrated over a range of test- 
section airspeeds from 50 to 350 mph and drop median 
volume diameters from 14 to 50 pm. Two sets of spray 
nozzles, the Mod-1 and Standard nozzles, are used to 
provide different ranges of liquid-water content. 
Depending on the nozzle set, the airspeed and the drop 
diameter, the test-section liquid-water content can be 
controlled from less than 0.2 to over 5 g/m 3 . 
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Figure 4.1.1 gives the IRT Appendix-C cloud calibration 
envelope for both nozzle sets for a test-section velocity 
of 200 mph. The L WC range increases for lower veloci- 
ties and shrinks for higher velocities. The solid line indi- 
cates the limiting envelope for the Mod-1 nozzles and the 
dashed line that for the Standard nozzle set. These enve- 
lopes illustrate the limitations that can be placed on 
available test conditions in the IRT. For example, test 
conditions that include a velocity of 200 mph, a drop 
MVD of 30 pm and an LWC of 0.4 g/m 3 can be seen 
from figure 4.1.1 to be unavailable. Similar limitations 
apply to all facilities, although the specific boundaries of 
operation may be different from those illustrated. Scal- 
ing methods can be used to determine appropriate condi- 
tions within the envelope that will simulate the results of 
tests at conditions outside the calibrated envelope. These 
methods will be discussed in section 5.2, and sample 
calculations will be presented in Appendix B. 

4.2. Tunnel Effects 

In addition to the limits inherent in the cloud calibration, 
several tunnel effects need to be recognized when per- 
forming scaling studies or evaluating results of scaling 
tests. These effects are common to all facilities, although 
some facilities may be more susceptible than others. 

LWC is typically determined without the test model in 
place. Anything that changes the flow in the tunnel, such 
as blockage that’s not centered or models that generate 
lift, will shift the position of the cloud. Because the 
cloud in an icing tunnel may only be uniform to within 
±20%, movement of the cloud will alter the centerline 



Drop Median Volume Diameter, pm 

Figure 4.1.1. — IRT Calibrated Envelope for 200 mph 
(89 m/s). 


LWC from the calibrated value. Furthermore, cloud size 
and L WC uniformity may vary with airspeed. Clearly, a 
non-intrusive method is needed to measure LWC at the 
test section with the model installed at the desired angle 
of attack. Until such a method is available users need to 
assure themselves of the calibration validity for their 
model and test conditions. Some ways to do this are 
outlined in the following paragraphs. 

The consistency of the cloud calibration over a range of 
conditions should be checked with the model installed. 
This can be done by testing at rime conditions where 
proper scaling can be assured simply to matching A c and 
po between a reference and scale test, with both sets of 
conditions satisfying n 0 = 1. The non-dimensional ice 
shapes from both tests should match if the tunnel calibra- 
tion is consistent for the two model sizes and sets of con- 
ditions. 

For example, figure 4.2.1 reproduces rime test results 
from reference 35. The ice shapes recorded for three 
scale tests (solid lines) are compared individually with 
the same reference test shape (shaded). The spray-bar 
conditions for these tests were chosen to produce the 
desired LWC based on the tunnel calibration. The A, 
reported in the table accompanying the figure was based 
on this calibrated value of L WC. The three portions of 
the figure show scale tests performed at three velocities: 
figure 4.2.1 (a) shows results for V s = 171 mph, (b) for 
200 mph, and (c) for 264 mph. 

Good consistency of the LWC calibration was seen for 
reference and scale test conditions in figure 4.2.1 (a) and 
(b), but the smaller scale ice shape for part (c) indicated 
that at the relatively high velocity for that test, the true 
LWC was apparently lower than indicated by the tunnel 
calibration. Testing with LWC’s significantly different 
from assumed values affects both the quantity of ice 
accreted and, for glaze ice, the freezing fraction. The 
resulting ice shapes can lead to incorrect interpretations 
of scale test results. 

The actual LWC for each test shown in figure 4.2.1 could 
have been determined by measuring the leading-edge 
thickness for each ice shape and applying equations 
(3.29) and (3.26) with the known value of / 3 0 and n 0 = 1. 

Another tunnel effect which can affect ice accretion and 
influence scale test results is the possibility that water 
drops have frozen before impact on the model. For the 
IRT this can be avoided by insuring that the temperatures 
of the air and water in the spray bars are high enough. 
Tests of ice accretion with varying spray -bar tempera- 
tures have been used to establish proper operating tem- 
peratures. A related problem is the failure of water drop 
temperatures to reach the air temperature. This is a 
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Figure 4.2.1. — Rime Scaling Tests to Check Tunnel LWC Calibration Consistency. 35 NACA 0012 Airfoils at 0°AOA. 


particular concern in small tunnels for which drop 
residence time may be too short for ambient temperatures 
to be reached. 

All of these effects need to be carefully checked; other- 
wise, it’s difficult to know if a poor scaling result dem- 
onstrates the inadequacy of the scaling method under 
evaluation or a tunnel effect. 

When performing scaling tests, it is always wise to per- 
form the scale and reference tests during the same test 
entry whenever possible. In addition, the evaluation of a 
particular scaling method should be based on a large 
number of tests including both repeated conditions as 
well as a range of reference test conditions. These pro- 
cedures avoid unwarranted conclusions due to random 
tunnel repeatability differences and help to define types 
of situations for which a scaling method may or may not 
be valid. 

4.3. Tunnel Comparisons 

Phenomena such as water-drop freeze-out or failure to 
reach ambient temperature can affect ice shapes. These 
effects can differ from one facility to another. Some of 


these phenomena are influenced by the physical size and 
features of a tunnel. In addition, the calibrations of L WC 
and 8 arc not presently subject to any standard methods 
and are not based on absolute reference values. Conse- 
quently, ice shapes generated for a given set of nominal 
conditions cannot be expected to reproduce from one 
tunnel to another. Illustrations of how shapes can differ 
due to calibration differences will be given in this 
section. 

In 1996, Anderson and Ruff ran a series of NACA 0012 
icing tests in the IRT. 35 The indicated test conditions 
were repeated with the same models in a second tunnel. 
These “facility-B” test results have not been published 
previously. Recently, results of tests at a third icing tun- 
nel (‘‘facility C”) were published,’ 0 and these used nomi- 
nal conditions close enough to previous IRT test condi- 
tions to provide additional data for facility comparison.. 
Facilities B and C both had a smaller test section than the 
IRT. 

Figures 4.3.1 and 4.3.2 compare the ice shapes from each 
of the facilities. The facility-B-IRT comparisons are 
given in figure 4.3.1, and the facility-C-IRT comparisons 
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x/c 

Facility B: 8-21-96 Run 7 

IRT: 3-26-98 Run 4 

(a) Same Nominal Conditions for Both Facilities. 


Fac. 

Date/Run 

c. 

hh 

v, 

s. 
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T, 

no 



in. 

°F 

mph 

pm 

g/m 

min 

B 

8-21-96/7 

21 

17 

150 

40 

0.65 

11.2 

0.47 

IRT 

3-26-98/4 

21 

17 

150 

40 

0.61 

11.2 

0.48 


x/c 

Facility B: 8-21-96 Run 7 

IRT: 3-26-98 Run 2 

(b) IRT Tests at Higher Temperature than Facility B. 


Fac. 

Date/Run 

c, 

*sh 

V, 

s. 

LWC, 

i 3 
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n 0 



in. 

°F 
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pm 

g/m 

min 

B 

8-21-96/7 

21 

17 

150 

40 

0.65 

11.2 

0.47 

IRT 

3-26-98/2 

21 

22 

150 

40 

0.61 

11.2 

0.30 


Figure 4.3.1. — Comparison of Ice Shapes Accreted in Facility B and IRT. ' 5 All Tests with NACA 0012 at 0°AOA; 
Pth 86%; H c , 1.90. 


in 4.3.2. The shaded shapes in the two figures represent 
those obtained in facility B or C, respectively, and the 
shape with the solid line in each portion of the figures 
was recorded in the IRT. The shapes in figure 4.3.1 (a) 
were obtained with the same nominal test conditions in 
both tunnels. The calculated freezing fraction, based on 
the indicated test conditions, was 0.47 for each. The 
facility B accretion is significantly larger than that from 
the IRT and has an overall shape that is characteristic of 
a lower freezing fraction. Both the increased size relative 
to the IRT accretion and a lower freezing fraction would 
result from a higher L WC compared with the IRT value. 

Figure 4.3.1 (b) compares the same facility B shape with 
one from the IRT at a reduced freezing fraction. For this 
illustration, the IRT freezing fraction was decreased by 
increasing the static temperature by 5°F (3°C) rather than 
by increasing L WC; however, section 3.5.13 demon- 
strated that the effects of LWC changes on ice shape 
characteristics can also be produced by changing tem- 
perature when the freezing fraction is constant. Although 
the quantity of ice is still much greater for the facility-B 
accretion, the features of the two test results corre- 
sponded well. These results indicate that the LWC cali- 
bration for facility B gave a significantly higher LWC 
than that for the IRT. 


Similar results were obtained when IRT shapes were 
compared with those from facility C in figure 4.3.2. 
Again, the nominal conditions in 4.3.2 (a) were the same 
for both facility C and the IRT. As for the facility-B 
shape, the facility-C shape was not only larger than that 
from the IRT, but had the characteristics suggestive of a 
lower freezing fraction. In 4.3.2 (b) the facility-C shape 
is compared with one from the IRT obtained at a lower 
freezing fraction (higher temperature), and the features, 
but not the size, of the two accretions were remarkably 
similar. Apparently, the true LWC for facility C, as for 
B, was significantly higher than that of the IRT for the 
same reported cloud conditions. 

It is clear from these comparisons that the cloud LWC 
calibrations cannot be the same in all three facilities. 
Because of the differences, different conclusions will be 
drawn about what shape will be achieved with a particu- 
lar freezing fraction. However, if these calibrations 
differ simply by a constant factor such that the LWC is 
consistent over the full range of conditions, evaluations 
of scaling tests should be valid in any of the facilities. 
Nevertheless, the true LWC must be validated for each 
facility to insure that users obtain realistic ice shapes and 
quantities. Presently, each facility is calibrated by meth- 
ods used traditionally by personnel in that establishment. 
Standardized calibration instrumentation and methods 
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Facility C: Ref-1 Facility C: Ref-1 

IRT: 11-13-00 Run 3 IRT: 11-13-00 Run 1 


(a) Nominal IRT and Facility C Conditions Almost Same. (b) IRT Tests at Higher Temperature than Facility C. 


c D c, t sl , V, 8 LWC, x, 

Fac. Run . 0 3 . n 0 

in. F mph pm g/m min 

C Ref-1 19.7 14 150 30 1.00 7.3 0.44 

IRT 11-13-00/3 21.0 14 149 38 1.00 7.3 0.40 


c, t s „ V, 8 LWC, x, 

in. °F mph pm g/m 3 min n " 


C Ref-1 19.7 14 150 30 

IRT 11-13-00/1 21.0 19 149 38 


Figure 4.3.2. — Comparison of Ice Shapes Accreted in Facility C 30 and in IRT. 35 NACA 0012 Airfoil at 0° AOA for all 
Tests; p 0 \ 86%; A c , 1.90. 


need to be established, and the resulting calibrations must 
be checked against rime ice shapes as described in 
section 4.2. Documentation of the effects of tunnel flow 
(angularity and turbulence) and blockage on local cloud 
is also needed. 

Efforts are being made to understand tunnel differences 
and improve calibration standards among facilities 
through the SAE AC-9C Icing Subcommittee. 63 Recent 
re-calibrations of the facilities compared in figures 4.3.1 
and 4.3.2 could significantly change the relative shapes 
and sizes of the accretions shown in those figures. 

4.4. Repeatability 

There is a random aspect to ice accretion. Each cloud of 
a given 8 and L WC includes a distribution of drop sizes. 
The size of a water drop hitting a particular spot on the 
model at a specific time into the spray may be different 
for each repeated test run. Small variations in local tem- 
perature, 8 and LWC in the cloud occur with position and 
with time. At the same time, there is some (usually 
small) variability in the test conditions from run to run. 
Because of the random nature and test variability, some 
differences in accreted ice shape from one run to another 
and from one test entry to another can be expected. 

Testing of scaling methods requires comparing ice 
shapes recorded during separate tunnel tests, one using a 
reference model with reference test conditions and one 


using the scale model with scale conditions. It is neces- 
sary, then, to document the extent of the shape variation 
from run to run to determine how closely scaled shapes 
can be expected to agree. Whenever tunnel physical 
features or spray-bar characteristics are changed, the 
tunnel should be re-calibrated. The new calibration 
should permit the V, 8 and L WC at the test-section center 
to reproduce values previously obtained such that tests 
performed with the new calibration will give ice shapes 
in reasonable agreement with those found up to that time. 

Figure 4.4.1 compares ice shapes taken from 1995 
through 1998 in the NASA Glenn IRT using the same 
model and the same test conditions each time. In figure 
4.4.1 (a) a March 1998 shape (shaded) is compared with 
one taken in October 1995. The gross features, including 
the leading-edge thickness, the horn size and the horn 
positions, of the main ice shape repeated well with only 
minor differences,. Apparently, the feathers were not 
traced in the October 1995 test. In general, this compari- 
son shows very good repeatability of the shapes. Good 
repeatability of shapes in the IRT is typical. It suggests 
that the random aspects of icing cited earlier have only a 
small effect on the major features of ace accretions. 

Occasionally, ice shapes do not repeat well from one 
entry to another, however, and the reasons are not well 
understood. This situation is shown in figure 4.4.1 (b) 
where a December 1995 test produced significant 
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x/c 


x/c 


3-26-98 Run 2 
10-26-95 Run 3 


3-26-98 Run 2 
12-11-95 Run 2 


(a) October 1995 Compared with March 1998. (b) December 1995 Compared with March 1998. 

Figure 4.4.1. — lee Shape Repeatability in 1RT. 21-in- (53.3-cm-) Chord NACA 0012 Airfoil at 0°AOA Tested with 
same Conditions on Three Entries. t st , 22°F (-6°C); V, 150 mph (67 m/s); 8, 40 pm ;LWC, 0.65 g/m 3 ; r, 11.2 min. 


differences from that of March 1998. While the leading- 
edge thickness matched, the horn size and angle were 
different for the two shapes. If the accretions in figure 
4.4.1 (b) had been produced from a test of a particular 
scaling method, the conclusion might well be that the 
shapes did not agree sufficiently well to validate that 
method. When evaluating scaling methods, it is impor- 
tant to repeat tests several times with multiple tunnel 
entries to avoid mistaken conclusions. 

Another issue that arises with respect to comparing 



-0.08 H 1 1 1 1 1 1 1 1 1 1 1 1 

-0.04 0.00 0.04 0.08 0.12 0.16 

x/c 


3-26-98 Run 2 CL 
3-26-98 Run 2 CL + 2in 


shapes recorded by hand tracing is the variability of the 
recorded shape due to the inaccuracy of the method. 
Figure 4.4.2 is a comparison of the mid-span tracing of 
an ice accretion with tracings made 2 in (5 cm) above 
mid-span (fig. 4.4.2 (a)) and with one made 1 in (2.5 cm) 
below mid-span (fig. 4.4.2 (b)). The March 1998 ice 
shape from figure 4.4.1 (mid-span position) has been 
used again for this illustration. For all three tracing sites 
the leading-edge thickness and most of the features of the 
main ice shape were reproduced very consistently. The 
horn tips varied somewhat from one tracing to another; 

1 (b) 



— i — i — i — i — i — i — i — i — i — i — i — i 

-0.04 0.00 0.04 0.08 0.12 0.16 

x/c 

3-26-98 Run 2 CL 
3-26-98 Run 2 CL - 1 in 


(a) Tracings at Mid-Span (CL) and 2 in Above. (b) Tracings at Mid-Span (CL) and 1 in Below. 

Figure 4.4.2. — Comparison of Tracings at Different Span-wise Locations in 1RT. 21-in.- (53.3-cm-) Chord NACA 
0012 Airfoil at 0°AOA; t s „ 22°F (-6°C); V, 150 mph (67 m/s); 8 40 pm; LWC, 0.65 g/m 3 ; r, 11.2 min. 
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some of this variation may have been due to actual 
differences over the span of the model and some to the 
pencil tracing technique itself. The feather size and loca- 
tion show the greatest diversity, which reflects the reality 
that feather formation is generally the result of random 
impacts of drops or possibly of water shed from the tips 
of the horns. Thus their size and density vary randomly 
along the span, and it is not realistic to expect scaling 
methods to precisely reproduce the specific location of 
feathers from a reference test, although the approximate 
feather size in a general location on the model should be 
simulated. 

4.5. IRT Scaling Models 

Scaling studies in the IRT have used cylinders, GLC 305 
airfoils and NACA 0012 airfoils made from aluminum, 
fiberglass and wood. Airfoil chords have ranged from 5 
to 36 in. All models were mounted vertically in the cen- 
ter of the IRT test section. Figure 4.5.1 shows typical 
NACA 0012 models in the IRT test section. In figure 
4.5.1 (a) a 6-ft-span, 36-in-chord fiberglass airfoil section 
can be seen, and 4.5.1 (b) shows a 2-ft-span, 21-in-chord 



(a) 91.4-cm-Chord NACA 0012 Model Installed in IRT 
Test Section. 

Figure 4.5.1. — NACA 0012 Airfoil Models in IRT Test 
Section. 



(b) 21-in- (53.3-cm-) Chord NACA 0012 Model In- 
stalled in IRT Test Section. 

Figure 4.5.1. — concluded. 


aluminum model. The shorter-span models were placed 
between end plates. The models’ mid spans were located 
on the tunnel centerline (midway between floor and ceil- 
ing of the test section). Three to five horizontal lines 
were typically marked around the leading edge of each 
airfoil to indicate the mid-span and of positions 1 and 2 
in above and below mid-span for locating ice-tracing 
templates. Shapes were recorded by hand tracing. 

4.6. IRT Scaling Test Procedures 

The typical procedure involved setting the temperature, 
airspeed and spray-bar air and water pressures for the 
test. Once these pressures stabilized, the water valves 
were opened to initiate the spray, and the spray timer was 
started. Before 1997, once the water valves were 
opened, water and air pressures sometimes took a minute 
or more to stabilize. Because spray duration is timed 
from the opening of the valves, when testing with the 
earlier system it was necessary either to shield the model 
until the spray had stabilized or to add time to correct for 
the ramp-up period. This correction was typically 20 sec. 
Spray times reported for pre-1997 tests were the desired 
exposure times before this spray extension was added. 
For icing exposures of about 10 minutes or more, errors 
in the correction applied to the spray time could be on the 
order of 30 sec without having any noticeable effect on 
the total quantity of ice accreted. A correction to the 
spray time was applied by Olsen and Newton 1 ' in their 
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scaling tests in the 1980’s and by Bilanin and Anderson 42 
and Anderson 33 ’' 4 64 ' 65 for models of 21 -in chord and lar- 
ger in scaling tests before 1996. 

1RT scaling tests in 1995 employed a metal shield in 
front of models smaller than 21 -in chord. When the 
cloud had stabilized, the shield was raised and the spray 
timing initiated. A complete description of this tech- 
nique is given in references 34 and 65. Although the 
spray conditions were steady throughout the test using 
this shield, the test-section blockage changed when the 
shield was raised, and some transients in the airspeed 
were experienced. 

An improved shield system was used for tests in 1996. A 
fabric patch was attached to the front of the model using 
Velcro strips. This shield was large enough to cover the 
leading-edge region of the model back to at least 20% 
chord and 16 in of the span around the centerline of the 
tunnel. Because the shield conformed to the shape of the 
model, it did not change the tunnel blockage as the metal 
shield had done. A rope was attached to the shield and 
passed through a port in the ceiling of the tunnel. When 
the cloud was stabilized, the shield was manually pulled 
free from the model and through the port. 

The IRT spray system was rebuilt in 1997. 58 The 10-bar 
spray system now recirculates water to permit the stabili- 
zation of air and water pressures prior to opening the 
water valve at each nozzle. This feature virtually elimi- 
nates start-up transients. Consequently, for tests made 
after this upgrade neither shielding nor corrections to the 
spray time have been needed. Comparisons of ice shapes 
obtained before and after the new spray-bar system was 
installed showed good repeatability. These comparisons 
tend to validate the methods used to account for spray- 
bar transients with the old system. 

At the completion of the icing spray time, the spray was 
turned off and the tunnel fan was stopped. Personnel 
then entered the test section, and a thin heated plate was 
used to cut horizontal slices into the accreted ice at the 
desired span-wise locations. The ice shape cross- 
sections were traced onto cardboard templates and later 
digitized for computer storage. Finally, the model was 
cleaned and the procedure repeated for the next test. 

5. Description and Evaluation of Scaling 
Methods 

The evaluation of scaling methods involves choosing a 
series of reference test conditions, applying either size 
scaling or test-parameter scaling to develop correspond- 
ing scale conditions and performing both the reference 
and scale tests in an icing wind tunnel. Ice shapes are 
recorded for both reference and scale tests and compared. 


At the present time, quantification of these ice shapes has 
not progressed sufficiently to permit objective, numerical 
comparison of ice shapes, although some approaches 
have been proposed and evaluated. 66 Consequently, it is 
necessary to plot ice shapes and visually match them to 
make subjective judgments of whether two tests have 
produced the same shape. In this section, scaling studies 
will be reviewed and ice shapes compared from a num- 
ber of studies. The purpose of these comparisons is to 
assess the effectiveness of scaling methods and, when 
possible, to deduce how significant each of the scaling 
parameters identified in section 3 is to determining the 
ice shape. This kind of sensitivity analysis is important 
to developing improved scaling methods, but it can also 
lead to a better understanding of the basic physics of ice 
accretion. The insight gained can be used to refine ice- 
accretion prediction codes as well. 

5.1. Basic Approach to Scaling Method Development 

Section 3 gave the definition and derivation of several 
similarity parameters. The matching of scale and refer- 
ence values of each similarity parameter provides one of 
a set of equations that can be solved for the scale test 
conditions. For sea-level tunnels, after the scale model 
size has been chosen, five scaled test conditions (tem- 
perature, airspeed, drop size, cloud liquid-water content 
and exposure time) need to be determined. For altitude 
tunnels, the scale pressure adds a sixth variable. Thus, a 
system of only five or six equations is needed, and a 
choice must be made of those similarity parameters most 
likely to affect ice shape for the conditions of interest. 
For some scaling methods, various test conditions them- 
selves are matched to the reference values so that even 
fewer similarity parameters are used. In some cases, it is 
not possible or practical to match certain parameters; for 
example, in section 3.2, it was noted that Re a and M a 
cannot both be matched if scale model size differs from 
the full scale. 

The scale test is defined first by the model size, which is 
selected by the test designer. For test-condition scaling, 
the model size is equal to the reference. For both test- 
condition and size scaling the scale model geometry and 
AOA must be the same as the reference. Two test pa- 
rameters can be found by applying similarity of drop 
trajectories and water catch. The remaining parameters 
can be computed from the equations formed by matching 
various terms in the energy balance and expressions re- 
lated to surface phenomena. The scaling methods that 
have been proposed differ primarily in which of these 
energy terms or surface parameters are matched. 

Table I compares scaling methods that have been evalu- 
ated by various studies in the past 30 years. Earlier 
methods described by Sibley and Smith, 16 Dodson, 17 and 
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Table I. — Comparison of Some Published Scaling Methods 


Method 

T est Conditions 

Similarity Parameters 

d 

tst 

Pst 

V 

8 

LWC 

T 

K 0 

A c 

n n 

b 

r A 

0 

9 

Ca 

Re a 

We s 

M a 

L WC x time 68 

M 

M 

M 

M 

M 

s 

calc 

(m) 
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* AEDC Method 

M matched to reference value 

(m) coincidentally matched to reference value 

s specified by user 

s 1 one of the test conditions t s „ p st , V, 8 or L WC 
can be selected 


a determined by ambient conditions and test 

parameters 

calc calculated from matched parameters 


Hauger and Engler 67 will not be discussed here. Methods 
proposed in the past few years that involve similarity 
parameters other than those in the table have not been 
included but will be discussed in section 5.3. Ruff 7 
evaluated four approaches to scaling, and all four are 
listed in the table. The “modified-Ruff method” shown 
is an adaptation of his method 4 for sea-level facilities — 
that is, facilities with no ability to control test-section 
pressure. Because this facility limitation exists for the 
IRT as well as many other facilities, the “modified-Ruff 
method” will be simply called the “Ruff method” in the 
remainder of this manual. 

Table 1 includes examples of both test-condition and 
size-scaling methods, indicating the similarity parameters 
that are matched between scale and reference. The gen- 
eral way scale test conditions are determined is shown: 
whether they are matched to the reference, calculated by 
matching scaling parameters or determined by ambient 
conditions. For example, in the ONERA method the user 
first specifies the scale model size, d. Then the test con- 
ditions t st , V, 8, L WC and r are calculated by matching 
the similarity parameters K n , A c , and b. Because the 


ONERA tunnel at Modane, for which this method was 
developed, has no control of pressure or temperature, the 
static values of these conditions are determined by the 
environment. In this way, all the conditions for the scale 
test are established. 

Note that the Ruff method 2 reduces to the Olsen method 
if the scale model size, velocity and drop size are chosen 
to be the same as the respective reference values. Both 
the ONERA and Ruff methods will be described in more 
detail in section 5.3. 

5.2. Methods to Scale Test Conditions 

In section 4.1 it was noted that the calibrated envelope of 
an icing facility can leave some combinations of MVD 
and LWC unobtainable. Facilities are also limited with 
respect to the velocity which can be tested. Furthermore, 
substitutions for temperature or pressure are sometimes 
necessary. By applying test-conditions scaling, alterna- 
tives for each of these test conditions can often be found 
that fall within the facility capability yet provide an ice 
shape that simulates that of the reference (desired) condi- 
tions. Test-condition scaling methods will be discussed 
here. 
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5.2.1. LWC 

5.2.1. 1. LWC x Time = Constant 

The simplest scaling method is used for L WC scaling 68 
and explicitly satisfies only water-catch similarity 
(matching the accumulation parameter.) For this applica- 
tion, scale model size, temperature, pressure, airspeed 
and drop size are matched to the corresponding reference 
values: 


C S - C R 

(5.1) 

^ st,S ~ ^ st,R 

(5.2) 

Pst,S — Pst,R 

(5.3) 

ii 

(5.4) 

S s - S R 

(5.5) 


It is understood, of course, that scale model geometry 
and AOA are the same as the reference. When the scale 
and reference values for the accumulation parameter, 
equation (3.26), are equated and constants cancelled, the 
product of the liquid-water content and time is seen to 
match: 

LWC s z s = LWC r t r (5.6) 

The user chooses the scale LWC and from equation 
(5.6), the scale accretion time can be found. The lead- 
ing-edge heat balance is assumed to be satisfied simply 
by matching the static temperature. The modified inertia 
parameter, K„, is matched because the model size, veloc- 
ity, cloud drop size and ambient pressure have been 
matched. The effect of liquid-water content on the heat 
balance is ignored. 

In addition, because the model size, airspeed and drop 
size are maintained at reference values, Re a , M a and vari- 
ous Weber numbers also match coincidentally. From 
equations (3.56) and (3.57) it is apparent that equating 
the scale and reference static temperatures, pressures and 
velocities insures that scale and reference values of both 
<f i and 9 will also be equal. Thus, this method satisfies 
equality of most of the parameters listed in section 3.8; 
the scale values of only n 0 , b , r A , and hfaJd differ from 
the reference values. 

Evaluations of the ‘ LWC x time = constant’ method have 
been made in the NASA-Glenn IRT 33 ' 34 ’ 64 us i n g both 
cylinders and a 21-in- (53.3-cm-) chord NACA 0012. 


More details about the models and test procedures were 
given in the references. Figure 3.5.2 (a) and (b) showed 
one set of results for which the L WC was decreased from 
1.4 to 1.0 then to 0.8 g/nr’ with the time increased to sat- 
isfy equation (5.6) and all other test conditions constant. 
Although the total quantity of ice accreted appeared to 
be constant, the increasing freezing fraction caused a 
decrease in the horn angle. The leading-edge ice thick- 
ness would also tend to increase with decreasing LWC 
(increasing no). For rime conditions, this method would 
be effective providing the temperature was low enough 
that rime was maintained for all LWC’s. For glaze ice at 
very small accumulation parameters (very short accretion 
times), the change in horn angle might not be noticed. 
But in general, figure 3.5.2 shows that this method fails 
because it does not maintain freezing fraction constant. 

The ‘LWC x time = constant’ scaling method is simple to 
apply and convenient to use when data are needed at 
liquid-water contents not available in flight or not 
achievable in an icing tunnel. Unfortunately, as figure 
3.5.2 demonstrates, simply matching static temperatures 
is not an adequate way to account for the leading-edge 
energy balance. 

5.2. 1.2. Olsen Method 

A refinement of the ‘L WC x r = constant’ method was 
suggested by Olsen and Newton. 1 ’’ 68 As with the previ- 
ous method, the user selects the scale L WC, while main- 
taining chord, velocity and drop size unchanged (i.e., 
equations (5.1), (5.4), and (5.5) are satisfied.) Since 
model size, airspeed and water drop size are unchanged 
from the reference values, K 0 ,s = Knjt- The scale icing 
time is once more found from equation (5.6) to insure 
that A C: s = A, :R . But, in the Olsen method, the freezing 
fraction is matched instead of the static temperature. The 
scale static temperature is then calculated from the scale 
freezing fraction. This constraint of calculating a static 
temperature to insure a specific freezing fraction makes 
the Olsen method less convenient to use that the ‘LWC x 
r = constant’ method, but gives more confidence that the 
scale ice shape will simulate the shape which would have 
accreted if the desired conditions could have been tested. 

Olsen and Newton proposed the Olsen method for use in 
size scaling as well as test-condition scaling. Flowever, 
we have subsequently learned that when size is scaled, 
the scale and reference velocity cannot match if surface- 
water phenomena are properly simulated. Therefore, the 
Olsen method can only be applied when the model size 
of scale and reference cases is the same. 

Experimental evaluations of this method have been re- 
ported in references 33 and 34. The scale results were 
compared with the same reference cases as used in the 
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NASA Glenn ‘ LWC x time = constant’ study discussed 
above. Typical results were shown in figure 3.5.3. It is 
apparent from that figure that the scale ice shapes using 
the Olsen method are significantly better simulations of 
the reference shapes than were those from the ‘ LWC x 
time = constant’ method (fig. 3.5.2). The Olsen method 
results in the correct size, location and angle of glaze 
horns, and the feathers aft of the horns also appear to be 
well matched in size to those on the reference accretion. 

Ruffs 7 method 2 was identical to the Olsen method. 
However, in his tests he found that ‘‘‘...the size and loca- 
tions of the glaze horns are not adequately reproduced.” 
He concluded that this method can only be used to simu- 
late the type of ice but not the shape. It is possible that 
Ruff s tests failed to give the good match of shapes dem- 
onstrated in the IRT because of problems with his facility 
L WC calibration. This speculation is also consistent with 
Bilanin’s 39 inability to calculate experimental freezing 
fractions from Ruff s ice shapes that were consistent with 
analytical values (see section 3.5.14). 

The success of the Olsen method over the ‘ LWC x r = 
constant’ method in simulating the reference ice-shape 
features is a result of requiring that the scale freezing 
fraction match the reference. Glaze ice shapes appear to 
be more sensitive to changes in n„ than to variations in 
the other heat -balance parameters, b , <j>, 6 or r A . The Ol- 
sen-method results show once again that the freezing 
fraction is one of the fundamental similarity parameters 
that must be matched for good scaling. The Olsen 
method is based on a fundamentally sound approach, but 
has been experimentally validated over only a specific 
range of test conditions. Therefore, some caution may be 
warranted for application to other conditions. 

5.2.2. Temperature 

Figures 3.5.4 and 3.5.5 showed that when the total tem- 
perature is higher than about 27°F ice shape characteris- 
tics can be quite different from those at lower tempera- 
tures. The size-scaling method recommended below in 
section 5.3.7 will give scale total temperatures higher 
than the reference, and it is possible for some sets of 
conditions for the resulting scale total to exceed the 27°F 
limit. If the original reference total temperature was be- 
low 27°F, the scale temperature should be scaled to a 
value below this limit as well. 

With model, velocity and drop size fixed, temperature 
scaling can be accomplished using the Olsen method 
described in section 5.2. 1.2. The premise of the Olsen 
method is that with c, V and S constant, scale and refer- 
ence ice shapes will match if the scale no matches the 
reference. To apply the Olsen method for temperature 
scaling, one first selects a scale static temperature so that 


the scale total temperature will be below 27°F. Then the 
scale L WC is calculated so that scale freezing fraction is 
the same as the reference value. Finally, the accretion 
time is determined by matching scale and reference A c . 

5.2.3. Drop Size 

Figures 3.3.6 and 3.3.7 showed that drop size variations 
from 15 to 20, from 25 to 40 and from 20 to 55 pm had 
no apparent effect on the main ice shape for the condi- 
tions and models tested. These observations imply that 
for Appendix-C conditions, if a test is needed with a 
value of § outside the tunnel capability, it should be pos- 
sible to substitute a value within the tunnel operating 
map to simulate the desired ice shape. With model size 
and velocity unchanged from the reference, temperature 
would have to be adjusted to maintain n 0 because S influ- 
ences n 0 through J3 0 . Time would also need to be set so 
that the product /3qA c was maintained the same as for the 
reference conditions. 

There are some limitations to this approach, however. 
First, very small drops cannot be used to simulate accre- 
tions with larger drops. The discussion following equa- 
tion (3.8) noted that when the inertia parameter, K, is less 
than 1/8, drops will not impinge on the model. There- 
fore, there must be a lower limit of drop MVD, depend- 
ent on model size and velocity, below which ice shapes 
will differ from those produced with larger drops. Fur- 
thermore, the lack of effect of MVD on ice shape has not 
been tested for drop sizes below 15 pm, and it cannot be 
assumed that the conclusions of figures 3.3.6 and 3.3.7 
apply to drops smaller than 1 5 pm. 

Second, changes to 5 also affect the impingement limit. 
The impingement limit is dependent on j3 0 ', thus, it should 
be possible to establish that limit experimentally with 
any combination of velocity and drop size that produced 
the reference flo for the given model. For any test, it is 
always desirable to determine test conditions such that 
the p n is as close a match to that for the conditions to be 
simulated as possible. 

The insensitivity of ice shapes to drop size as shown in 
figures 3.3.6 and 3.3.7 were for Appendix-C conditions 
and cannot be expected to be valid in general for SLD 
conditions. Even within Appendix C the observation that 
drop size has little effect on ice shape comes from only 
limited testing. Tests over a range of speed, model size 
and cloud conditions are needed before a method of 
drop-size substitution can be recommended with any 
confidence. Some such tests are currently being con- 
ducted at the NASA Glenn IRT for both Appendix C and 
SLD conditions to try to define the limits of conditions 
for which an available MVD can be used in an icing test 
to simulate results with another. 
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5.2.4. Airspeed 

If a test is required to simulate airspeed outside the avail- 
able range for a facility, the user would start by selecting 
a more suitable value. Typically, the desired (reference) 
velocity would be higher than the facility permits, so the 
scale velocity will be lower than the reference. No stud- 
ies have been made specifically to evaluate ways to scale 
airspeed. However, based on knowledge of what similar- 
ity parameters are most important to ice accretion, it is 
possible to propose a method. 

At the present time, Wei appears to be one of the leading 
candidates of similarity parameter related to surface- 
water effects. In order to match scale and reference Wei, 
equation (3.71) shows that scale model size must be 
greater than the reference if Vs <V R . This fact will limit 
the usefulness of this method. For situations for which 
this approach can be applied, however, LWC S can be 
chosen arbitrarily. Then by setting K 0 ,s = K 0fi , A CtS = A CyJi 
and n 0 s = rioji, the scale 8, rand t sl can be found, respec- 
tively. Additional testing is needed to verify this 
approach, so it should be used with caution. 

For very short accretion times, for which features like 
horns have not yet developed, it should be possible to 
ignore We R , match the model size and choose the scale 
velocity arbitrarily, but the limits of this approach have 
not been defined. 

5.2.5. Pressure (Altitude) 

Many icing wind tunnels do not have the capability to 
simulate altitude by controlling test-section pressure. 
The scaling method evaluations shown in this manual 
were obtained primarily in the IRT, which is a sea-level 
facility. Yet, aircraft are flown in icing conditions at 
altitudes of 20,000 ft (6000 m) and more, and it is often 
necessary to simulate the altitude ice accretion in sea- 
level facilities. Thus, the question of what effect pres- 
sure has on ice shapes is an important one and deserves 
review in some detail. 

This issue has been discussed by Bartlett. 43,44 ' 69 In refer- 
ence 44 he reported, “The data base at AEDC indicates 
that pressure changes have an almost insignificant effect 
upon the formation of ice.” He also presented 69 limited 
ice shape data taken on a 2-in-diameter cylinder in the 
AEDC Icing Research Test Cell. The data were for tem- 
peratures of 4 and 23°F, a velocity of 250 ft/sec, MVD of 
28 pm, an LWC of 0.8 g/m J and spray times of 18 and 
19 min. The static pressures were 8, 10, 12 and 14 psia 
for each of the two sets of tests (one set at each tempera- 
ture). Within the tunnel ice-shape repeatability, there 
was no effect of pressure. Bartlett concluded, “The 
effects of pressure (altitude) on icing similitude in the 
range studied herein are negligible. . .” 


The National Research Council, Canada and the Centro 
Italiano Ricerche Aerospaziali (CIRA) jointly conducted 
a series of tests at the NRC Altitude Icing Wind Tunnel. 
Oleskiw, et al. 45 reported that when test conditions were 
repeated with only pressure varied, pressure produced 
“only relatively small changes” in the main ice shape. 
The published ice shapes showed differences no greater 
than those typically observed when tests were repeated. 

De Gregorio, et al 70,71 reported more of the results from 
the same NRC/CIRA test series. Within ice-shape re- 
peatability, pressures of 14.7, 12.7 and 8.2 psia (101, 87 
and 56 kPa) had no effect on ice shape at a temperature 
of 14°F (-10°C), V of 212 mph (95 m/s), 8 of 23 pm, 
LWC of 0.85 g/m 3 and spray time of 10 min. All but one 
of the other sets of data comparing accretions at high and 
low pressure showed evidence of drop freeze out for the 
low-pressure test, so no conclusions about the effect of 
pressure could be reached. 

These experimental results lead to the conclusion that to 
scale altitude test conditions for a sea-level tunnel test, 
the correct ice shapes can be achieved by simply match- 
ing the desired model size and temperature, velocity, 
drop size, liquid-water content and spray time while 
using the available tunnel pressure. 

The observed lack of a pressure effect on ice shapes was 
also cited in section 3.6 to help deduce the importance of 
some of the scaling parameters with regard to their role 
in ice-accretion physics. We can deduce that parameters 
that are strongly affected by pressure cannot have a ma- 
jor effect on ice shape, while those parameters that most 
strongly influence the ice shape must be independent of 
pressure or have only a weak dependence. 

Table II is a compilation of the calculated effects of pres- 
sure on each of the similarity parameters. The pressure- 
altitude relation is from the U.S. Standard Atmosphere. 72 
The values were based on an NACA 0012 with a chord 
of 21 in. The parameters K n , fin, A„ n n , </>, We& We L , Ca 
and r A do not vary with pressure, or the change is less 
than 10% over the pressure range shown. The lack of an 
experimental effect of pressure on ice shape is thus con- 
sistent with what has already been assumed or shown 
about the importance of K 0 (J3 0 ), A c and n 0 on ice shape. 
Even for a chord of 72 in, the effect of pressure on fi lh 
over the range of pressure shown and for the velocity 
given, is less than 10%. For higher velocities, the varia- 
tion of p n with pressure is even less than that shown. The 
matching of K 0 between scale and reference cases insures 
matching of /? 0 , but K 0 in Table II increased by 30% as 
pressure dropped over the range given; thus, while size- 
scaling methods commonly include this parameter, its 
change with pressure suggests that it does not need to be 
matched rigorously. Although drop MVD (and therefore 
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Table II. 

Effect of Pressure (Altitude) on Similarity Parameters 
c, 21 in; t s „ 20°F; t lo „ 24°F; V, 150 mph; 8, 30pm; LWC, 
1 g/m 3 ; time, 7.3 min. NACA 0012 Airfoil. Pressures 
from U.S. Standard Atmosphere, 1976. 12 


Pst, psia 

14.7 

12.2 

10.1 

8.3 

6.8 

Altitude, ft 

0 

5000 

10000 

^5000 

20000 

K 0 

3.58 

3.84 

4.13 

4.43 

4.78 

Po, % 

80.0 

80.8 

81.8 

82.7 

83.6 

A c _ 

1.90 

1.90 

1.90 

1.90 

1.90 

n 0 

0.27 

0.27 

0.27 

0.27 

0.27 

b 

0.54 

0.60 

0.66 

0.74 

0.83 

& °F 

11.0 

11.0 

11.0 

11.0 

11.0 

a °f 

14.8 

16.1 

17.9 

20.0 

22.8 

Ai 

0.24 

0.24 

0.25 

0.25 

0.26 

Re a , JO 4 

8.95 

7.45 

6.16 

5.05 

4.11 

Wes, 10 3 

2.07 

2.07 

2.07 

2.07 

2.07 

We c , 10 3 

1.54 

1.28 

1.06 

0.87 

0.71 

We L , 10 6 

1.16 

1.16 

1.16 

1.16 

1.16 

M a 

0.20 

0.20 

0.20 

0.20 

0.20 

Ca 

1.84 

1.84 

1.84 

1.84 

1.84 


K 0 or Pd) are important for determining impingement 
limits, figures 3.3.6 and 3.3.7 demonstrated that it does 
not have a significant effect on the main ice shape. The 
studies of the effects of pressure cited above concen- 
trated on main ice shape effects. 

Both 0 and h change significantly with pressure; there- 
fore, neither of these parameters can have a strong inde- 
pendent effect on ice shape. 

Since pressure has little or no discernible effect on ice 
shape, if the Weber number does affect ice shape, it must 
be in a form that is independent of pressure. Because of 
the change in air density, We c (based on model size and 
air density) decreases to about half its sea-level value 
when pressure is reduced to simulate a 6000-m altitude. 
Thus, it would seem that the appropriate Weber number 
must be one dependent on the density of water, rather 
than air; for example, Wegor We L (eqs. (3.70) or (3.69)). 

The effect of pressure on air density also reduces Re a to 
about half its sea-level value when pressure is reduced to 
the 20,000-ft-altitude level. This result suggests that if 
Re a has any effect on ice shape, it is likely to be a weak 
one. 

Thus, all the evidence indicates that pressure effects can 
be ignored when altitude ice accretion is to be simulated 
in sea-level tunnels. However, a rigorous pressure- 
scaling method can be devised based on matching those 


similarity parameters that have been shown to be impor- 
tant. This scaling method has not been tested, and it 
needs to be validated in a pressure-controlled tunnel be- 
fore it can be recommended. With model size, velocity 
and L WC matched and the reference and scale pressures 
defined, the scale drop size is calculated so that the scale 
and reference K n match. The resulting scale drop size 
could differ from the reference by 15 to 20%. Because 
drop size has a minor effect on freezing fraction, a new 
temperature should be calculated such that the scale 
freezing fraction matches the reference. The scale 
temperature will differ from the reference by only a small 
amount, perhaps on the order of 1°F. Because d, V and 
L WC were all matched to the reference, the accumulation 
parameter will match, and the scale and reference icing 
time will be the same. 

5.2.6 Summary of Methods to Scale Test Parameters 

Table III summarizes the methods discussed above for 
scaling test parameters. The Olsen method is recom- 
mended for scaling either LWC or t sl , using the 
approaches described in sections 5. 2. 1.2 and 5.2.2 and 
outlined in Table III. Methods to scale drop size, air- 
speed and pressure are proposed, based on matching the 
similarity parameters of importance to scaling ice shapes. 
These proposed methods for scaling 8, V and p s , have not 
been tested, and thorough testing is required before those 
methods can be recommended. 

5.3. Methods to Scale Model Size 

A common application of scaling is for testing situations 
in which the facility restricts the size of model that can 
be used. Scaling methods are applied to find test condi- 
tions for the subscale model that will produce ice shapes 
similar to those that would accrete on a full-sized, or 
reference, model. A number of methods have been 
developed through the years but only three will be dis- 
cussed here to illustrate some of the different approaches 
which have been used. The recommended method to 
scale size will be given in section 5.3.7. 

Table IV compares the way the scaled test parameters are 
determined in the three scaling methods. For all three, 
with scale model geometry and AOA the same as the 
reference, the user selects the scale model size. The test 
conditions are then found by different strategies. In the 
following sections, these approaches to finding the 
appropriate scale conditions will be discussed further. 

5.3.1. Ingelman-Sundberg 

The Swedish-Soviet Working Group on Aircraft Safety 
described a size-scaling method 73 that relied on matching 
only the similarity parameters K 0 and A, and test condi- 
tions t s , and LWC (see Table I.) The user specifies the 
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Table 111. — Methods for Scaling Test Conditions 


Test Condition 
to be Scaled 

T est Conditions 

Similarity Parameters 

d 

1st 

Pst 

V 

5 

LWC 

T 

K 0 

A c 

n 0 

b 

<t> 

e 

Re a We s We L M a 

LWC* 

M 

calc 

M 

M 

M 

s 

calc 

(m) 

M 

M 




(m) (m) (m) (m) 

t st * 

M 

s 

M 

M 

M 

calc 

calc 

(m) 

M 

M 




(m) (m) 

5 

M 

calc 

a 

calc 

s 

M 

calc 

M 

M 

M 





V 

calc 

calc 

a 

s 

calc 

M 

calc 

M 

M 

M 




M 

Pst 

M 

calc 

s 

M 

calc 

M 

M 

M 

(m) 

M 




(m) 


* Olsen method 


M matched to reference value 

(m) coincidentally matched to reference value 

s specified by user 


a determined by ambient conditions and test 

parameters 

calc calculated from matched parameters 


scale model size and velocity. With model size and ve- 
locity given, the drop size can be determined by match- 
ing scale and reference K n , and, with LWC S and V s 
known, icing time can be found by matching A c . If scale 
and reference velocities are the same, the freezing 
fractions will also match, but, other than for this special 
case, this method does not recognize the importance of 
the freezing fraction. Furthermore, as we now know, 
velocity has too great an effect on ice shape to be chosen 
arbitrarily. Tests in reference 64 showed that this 
method is inadequate in general to provide a good simu- 
lation of glaze ice shapes. However, for rime ice or for 
very short accretion times, for which glaze features have 
not yet formed, the Ingelman-Sundberg method should 
be adequate. 


Table IV. 

Methods to Determine Scale Test Conditions 
For Size Scaling 


Test 

Parameter 

Ingelman- 

Sundberg 7 ' 

ONERA 14 ' 15 

Ruff 7 

c s 

user selects 

user selects 

user selects 

tst.S 

tst.R 

t lolt s known 

II 

Co 

V s 

user selects 

no,s = no.R 

user selects 
(typical) 

MVD S 

Ko.s = Ko.r 

Ko.s = Ko.r 

Ko.s = Ko.r 

LWC S 

= LWC r 

b s =b R 

nojs = no. R 

T S 

II 

Co 

o 

A c ,s = A Cy R 

A c ,s = A c ,r 

Pst,S 

Ptot,S Ptot,R 

Ptots known 

II 


5.3.2. ONERA 

Charpin and coworkers 14 15 of ONERA described a scal- 
ing technique in which the parameters K n and A c and two 
terms from the energy balance, the freezing fraction, no, 
and the relative heat factor, b, were matched to the refer- 
ence values. If the surface temperature is assumed to be 
at the freezing point of water (t s = tf= 0°C), and the 
properties of air and water are substituted into equation 
(3.54), the resulting equation is ONERA’s form of the 
energy equation: 


1.058x10 


- = /* (1 + 6) + (1732 K)i^ 

Pst 


cal V 2 

+79.7 — n„b + (3.6458 +b) 

g 8373- 


s K 


(5.7) 


To develop equation (5.7), the following properties of 
air and water at 0°C have been used: ha/h c = 2.9 gK/cal, 
A v = 597.3 caEg,/),™ = 610.8 nt/m 2 , A f = 79.7 caEg, c PiWS 
= 1 caEgK, r = .875, and c p = .24 caEgK. Note that the K 
appearing in equation (5.7) is the absolute temperature 
unit, Kelvin, not the inertia parameter. 

The procedure to find scale test conditions is as follows: 
Knowing the reference test conditions, b K can be calcu- 
lated from equation (3.55), then equation (5.7) is solved 
for n IIR . b s and n 0 ,s can now be set equal to these respec- 
tive reference values. Substituting these into equation 
(5.7) gives the energy equation for the scale case, which 
now relates four of the scale test conditions: velocity, 
static temperature, static pressure, and vapor pressure. 
The Modane tunnel, for which the ONERA method was 
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developed, does not permit control of either static tem- 
perature or static pressure. The total temperature and 
pressure in the settling chamber can be measured (they 
are, of course, the ambient temperature and pressure at 
the time of the test.) The static temperature in the test 
section is a function of velocity and the known total 
temperature. Likewise, static pressure is a function of 
velocity and total pressure. Finally, the vapor pressure in 
the test section is a function only of static temperature. 

Thus, including equation (5.7), we have a system of 
four equations that involve only the unknowns, V s , t slA 
p st s and p w- s. These equations can be solved by using a 
simple iterative approach. For example, one could start 
by assuming a value for V s , perhaps setting V s = V R as a 
first guess. Then t st ,s and p uS can be calculated, and, 
knowing t st ,s, p w ,s can be determined. Finally, equation 
(5.7) can be solved for a new estimate of V s . The proc- 
ess is repeated until the test conditions found in successive 
iterations match. With V s , t stS and p stS thus established, 
the scale drop size can be found from equation (3.18). 
Finally, LWC S is calculated by matching the expressions 
for b s and b R , using the definition in equation (3.55), 
along with definitions of f} 0 from equation (3.13) and h c 
from equation (3.32). 

As discussed following equation (3.32), Charpin and 
coworkers used Nu a oc Re a 0 ' 8 to determine h c rather than 
Nu a oc Re a 05 (eq. (3.33)). The 0.8 exponent is not consis- 
tent with measured heat-transfer coefficients in the 1RT, 
but it may reflect the situation in the Modane tunnel. It is 
not recommended for other facilities without further 
study of heat transfer coefficients in those tunnels. 

The ONERA method was evaluated in the Modane tun- 
nel with scales of 1/6 and 1/12, and it was reported that 
this approach was effective. 141 ~ Tests of a modified form 
of this method were reported in references 33 and 64, but 
rather than calculating the scale velocity as specified in 
references 14 and 15, a value was chosen arbitrarily. 
With b s = b R and n os = n 0R , the scale temperature was 
then determined from the energy balance, equation (5.7). 
This approach did not adequately match scale and ref- 
erence ice shape features for size-scaled tests. Tests 
reported in reference 33 using a convective heat-transfer 
film coefficient dependency of Re a °' 8 instead of Re a 0 ' 8 
seemed to improve the ability of the method to scale the 
ice shapes somewhat, but features still did not compare 
as well as one would want. None of the tests of refer- 
ences 33 and 64, however, applied the ONERA method 
as used in the Modane tunnel; therefore, these assess- 
ments only apply to a method for which the same 
similarity parameters were matched. 


The ONERA method typically results in scale velocities 
less than the reference. We now know that scale veloci- 
ties need to be higher than the reference for good simula- 
tion of reference ice shapes using sub-scale models. 
Thus, the ONERA method would almost certainly be 
improved by replacing b with We L in the list of four simi- 
larity parameters to be matched, then modifying the 
procedure as follows: First find the scale velocity by 
matching the We L , then static temperature and pressure 
can be determined from the totals, vapor pressure can be 
calculated from static temperature, and the energy equa- 
tion solved for bs with scale and reference freezing 
fraction matched. Next, determine scale LWC from b S - 
The remaining test conditions would be found as they are 
in the present ONERA method. Such an approach needs 
to be evaluated before it can be recommended, however. 
To date, no such evaluation has been performed. 

5.3.3. Ruff(AEDC) 

Ruff 7 was the first to systematically evaluate scaling 
using various combinations of the similarity parameters, 
K n , A,, no, b, (/> and 9. Fie tested four methods (see Table 
I.) The first matched only K» and A c , the second K n , A c 
and no, the third, K 0 , A c , n 0 and b, and finally, the fourth 
matched K n , A c , n 0 , ^and 9. Note that Ruffs third 
method involved the same similarity parameters as the 
ONERA method, but unlike that approach, Ruff assigned 
a scale velocity and calculated a scale static temperature 
by matching freezing fraction. Ruff also used equation 
(3.33) for the Nusselt number. 

Ruff achieved the best agreement of scale and reference 
ice shapes when K 0 , A c , n 0 , (p and 9 were matched. It can 
be seen from equation (3.59) that when these last three 
parameters are matched, so too is b. This approach has 
been widely adopted, and is often referred to as the 
AEDC Method. It applies to tunnels with altitude- 
simulation capability, because it permits determining the 
value of the static pressure. The matching of five simi- 
larity parameters gives five equations that are solved for 
five of the six scale test conditions t s „ p st , V, 5, L WC and 
r (see Table IV). This method permits the user to 
arbitrarily assign a value to one of the test conditions. 
Typically, V s would be selected, and this is the way Ruff 
conducted his tests. For facilities with limited airspeed 
capability, setting a value for scale velocity seems rea- 
sonable, even necessary; however, it totally ignores the 
physics behind the large effect that velocity has on ice 
shape (see figure 3.6.1). The effect of velocity, demon- 
strated in references 18, 35 and 46, was not appreciated 
until the late 1990’s. 

In Ruff s method 4 with model size selected, suppose a 
value for V s is arbitrarily picked. Then t stS and p,, s can 
be determined by solving the two equations formed by 
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matching <f> s and 6 S to their respective reference values. 
K 0 ,s = K 0 r can now be solved for S s , and L WC S found by 
matching n 0 ,s with n 0R . Finally, z s is established by 
matching scale and reference A c . 

For tunnels that cannot control the test-section pressure, a 
modified Ruff Method has been used (refs. 33 and 34) in 
which, typically, the air-energy transfer parameter is 
ignored and a value for scale velocity can again be 
selected. For a tunnel like the IRT, for which the settling 
chamber upstream of the test section is vented to atmos- 
phere, the total pressure is the ambient pressure. Then, 
the scale test-section static pressure is determined from 
this total pressure and the scale velocity. The static tem- 
perature is found from matching </> s to <f fo, as shown in 
Table IV. The scale drop size is established by matching 
scale and reference K n , scale L WC by matching no and 
scale icing time by matching A c . 

A variation of the modified Ruff method 46 for sea-level 
tunnels is to determine scale temperature by matching 0 
while ignoring </>. This approach results typically in 
lower temperatures than by using constant (j> . Reference 
46 concluded from limited testing that it made no differ- 
ence to the ice shape whether temperature was found by 
matching <f> or 0. Reference 49 also showed scale ice 
shapes that matched the reference when temperature was 
found by applying constant 0 with (f) unmatched. These 
results are not surprising since neither (/) nor 0 has been 
shown to have an effect on ice shape independent of the 
freezing fraction. Thus, it is highly probable that the 
scale temperature (or, alternatively, L WC) can be chosen 
arbitrarily with the L WC (or temperature) determined by 
matching scale and reference no. Anderson and Tsao 74 
demonstrated good matches of scale and reference 
shapes for conditions for which neither (j) nor 0 matched 
reference values. Additional testing is needed to confirm 
this result, but if validated, it would provide useful flexi- 
bility in the definition of scale test conditions. 

5.3.4. Methods to Choose Scale Velocity 

The Ingelman-Sundberg size-scaling method discussed 
in section 5.3.1 left the choice of scale velocity to the 
user, as though any velocity would be equally appropri- 
ate. The Ruff method, too, permits the user to arbitrarily 
select a scale velocity. The ONERA method calculates 
V s by matching freezing fraction; this approach gives 
scale velocities that are typically lower than the refer- 
ence. In this section we shall look at alternate ways to 
select scale velocity and see how their choice affects ice 
shape. A recommended method will be given. 

The simplest way to select the scale velocity is to equate 
it to the reference value, but when high airspeeds need to 
be simulated, facility limitations may make it tempting to 


use a scale velocity that is less than the reference. 
Matching velocities is also equivalent to approximately 
equating the scale and reference Mach numbers, so if M a 
could be shown to be important in the physics of icing, 
this is a valid approach. However, most icing encounters 
take place at low Mach numbers where compressibility is 
not important, and reference 35, in limited testing, found 
no effect of Mach number up to M a = 0.5. 

If M a does not have to match and, as we have seen, 
velocity has a significant effect on ice shape independent 
of the parameters K 0 , A c , n 0 , b, <j> and 6, it is important to 
determine the velocity that gives the most faithful simu- 
lation of reference ice shapes. Tests to evaluate a variety 
of ways to choose scale velocities were reported in refer- 
ences 35 and 46. Figure 5.3.1 reproduces results from 
reference 35 at a nominal freezing fraction of 0.5. The 
figure shows how scale ice shapes vary with different 
scale velocities. In each part of the figure the shaded 
shape is the reference, recorded on a 21 -in-chord NACA 
0012 airfoil with a velocity of 150 mph. The solid line is 
the 10.5-in (26.7-cm) scale result. The changes in horn 
angle as the scale velocity increased is entirely consistent 
with what was shown earlier in figure 3.6.1. 

The scale velocity was calculated using different meth- 
ods for each portion of the figure. Figure 5.3.1 (a) gives 
the scale ice shape for a velocity of 127 mph (57 m/s), 
just 81% of the reference value. This velocity resulted 
from matching the scale and reference film thickness 
from equation (3.73). While the reference horns pro- 
jected slightly forward, the scale horns were swept back, 
comparing poorly with the reference. Note that this 
correlation for hf,i m was based on heavy-rain studies 
rather than Appendix-C icing conditions. A later study 
by Feo produced equation (3.74). When this latter water- 
film thickness was matched between scale and reference, 
the resulting V s was significantly higher than V R , and the 
resulting scale ice shapes were a better match to the ref- 
erence 49 than that shown in figure 5.3.1 (a). Results from 
tests using this similarity parameter will be discussed 
later in this secton. 

When the scale velocity was increased to match the ref- 
erence, 150 mph (67 m/s), the scale shape in figure 5.3.1 
(b) resulted. While the horns moved forward slightly, the 
change was not sufficient to provide a good match of the 
reference shape. 

For figure 5.3.1 (c) the scale velocity was increased to 
195 mph (87 m/s), or 30% higher than the reference. 
This velocity was that obtained by matching the scale 
Wes to the reference value. Although the scale horn an- 
gle did not precisely match the reference, the shapes 
probably agreed within typical repeatability. 
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Figure 5.3.1. — Effect of Scale Velocity on Scaling of a Reference Ice Shape 35 . Reference Test, 3-26-98 Run 4. hjuJd 
From Equation (3.73); Values for hjuJd are Relative to Reference. 
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Test results for one more incremental increase in scale 
velocity are shown figure 5.3.1 (d). This velocity was 
the average of that determined by matching Wes and that 
from matching Re a . For the tests reported in reference 
35, the best scaling resulted using this average -velocity 
approach. 

Although not shown in figure 5.3.1, tests reported in ref- 
erence 35 were also made with a scale velocity found by 
matching scale and reference Re a . This scale velocity 
was 2.1 times the reference. The angle between the 
horns for the scale shape was more acute than that of the 
reference shape; therefore, this scale velocity appeared to 
be higher than necessary for good simulation of the ref- 
erence shape. 

The results shown in figure 5.3.1 are for a freezing frac- 
tion of about 0.5. Tests in references 35 and 46 were for 
either 'A- or 2/3 -size scaling and were conducted at addi- 
tional nominal freezing fractions of 0.3, 0.6 and 0.8; the 
general conclusion of all these tests was that the average- 
ly method of choosing scale velocity gave the best match 
of scale and reference ice shapes. This success of the 
average - V method was interpreted to be an indication 
that both Wes and Re„ are important to the glaze icing 
process, and this conclusion is consistent with the belief 
that water-film phenomena play some part in determin- 
ing ice shapes. 


Feo 47 proposed matching the water-film thickness of 
equation (3.74) to find the scale velocity. This thickness 
expression included terms only for LWC and Re a , 
omitting the Weber number. Thus, it seemed to ignore 
some of the physical dependencies known from earlier 
tests to be important (see figure 3.6.1 and associated dis- 
cussion). Nevertheless, 1/2-size scaling tests were con- 
ducted 49 for which the film thickness of equation (3.74) 
was included as a similarity parameter to supplement the 
Ruff method. Freezing fractions of 0.3, 0.4, 0.5, 0.8 and 
1 .0 were tested. The resulting scale velocities were close 
to those of the average - V method; thus, it was not sur- 
prising that scale ice shapes were again a good match for 
the reference. Water-film parameters other than hjuJd 
have also been used 49 ’ 51 to find scale velocity. These 
include Weber numbers based on water-film thickness 
using either free-stream or water-film-edge velocities. 

Additional examples of scaling with different scale ve- 
locities comes from IRT tests by Chen 18 with GLC 305 
airfoils. Figure 5.3.2 gives two sets of comparisons for a 
freezing fraction of 0.3. The reference shape for each 
pair is shown shaded and the scale is represented with a 
solid line. The reference models had a chord of 36 in and 
the scale, 18 in. 

The scale test in figure 5.3.2 (a) used a scale velocity that 
permitted a near match of the Weber number based on 
drop size, Wes The resulting scale shape agreed with the 
reference within normal ice-shape repeatability. Even 



x/c 


x/c 


c = 36 in. Run 736.11 
c = 18 in. Run 718.31 


c = 36 in. Run 736.21 
c. = 18 in, Run 718.21 


(a) 

Constant We; Wes/We R = 

1.04. 



(b) 

Nearly Constant We L ; We LS /We LR = 1.12. 


Run 

c, 


V, 

5, 

LWC, 

T, 

flo, 

A 

n 0 

h 

<f. 

e. 

Re a , 

We s 

We a 

We h 

in 

°F 

mph 

pm 

g/m 3 

min 

% 

4i c 

u 

°F 

°F 

10 4 

10 3 

10 3 

10 6 

(a) 736.11 

36 

16 

148 

50 

1.07 

13.8 

85.4 

2.61 

0.30 

0.75 

14.8 

20.6 

12.4 

3.35 

2.10 

1.64 

718.31 

18 

14 

198 

29 

1.34 

4.0 

85.5 

2.55 

0.30 

0.78 

16.3 

20.5 

8.2 

3.49 

1.86 

1.48 

(b) 736.21 

36 

16 

158 

50 

1.01 

13.5 

85.9 

2.58 

0.30 

0.74 

14.7 

20.1 

13.3 

3.84 

2.40 

1.88 

718.21 

18 

12 

237 

27 

1.35 

3.3 

85.6 

2.53 

0.29 

0.86 

17.8 

20.5 

9.7 

4.65 

2.61 

2.11 


Figure 5.3.2. — Scaling Results with Alternate Scale Velocities. GLC 305 Airfoil at 0°AOA. Ice-Shape Data from 
March 1998 IRT Tests by Chen. 18 
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better agreement of shapes was evident in figure 5.3.2 
(b), for which a higher scale velocity was tested. For all 
the tests represented in figure 5.3.2, scale and reference 
We L matched within about ±10%. 

Limited evaluation of scaling with matched We L using 
NACA 0012 airfoils has been published 74 recently. 
These tests were made using SLD reference conditions. 
Figure 5.3.3 gives results for freezing fractions of 0.3 
(fig. 5.3.3 (a)) and 0.5 (fig. 5.3.3 (b)) for a scale-to- 
reference size ratio of 1:1.7. Excellent agreement be- 
tween scale and reference shapes is evident. The ex- 
perimental results of figures 5.3.2 and 5.3.3 are consis- 
tent with the evidence in section 3.6 that suggested that 
We L is an important similarity parameter. Flowever, as 
we have seen, other methods discussed in section 3.6 to 
find scale velocity have also been used with good ex- 
perimental results. 

The reason that a number of methods of finding scale 
velocity give similar scaling results is that the resulting 
scale velocities are not significantly different. Table V 
shows sample scale velocities found by various methods 
for a reference velocity of 150 mph and for scale-to- 
reference size ratios of 1/2 and 1/3. For this reference 
velocity, for the 1/3 -scale case, no solution could be 


found for some of the methods; these situations are 
indicated with asterisks. For half-size scaling, several 
methods give scale velocities in the range between that 
given by constant Wes and that from the average- V ap- 
proach. Within this range of scale velocities, it is diffi- 
cult to discern significant differences in ice shape. Thus 
icing tests cannot help us identify the most important 
water- film similarity parameter. 

Figure 5.3.4 graphs some of the information in Table V 
along with additional data. Flere, the scale-to-reference 
velocity ratio is given as a function of the reference -to- 
scale size ratio for a range of size ratios from 1.7 to 4. 
The shaded band represents scale velocity ratios in the 
range of ±15%. The velocity ratio for most of the meth- 
ods was calculated for a reference velocity of 150 mph. 
At this reference velocity, however, the routines used to 
calculate scale velocity failed for the average-velocity 
and constant- iffy methods at scale ratios greater than 2; 
therefore, to obtain values over the complete range of 
scale ratios, the reference velocity used for these methods 
was 115 mph. It is evident that even with a scale ratio as 
high as 4 there are a number of methods whose scale 
velocities differ by no more than ±15%. This is not a 
sufficient difference to produce significant variations in 
ice shape. 
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Figure 5.3.3. Examples of Scaling for which Scale Velocity was Found by Matching We/. NACA 0012 Airfoils at 
0° AO A. 74 
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Table V. 

Scale Velocity for Different Methods. 

NACA 0012 Airfoil; t slR , 6°F; S R , 40 pm; LWC R , 1.0 
g/m 3 . 


Method 

c , in 

V, mph 

Reference 

36 

150 

Constant V 

18 

150 

Constant We w 

18 

175 

Constant We/, 

18 

193 

Constant Wes 

18 

197 

Constant Wet 

18 

212 

Constant Wei 

18 

212 

Constant We, 

18 

218 

Constant We,, + Const. 

18 

238 

Constant hf,i m /d (eq. (3.74)) 

18 

253 

Average V 

18 

265 

Constant Re a 

18 

332 

Constant V 

12 

150 

Constant We w 

12 

191 

Constant We ,, 

12 

224 

Constant Wes 

12 

230 

Constant Wet 

12 

* 

Constant WeL 

12 

260 

Constant We, 

12 

274 

Constant We,, + Const. hfn„/d 

12 

312 

Constant hf,, m /d (eq. (3.74)) 

12 

318 

Average V 

12 

* 

Constant Re a 

12 

* 


* No solution 


As Table V and figure 5.3.4 show, for half-size scaling, 
the scale velocity resulting from matching Wei falls in 
the range between constant Wes and average V that has 
shown reasonable scaling results in past studies. At this 
time, WeL seems to have greater experimental justifica- 
tion than other similarity parameters. The practical ad- 
vantage of using constant Wei is that the scale velocity- 
to-reference velocity ratio is simply the square root of the 
size ratio (eq. (3.71)). Furthermore, it leads to a more 
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Figure 5.3.4. Comparison of Velocity Ratio for Ten 
Methods to Find Scale Velocity. NACA 0012; c R , 36 in; 
t stl R, 6°F; S R , 40 pm; LWC R , 1.0 g/m 3 . 
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* No solution for some scale ratios with V R = 150 mph 


moderate scale velocity than some of the similarity pa- 
rameters that might be used. 

5.3.5. Size Scaling for Intercycle Ice Studies 

It is clear from the results of figure 5.3.1 (b) that for 14 - 
size scaling, the scale velocity cannot be set equal to the 
reference for adequate simulation of the horn angle or 
position for glaze ice. Flowever, these conclusions apply 
to significant ice-accretion times. For short accretion 
times, for which horns or other features have insufficient 
time to develop, constant- V determination of the scale 
velocity may be entirely adequate. Such is the case for 
intercycle-ice accretion, as one example. 

Intercycle ice is the ice that forms between cycles of an 
impulse deicing ice-protection system. The most com- 
mon form of impulse system is the pneumatic-boot de- 
icer in which pneumatic tubes inflate and deflate rapidly 
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at prescribed intervals to break ice from the leading-edge 
region of airfoils, tailplanes, or other critical surfaces. 
Aerodynamic forces then cause the ice to shed. Typi- 
cally, in such a process all the ice does not shed, and the 
surface will be left with some residual ice. After the boot 
deflates, ice will begin to form both on the clean portions 
of the surface and on the residual ice. The accretion that 
has formed by the next boot inflation cycle is termed the 
intercycle ice. 

Intercycle ice represents significant roughness on critical 
surfaces, and it is important to determine the aerody- 
namic effects of the roughness. Wind-tunnel testing is a 
relatively inexpensive and safe way to record intercycle 
ice shapes. If the desired wing section is too large for an 
icing tunnel, it is possible to perform a scale test to estab- 
lish the intercycle ice sizes and shapes. As for other ic- 
ing situations, shapes recorded by molding or by pencil 
tracing can be transformed into castings to be applied to 
wing sections for aerodynamic testing. 

In 2000 a preliminary study 75 was made in the Goodrich 
Icing Wind Tunnel (IWT) to determine if the size-scaling 
methods discussed in preceding sections of this manual 
could be applied to intercycle ice accretion. Intercycle 
ice shapes were recorded for a reference model and com- 
pared with shapes from a '/ 2 -size scale test. The desired 
reference model was a 72-in-chord NACA 23012 airfoil 
section, and the scale was a 36-in-chord 23012. Because 
the full 72-in-chord model would have been too large for 
the IWT, a hybrid-scaled 2,3 model was used in its place 
for the reference tests. The hybrid model consisted of a 
nose section with 72-in-chord NACA 23012 coordinates 
and a truncated afterbody. The resulting chord was 36 
in. The afterbody was designed so that the resulting hy- 
brid airfoil would have the same collection-efficiency 
curve at 2° AOA as the full 72-in-chord NACA 23012. 
This objective was met. The hybrid airfoil included a 
trailing-edge flap that was adjusted to correct the //curve 
for off-design angles of attack. Drop-trajectory and 
flowfield-analysis codes were used to determine the af- 
terbody coordinates at 2° AOA and to establish appropri- 
ate flap deflections for other angles of attack. 

The pneumatic deicer was mounted flush with the nose 
section of both models. It consisted of two 1.25-in-wide 
inflatable tubes at the leading edge of each model and 
one-in-wide tubes immediately aft on the upper and 
lower surfaces. The hybrid (reference) model used three 
1-in tubes on the upper surface and five on the lower 
surface, while the scale model had only one 1-in tube on 
the upper surface and two on the lower. The active por- 
tion of the deicers extended to the same non-dimensional 
chord-wise extent on both models. Thus, while the over- 


all active portions of the deicers were geometrically 
scaled, the individual inflatable tubes were not. 

Test conditions for the scale model were determined 
from the reference conditions using the Ruff method, 
described in section 5.3.3, with scale and reference ve- 
locities equated. The ratios of scale -to-reference boot 
cycle times were the same as the scale -to-reference total 
icing time. 

Two examples of the intercycle-ice scaling results from 
reference 75 are given in figure 5.3.5. In general, the 
scale tests reproduced the non-dimensional ice thickness 
of the reference tests at and near the leading-edge region. 
The non-dimensional size of ridges and features aft of the 
leading edge were also determined accurately by the 
scale tests. However, the chord-wise positions of these 
aft features were not well simulated. These features 
tended to form along the boot stitching; therefore, the 
authors’ speculated that if the pneumatic tubes for the 
deicers had been geometrically scaled, there might have 
been better agreement between the non-dimensional 
scale and reference locations of these features. Although 
further study and testing of scaling for intercycle ice ac- 
cretion is needed, the results demonstrate how constant- 
velocity scaling can be applied effectively when ice ac- 
cretion times are short. 

5.3.6. Acceptable Parameter Variations 

Even for glaze ice, the similarity parameters may not 
always need to be matched precisely between scale and 
reference to insure good scale ice shapes. In fact, in 
practice, uncertainties in test conditions make it impossi- 
ble to specify precisely the scale parameter values. Thus, 
it is helpful to know what tolerances are acceptable. This 
topic was reviewed in reference 19. There it was rec- 
ommended that uncertainties in test parameters should be 
minimized such that the scale-reference match of no and 
of A c should be within ±10%. If scale and reference j3 n 
do not match, the product /3oA c should match within 
±10%. Ice shape is not as sensitive to changes in veloc- 
ity as to those in no and A c , as shown in section 5.3.4. It 
is recommended, however, that We L ,s, the similarity pa- 
rameter from which scale velocity is found, be main- 
tained within ±15% of WeL,R. Additional testing of the 
effects of variations in this parameter on ice shape is 
needed to better establish acceptable variations in We L . 

5.3.7. Recommended Method to Scale Size 

The Ruff method reviewed in section 5.3.3, supple- 
mented with the constant- Wei approach (section 5.3.4) to 
calculate scale velocity, is recommended for simulating 
the main ice shape with a model scaled in size. This 
method requires matching scale and reference model 
non-dimensional geometry and AOA, as well as scale 
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Figure 5.3.5. Results of Size-Scaling Tests for Intercycle Ice. 75 Reference Model: 36-in-Chord Hybrid with 72-in- 
Chord NACA 23012 Nose Section and Truncated Trailing Edge; Scale Model: 36-in-Chord NACA 23012. 


and reference values of the similarity parameters K 0 , A c , 
n 0 , </>, 0 and We, for tunnels with altitude capability or K 0 , 
A„ no. We, and the choice of either ^ or 9 for sea-level 
tunnels. At this time, matching of We, appears to have 
the strongest analytical basis for finding scale velocity. 
The velocity resulting from use of this method is of the 
same magnitude as several other proposed approaches 
listed in Table V, and experimental evaluation of this 
method seems to validate it. As greater understanding of 
water-film effects is gained, this recommendation may 
change. 

Scale drop size can be found either from equation (3.18) 
or by matching K n from equation (3.8). For simulating 
the main ice shape, a great deal of flexibility is possible 
in choosing scale drop size, because, as shown in figures 
3.3.6 and 3.3.7, for the conditions investigated to date, K 0 
has no measurable effect on main ice shape. Neverthe- 
less, it is recommended that K„ be matched whenever 
possible, because it does affect the impingement limit. If 
scale and reference K 0 can not be matched, instead of 
simply matching A c , the product j3gA c should be matched 
as closely as possible. Furthermore, if K 0 is not matched 


and impingement limits are important additional tests 
should be performed at the reference J3 0 , A c and n 0 to 
determine these limits. For scaled impingement-limit 
tests, the parameters <f i (or 9) and We,, necessary for 
faithful scaling of main ice shapes, can be ignored. 

For scaling tests to determine main ice shape, there is a 
high probability that the parameters (/> and 9 do not have 
to be matched to the reference values, but additional test- 
ing is required to confirm that this simplification is valid. 
Until such testing is completed, it is recommended that 
these parameters continue to be used. 

In some cases, the scale total temperature found by 
matching either ^ or 9 is above freezing. References 35 
and 46 showed that when total temperatures exceeded 
27°F (-3°C) ice shapes are distorted enough that a good 
match of the reference is not possible. Thus, when scale 
total temperature is greater than 27°F, (f> and 9 must be 
abandoned and a scale temperature selected such that the 
total temperature is below 27°F (-3°C). LWC must then 
also be adjusted to maintain constant n 0 . 
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For very short accretion times for which ice features do 
not have time to form, simplifications to the above rec- 
ommendations are possible. For example, a constant- 
velocity approach can be substituted for constant- WeL, as 
discussed in section 5.3.5. For short times (small values 
of A c ) the primary concern should be to match scale and 
reference (it A,- 

The matching of all the parameters required to find scale 
test conditions necessitates solving several equations 
simultaneously. The programming of solutions to this 
system of equations is not difficult, and is simplified by 
the use of mathematical software such as Mathcad or 
Mathematical. 

5.3.8. Practical Limitations to Size Scaling 

This section reviews the practical constraints to size scal- 
ing imposed by facility restrictions as well as by natural 
physical limitations. Both types of restriction impose an 
upper limit on the reference-to-scale model size ratio. 
Restrictions on test-condition scaling have not been ex- 
plored sufficiently to define those limits; therefore, this 
discussion will concentrate on size-scaling limits. 

Figure 5.3.6 shows the relationships between scale con- 
ditions and scale ratio. The values in this figure were 
calculated by applying the method recommended in this 
manual; that is, the Ruff method supplemented with con- 
stant WeL to determine scale velocity. The validity of 
this method has not been tested at scale ratios greater 
than 3.4, but there does not appear to be any fundamental 
reason that it can’t be applied to higher scale ratios. 
Testing scaling methods at higher scale ratios is difficult 
because, to perform such tests, both the reference and 
scale conditions have to fall within the tunnel calibrated 
operating map. 

Figure 5.3.6. (a) gives the effect of scale ratio on scale 
drop size for reference drop sizes of 20, 30 and 40 pm. 
The calculations were made for a reference velocity of 
150 mph, but the effect of velocity on scale drop size is 
not significant. It can be seen from the figure that for a 
facility with a minimum calibrated drop size of 14 pm, 
the maximum scale ratio is about 3.5 for a reference drop 
size of 40 pm. For reference drop sizes of 30 and 20 pm, 
the maximum scale ratio is 2.5 and 1.5, respectively. If 
the facility calibration is extended down to an MVD of 
10 pm, scale ratios of close to 5.5 could be tested with 
the 40-pm reference drop size. 

Figure 5.3.6. (b) shows the effect of scale ratio on scale 
velocity for reference velocities of 100, 150 and 200 
mph. For this illustration, the constant- WeL method of 
finding scale velocity was used. If the reference velocity 
is 200 mph and the maximum test-section velocity is 350 
mph, scale ratios only as high as 3 can be tested. Obvi- 




Scale Ratio, c F /c s 
(b) Effect on Scale Velocity. 

Figure 5.3.6. Effect of Scale Size Ratio on Scale Test 
Conditions. Scaling using Ruff Method with Constant 
We L . c R , 72 in; T R , 490 R. 

ously, lower reference velocities permit higher scale ra- 
tios for the same maximum test-section velocity. For a 
facility with a maximum test-section velocity of 200 
mph, the scale ratio is limited to approximately 2 when 
simulating a reference velocity of 150 mph and 4 for a 
reference velocity of 100 mph. 

There may be other physical limitations. The minimum- 
MVD problem might be eased somewhat by extending 
the calibration of spray systems to smaller drop sizes. 


NASA/CR— 2004-2 12875 


50 





Such calibration depends on accurately measuring small 
drops, which cannot yet be done with certainty. How- 
ever, even if possible, this approach is limited by the 
minimum drop size that will impinge on a model (see 
figure 3.3.2.) For smaller and smaller cloud MVD’s it is 
more and more likely that the smallest portions of the 
cloud will fall into the regime in which no impingement 
occurs. This situation makes it questionable whether a 
very small-MVD cloud can be used reliably to simulate a 
cloud with a larger MVD. Testing to determine practical 
minimum MVD’s has not been done and such testing is 
needed if large reference-to-scale model size ratios are 
desired. 

Velocities are limited by physical constraints, as well. 
As the velocity increases, the risk of aerodynamic dam- 
age to the scale ice shape during accretion increases. 
Frequent shedding of feathers and even of small portions 
of glaze horns has been observed in tests with velocities 
over 300 mph. Furthermore, as velocities increase above 
350 mph compressibility effects may need to be consid- 
ered, and the scaling assumptions that Mach number can 
be ignored may no longer be valid. Thus, it appears that 
a velocity somewhere around 350 mph is probably a 
practical maximum for scaling applications, even if the 
facility is capable of higher speeds. 

In summary, then, because of both facility and funda- 
mental physical constraints, scale model sizes cannot be 
much smaller than to Vi to 14 that of the reference, de- 
pending on reference conditions and facility capabilities. 

6. Recommended Future Studies in Scaling 

This manual concentrated on scaling for simple geome- 
tries and two-dimensional wing sections. Future studies 
need to address such applications as swept wings, rotor- 
craft and scaling for ice-protection systems. For small 
sweep angles, swept-wing accretions may be close 
enough to unswept that the scaling methods in this man- 
ual might be applicable. Testing is required to verify this 
approach, however, and to date none has been done. 


For rotorcraft, it may be possible to divide the rotor into 
several radial regions according to rotating speed. De- 
termination of the local ice shape would come from tests 
performed at various speeds to represent each radial loca- 
tion. Near the tip, speeds will be too high for icing facili- 
ties, and lower speeds could be found by matching a 
similarity parameter such as We L . Matching any of the 
Weber numbers to simulate a high speed with a lower 
speed necessitates the use of a larger chord. Therefore, 
this approach has practical limits. Furthermore, simula- 
tion of high-speed shedding is not addressed by this ap- 
proach. Finally, the Mach number may need to be in- 
cluded in the scaling methodology for regions at and near 
the tip. 

In addition to ice shedding caused by aerodynamic 
forces, shedding due to centrifugal forces can have a 
significant effect on ice shape for rotating blades. A 
shedding analysis needs to include issues such as varying 
adhesion of ice to the surface and local stresses on the ice 
for different conditions. 

Scaling for thermal ice -protection systems will require 
the addition of terms to the energy balance equation of 
section 3.5 to represent the ice-protection system. For 
such applications of scaling, there is no ice shape to be 
concerned about, so, in addition to thermal similarity 
parameters, it is only necessary to match the product PoA c 
to insure the correct water catch is simulated. This is not 
difficult to do, but scaling methods for thermal ice pro- 
tection need to be evaluated. 

Intercycle-ice scaling studies, initiated in the work of 
reference 75, need to be continued. Future experiments 
in this area should consider the effects of pneumatic-tube 
size relative to airfoil size, inflation-deflation schedules 
and inflation pressure. Reference 75 used a hybrid- 
scaled model as equivalent to the full-size reference; in 
future studies, if possible, a model with the full-size 
chord and geometry should be used for clear-cut com- 
parisons with the scale results. 
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Appendix A 

Resources for Scaling Calculations 


This Appendix contains a listing of all the properties and 
relationships needed to find the scaled test conditions 
that will simulate an icing encounter with given reference 
conditions. A sample size-scaling case will be solved 
step by step using the Ruff method with constant We L , 
and the solution to an L IFC-scaling case using the Olsen 
method will also be demonstrated. Finally tables of 
sample reference and scale test conditions will be given 
for both types of scaling. 

A.l. Properties of Air, Water and Ice and General 
Relationships 

This section gives the constants or expressions used to 
find the properties of air and water. These property val- 
ues are those used in scaling calculations at NASA 
Glenn. They will be listed alphabetically. Values or 
equations will be given in pairs, with the first of each pair 
in English units and the second in S.I. All symbols used 
here appear in the Nomenclature list, Section 1 of the 
Manual. References for this Appendix appear at the end 
of the Appendix. 

c r Specific Heat at Constant Pressure The spe- 

cific heat of air at a pressure of 14.5 psia over the range - 
40°F to 40°F is given by NIST 76 as 


c = . 240 —^^— 
p ' a lbmR 


= .240 


cal 

gK 


(A.l) 


At this pressure, c pa varies with temperature only in the 
5 th decimal place. The variation with pressure for these 
temperatures is less than .2% over the range 6 - 14.5 
psia. Therefore, a constant value for c pa can be used for 
all icing scaling situations. 

The ratio, y, of the constant-pressure specific heat to the 
constant-volume specific heat for air will be discussed 
below. 


Data for the specific heat of water over the range -40°F 
to 40°F were given by Pruppacher and Klett. 29 A fit to 
their data gives the expression used for water on the sur- 
face of the model: 


Rtn 

c DW =1.0074 + 

p ’ ws lbmR 

Rtn 

2.559xl0“ 5 j(T s -491. 67R) 2 


lbmR 


-1.0074-HV 

gK 
-5 cal 


(A.2) 


8.29x10 


gK 


(T s -273. 15 K) 2 


3 v s 


Note that for the first form of equation (A.2) the absolute 
surface temperature, T s , is in Rankine degrees, while in 
the second it must be in Kelvin. 


k„ . Thermal Conductivity of Air A fit to the NIST air 
thermal conductivity data 76 at a pressure of 14.5 psia for 
the range of temperature from -40 to 40°F is 


k n = -.008529-^ + . 001016 Bt “ - 


hr ft R 


= -12.69- 


cal 

hrmK 


- + 2.029- 


hrftR 
cal 


(A.3) 


hrmK 




film 


For this temperature range k a varies with pressure only 
by approximately .1% for pressures from 6 to 14.5 psia. 
Therefore, equation (A.3) can be used without regard for 
pressure for icing scaling calculations. As the units in 
equation (A.3) indicate, the absolute temperature must be 
in Rankine or Kelvin, respectively. As indicated, the 
thermal conductivity of air is evaluated at the film tem- 
perature, defined by equation (A. 13), when calculating 
the convective heat-transfer coefficient from equation 
(A.33). 

D „ Diffusivity of Water vapor The water- 

vapor diffiisivity in air was given by the following ex- 
pression in Pruppacher and Klett. 29 


A, 


= 2.271x10 


4 All 

( T ^ 

* film 

1.94 

^ 14.696 

s I 

v 491. 67R j 


V Pst , 


= .211 


cm 2 

f 't \ 

* film 

1.94 

^1.0132xl0 5 Pa" 

s 

v 273.15K y 


l Pst ) 


(A.4) 


Note that the diffusivity is evaluated at the film tempera- 
ture, equation (A. 13). 

p. Air Pressure For the IRT the settling chamber is 
vented to atmosphere, so the total pressure everywhere in 
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the tunnel is approximately atmospheric. For the altitude 
in Cleveland, Ohio, a typical value for atmospheric pres- 
sure is 14.5 psia. Thus, for scaling calculations, 


Ptot = 14.5psia 
= 9.997 xl0 4 Pa 


(A.5) 


is used for 1RT testing for both scale and reference con- 
ditions. 

The static pressure in the absence of compressibility is 


Pst = ' 


1 + 


Ptot 
V 2 

j 


The compressible form of equation (A.6) is 

Ptot 


Pst ~ ' 


1 + ^ 1 M 2>- 1 


(A.6) 


(A.7) 


The static pressures found from equations (A.6) and 
(A.7) differ by approximately 1.5% at an airspeed of 400 
mph. At 350 mph, a practical maximum airspeed for the 
IRT, the difference is less than 1%, and it continues to 
decrease with speed. Thus, for typical icing conditions, 
equation (A.6) provides an adequate approximation. 

For some tunnels, the test section, rather than the settling 
chamber, is vented to the atmosphere. In that situation, 
the test-section static pressure is the ambient pressure 
and the total pressure is calculated from equation (A.6). 

p„ „ Vapor Pressure of Water Prappacher and 

Klett gave a curve fit for the saturation pressure of vapor 
over water for the range of -58°F to 32°F (-50°C to 0°C). 
It is: 


Pw 

~ a 0 + 





r 

r 






a 3 + 

AT 

q.j + AT 

0-2 + AT 

AT 

' a 4 + 5 


V 

V 

LAI 

V 

AT(a 5 +ATa 6 ) / 




JJ) 


(A- 8) 


where 


AT = T -491. 67R 
= 7’-273.15K 


(A.9) 


a 3 = 6.5919xl0" 7 psi/R 3 
c/ 4 = 4.1880xl0" 9 psi/R 4 
a 5 = 1.5613xl0" n psi/R 5 
a ,5 = 2.6169x 10- 14 psi/R 6 


= 2.6506xl0" 2 Pa/K 3 
= 3.0312xl0" 4 Pa/K 4 
= 2.0341xl0" 6 Pa/K 5 
= 6.1368xl0" 9 Pa/K 6 


Equation (A.8) can be used to find both the vapor pres- 
sure at the surface, p, m ,, and that in the atmosphere, p m by 
substituting the appropriate absolute temperature, T, in 
equation (A.9). For glaze ice the surface is at the freez- 
ing temperature of water, so AT in either system of units 
is 0, and p mv = a 0 . To determine the atmospheric vapor 
pressure the ambient static temperature, T st , in Rankine 
or Kelvin units, respectively, should be used. 


In section 3.5 sublimation was neglected in the energy 
balance because it is small compared with the evapora- 
tion term. Flowever, for completeness, it will be dis- 
cussed briefly here. The driving force for sublimation is 
the difference between the partial pressure of vapor at the 
surface and that in the atmosphere. The atmospheric 
vapor pressure can be determined from equations (A.8), 
(A.9) and the coefficients an to a 6 . At the surface, the 
partial pressure of vapor is the saturation pressure of wa- 
ter vapor over ice. Pmppacher and Klett’s fit for the 
saturation pressure of vapor over ice for the range -58°F 
to 32°F (-50°C to 0°C) is: 


Pwi = b o + 



c 

f 

( h 3 + 

AT 

bj +AT 

b 2 +AT 

AT 

r b 4+ 

K AT(b 5+ ATb 6 1 


V 

V 

V 




(A. 10) 


with AT given by equation (A.9) and coefficients 


Z? 0 = .088606 psia 
bi = 4.0568x1 O' 3 psia/R 
b 2 = 8.4426x1 O' 5 psia/R 2 
Z> 3 = 1.0386xl0' 6 psia/R 3 
Z> 4 = 8.0476x1 O' 9 psia/R 4 
b 3 = 3.7141xl0" n psia/R 5 
be = 7.8412xl0" 14 psia/R 6 
R„ , Gas Constant for Air 


= 610.92 nt/m 2 
= 50.347 Pa/K 
= 1.8860 Pa/K 2 
= 4.1762xl0" 2 Pa/K 3 
= 5.8247xl0" 4 Pa/K 4 
= 4.8388xl0" 6 Pa/K 5 
= 1.8388xl0" 8 Pa/K 6 
The air gas constant is 


with coefficients: 
ao= .088586 psi 
ai = 3.5748x10° psi/R 
a 2 = 6.3964x 10" 5 psi/R 2 


= 610.78 Pa 
= 44.365 Pa/K 
= 1.4289 Pa/K 2 


R =53.34 


lbfft 

lbmR 


= .287 


ntm 


(A.11) 
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(A. 18) 


T, Absolute Temperature Absolute temperatures were 
calculated from the temperatures in °F and °C in the 
standard way: 

T = t + 459.67R =t + 273.15K (A.12) 

By convention, the degree symbol, °, is omitted from 
Kelvin unit, 77 K. To be consistent with this practice, the 
degree symbol is omitted from the Rankine unit, R, as 
well. 

Tfi/m, Film Temperature The film temperature is used 
to evaluate air properties for convective heat transfer 
parameters and for diffusivity. The film temperature is 
the average of the free-stream static temperature and the 
surface temperature: 

Tfi,m=-(T S +T s ,) (A.13) 


This temperature can be interpreted as approximately the 
mean temperature of the boundary layer. 


t a Surface T emperature For glaze ice the existence of 
liquid water on the surface requires the surface tempera- 
ture to be the freezing temperature of water. Thus, 


t s = 32°F 
= 0°C 


(A. 14) 


and 


T s = 491.67R 
= 273. 15K 


(A. 15) 


7 %, Total Temperature The total temperature is re- 
lated to the static temperature for incompressible flow by 
the following expression: 


c p,a 

r = — 

C v,a 

= 1.40 

From reference 76 at -40°F y varies from 1.402 to 1.404 
as pressure varies from 6 to 14.5 psia. At 6 psia y de- 
creases from 1.402 to 1.401 as temperature increases 
from -40 to 40°F. At 14.5 psia / decreases from 1.404 to 
1.402 for the same temperature increase. All can be 
rounded off to 1.40. Therefore, the value in equation 
(A. 18) is representative of yfor the range of temperatures 
and pressures of interest to icing. 

A, Latent Fleat From Pmppacher and Klett, the latent 
heat released when water freezes is 


A f =143.46 — + .485— (r-491.67R)- 
J lbm lbmR 7 ’ 

1.39xl0~ 3 BtU „ (T-491.67R) 2 


lbmR 

= 79.7 — + .485— (7-273. 15K)- 
g gK^ 7 

2.5x10“ 3 -^-(7-273.15K) 2 

gK“ 


(A. 19) 


For scaling calculations, the temperature in this expres- 
sion is the surface temperature in absolute units. Thus, 
incorporating equation (A. 15) results in 


A f =143.46 


Btu 

lbm 


.79.7A 


for glaze ice. 


(A.20) 


Pmppacher and Klett give the following expression for 
the latent heat of vaporization: 


T lo t=T st+: f- (A. 16) 

^ c p,a 

where c pjl is the specific heat of air, equation (A.l). For 
compressible flow, the relation is 


4, =1075 


= 597.3 


Btu f 491.67R 
lbm 1. 7 

cal f 273. 15K 

g l T 


(A.21) 


^=^(i+V m2 J (A17) 

However, there is less than .01% difference between the 
results obtained from equations (A. 16) and (A. 17) for 
speeds up to 400 mph. Note that equation (A. 16) can be 
applied to either absolute or ordinary temperatures, while 
equation (A. 17), in the form given, can only be used with 
absolute temperatures. 

r. Ratio of Specific Heats for Air The ratio of specific 
heats is 


Where the exponent E is also a function of temperature: 

E = .167 + 2.039 xlO -4 — 

l (A.22) 

= . 167 + 3. 670xl0“ 4 — 

K 

As in the evaluation of Af the temperature to be used in 
these expressions is the surface temperature. Thus, for 
glaze ice 

A v =1075— =597.3— (A.23) 

lbm g 
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Because the heat balance equations, in which Af and A v 
appear, are only of interest for glaze ice, equations 
(A.20) and (A.23) are not needed for rime scaling calcu- 
lations. 


u. Viscosity The viscosity of air as a function of 
temperature can be obtained from reference 76. Over the 
range -40°F to 40°F these values were fit to the following 
expression: 


Pa =■ 


10 


-5 


.18995 + 333.18 


lbm 
R'l fts 


10 


-4 


(A. 24) 


.12764 + 124.38 


K 


cms 


The air viscosity varies by less than .1% over the pres- 
sure range of 6 - 14.5 psia for the given temperature 
range, so the effect of pressure has been ignored. De- 
pending on the application the temperature T may be 
either the free-stream static, T st , or the film temperature 
from equation (A. 13). 


The viscosity of water used in these calculations was the 
value at the freezing point of water: 


ured using a standard technique on three occasions using 
different samples. 42 The value that has been used consis- 
tently in scaling calculations for the IRT is 

a wa =4.45x10 -3 ^- =65^ (A.29) 

it cm 

The two other values reported in reference 42 were 
3.33xlO' 3 lbfft and 4. 1 Ox 10 _3 lbf/ft (48.6 and 59.8 
dyne/cm). 

Surface tension changes when surfactant is added to the 
water. Reference 42 reported surface tensions of water 
measured with the addition of two surfactants at various 
concentrations. For all combinations, surface tension 
varied only from 1.94x10° Ibf/ft to 2.07xl0' 3 lbf/ft. For 
calculations of Weber numbers for IRT tests with surfac- 
tant addition, 2.06xl0" 3 lbf/ft (30 dyne/cm) was selected 
as representative, and this value has been used consis- 
tently. 

A.2. Similarity Parameters and Other Terms 

The definitions of the similarity parameters and other 
formulations needed to perform icing scaling calculations 
will be repeated here from the main text; they are listed 
alphabetically. 


p w = 1.20xl0“ 3 — = .01786—^— (A.25) 

fts cms 

A Density The density of air as a function of 

temperature and pressure was found using the gas-law 
expression 


LWCVt 
Pi d 

LWC VP 0 c pws 

K 


(A.30) 

(A.31) 


Pa 


Pst 

R a P s t 


(A.26) 


where p st and T st are the static pressure and temperature 
of the airflow. 


Ice density varies depending on the type of ice. For the 
scaling calculations a fixed value was assumed: 


-,-lbm 

Pi =51.2 — — =.917 


ft 3 ' cm 3 


(A. 27) 


There is an uncertainty of perhaps as much as 10% on 
this value, but as long as the same value is used for both 
reference and scale calculations, the scaling results 
should be insensitive to this uncertainty. 

The density of water at the freezing point was used. It is 

Pw =62A^ =1-0— (A.28) 

ft~ cm 


gro . Surface Tension The surface tension of water 
against air for the IRT demineralized water was meas- 


h = 


d = 2 r le 
k,Nu , 


1 film 

~d~ 


/ n.1/2 

r LWC^ 


Re: 


K 


'Pry 


Sc, 


p.a \ a J 


(A. 3 2) 
(A.33) 

(A.34) 

(A.35) 


When computing the values of h c and h G from equations 
(A.33) and (A.35) the film temperature, equation (A. 13), 
should be used to evaluate k a , Nu a , Pr a and Sc a . 


K = 


pjy_ 

18 (i p a 


(A.36) 
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(A.39) 

(A.40) 

(A.41) 

(A.42) 

(A.43) 
(A. 44) 
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Equation (A.51) is based on equation (0.18) of the main 
text with at= 3.38 substituted. 


F 2 A 

0 = T-T t + ^A v 

2C p,a K 


Pww P tot Pw 

T T p 

St tot Jr St 

1 P tot Pww 

, 0.622 T tot T st 


A tabulation of range parameter as a function of drop 
Reynolds number was given by Langmuir and Blodgett. 
A fit to their data has been used in scaling calculations 
for the IRT. It is 


0.8388 + 0.0014837^ +0.1847 ^Re~ s 

t/2 


tj) -Tf -T st - 


(A.54) 


Note that when computing the value of Re a for use in 
equation (3.33), the air density, p a , and viscosity, p a , 
should be evaluated at the film temperature, equation 
(A.13). 


Rex =- 


(A.45) 

(A.46) 
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We L =- 
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(A. 50) 
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Appendix B 

Sample Scaling Calculations 


B.l. Sample Calculations for Size Scaling (Ruff 
Method with Constant We L ) 

An example of calculating the scale test conditions to 
simulate a given reference icing encounter will be de- 
tailed in this section. This illustration uses the Ruff 
method supplemented with constant We L to find the scale 
velocity. 

B. 1 .1. Choose Reference Conditions The fol- 

lowing reference conditions will be used for an NACA 
0012 airfoil: 

<y; = 36in =91.4 cm 

t st , R = 5.8 °F = -14.6 °C 

p tot = 14.5 psia =1.00xl0 5 Pa 

F/;=150mph =67.1 m/s 

8 r = 40 pm 
LWC r = 1 g/m 3 
t r = 15.1 min 


jU aR = 1.104xl0" 5 lbm/ft s = 1.643xl0' 4 g/cm s 
RcgR = 213.5 
A r = 0.2592 

p w = 62.43 lbm/ft 3 = 1 .000 g/cm 3 

K r = 12.56 
K-o,r = 3.350 
Po.r = 0.789 

Values of the parameters ( b , <f> and 0) needed for the 
freezing fraction and the freezing fraction itself will now 
be found. 

c p ,wsr = 1.007 Btu/lbm R = 1.007 cal/g K 
Ttnm,R= 478.6 R = 265.9 K 

k„,R = 0.0137 Btu/hr ft R = 20.4 cal/hr m K 
Pr a . R - 0.712 

The Prandtl number was evaluated at the film tempera- 
ture. 


Calculate Values for Reference Properties and 


Parameters The leading-edge radius for the NACA 
0012 airfoil is given by Abbott and von Doenhoff 78 as 


1.58% of chord. Thus, with 


d R = 0.0316c ff 

= 1.138in =2. 890cm 


From equation (A. 12) the reference static temperature is 


p a ,R = 0.0816 lbm/ft3 = 1.308xl0 3 g/m 3 
Re aR = 1.542x10 s 

This is the model Reynolds number based on twice the 
leading-edge radius of curvature and properties at the 
free-stream static temperature. The following Reynolds 
number uses properties evaluated at the film temperature, 
as does the Nusselt number: 


T stR = 465.47 R = 258.6 K 

Substituting V R and T R into equation (A.38) gives the 
reference Mach number, 

M aM = 0.2080 

The total temperature and static pressure can now be 
determined from equations (A. 17), (A. 12) and (A.7), 
respectively: 

Tto,.R = 469.5 R = 260.8 K 

t, 0 ,R= 9.8°F = -12.3°C 

p stR = 14.07 psi = 9.700xl0 4 Pa 

The accumulation parameter is (eq. (3.26)) 

A c , r = 2.293 

The properties and parameters needed to determine K n 
will be calculated next, using the static temperature, T stR , 
for temperature-dependent properties: 

p aR = 0.08163 lbm/ft 3 = 1.308 g/m 3 


Rsfdm.R— 1.468xl0 5 
Nu r = 381.3 

h c , R = 55.10 Btu/hr ft 2 R = 2.690x 10 5 cal/hr m 2 K 
b R = 0.7134 

(j> R = 25.24 R = 14.02 K 

/7 Wj R = 0.0886 psia =610.8 Pa 

p w , R = 0.0288 psia =198.4 Pa 

D vR = 2.249xl0" 4 ft 2 /s = 2.089xl0" 5 m 2 /s 
Sc a , R = 0.6316 

h GR = 0.06909 lbm/ft 2 s = 337.3 g/m 2 s 
9 r = 34.62 R = 19.23 K 

A fJl = 143.5 Btu/lbm = 79.70 cal/g 
yio,R = 0.5 1 8 
We L . R = 2.00xl0 6 
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B.1.3. Calculate Scale Properties, Parameters and Test 
Conditions A scale model chord of 10.5 in was 
chosen for this example. 

d s - 0.03 16c s 

= 0.332 in =0.843 cm 

Matching the Wei between r eferenc e and scale gives the 
simple relationship V s = V R y Jc R / c s . Then 


V s = 277.7 mph = 124.2 m/s 

The scale Wei calculated from this velocity is, of course, 

We L , s = 2.00xl0 6 

The static temperature can be found in terms of <j) s by 
rearranging equation (A.54): 

T st,S ~Tf -0S ~Vs / 2c p,ws ■ 

Because t/>s = (j>R is one of the similarity equalities, this 
equation can be solved to give: 


T st , s = 463. 1 R = 257.3 K 

#** = 3.5 °F = -15.8 °C 

M a , s = 0.386 

T m ,s= 477.0 R = 265.0 K 


tm.s = 17.3 °F = -8.2 °C 

p stS = 13.1 psia =9.02xl0 4 Pa 


The scale drop size can now be determined. Equation 
(A.51) can be used for this purpose, or that equation can 
provide a first guess of 8 S to try to match K os with K 0R . 
Solving equation (A.51) for Ogives 


8 S = 14.5 pm 


K 0 ,s can now be calculated to compare with K 0R \ 
p a ,s= 0.0765 lbm/ft 3 = 1.23xl0 3 g/m 3 
JUq,s = 1 . 1 Ox 1 0" 3 lbm/ft s = 1.64xl0" 4 g/cm s 

Ress= 134.9 
^-Stokes = -3 14 


K s = 10.5 


K n .s= 3.39 

This value of K n , s is within about 1% of the K 0 ,r calcu- 
lated above. For practical applications, this is close 
enough. For the purposes of illustration, however, the 
calculation will be continued with a new estimate for 8 S : 

<%= (14.5 pm )(K or /Kos) = 14.3 pm 

Rcfts ~ 133.3 


K s = 10.3 


K 0 ,s= 3.33 


Finally, linear interpolation between the two previous 8s 
and their corresponding K os gives a third estimate of drop 
size: 


8 S = 14.5pm+(l4.3pm-14.5pm) 


3.35-3.39 

3.33-3.39 


= 14.4pm 


Rets* 133.9 
^-Stokes = 0.315 
K s = 10.4 
K„. s = 3.35 
Po,s= 0.789 

The final value of K 0 ,s matches K 0 r to two decimal 
places. Note that the final value of 8s is only .1 pm less 
than the original estimate from equation (A.51). This 
difference is much less than the experimental uncertainty 
in drop MVD, so it can be seen that equation (A.51) pro- 
vides entirely adequate values of scale drop size for this 
example. This adequacy of equation (A.51) has also 
been observed for a variety of scaling situations and ref- 
erence test conditions. 


The scale L WC can now be determined by matching the 
scale and reference freezing fraction, no. As with the 
solution for drop size, an iterative procedure is required. 
The first guess is to match L WC S and L WC R . Thus, 

LWC S = 1 g/m 3 


c P ,wsS = 1-007 Btu/lbmR 
Tf,i m ,s = 477.4 R 
k a , s = 0.01367 Btu/hr ft R 
7^ = 0.7118 
Pa,s~ 0.07626 lbm/ft 3 
Re a ,s = 7.81 lxlO 4 
Re/iim, s = 7.42 1 x 1 0 4 
Nu a ,s = 271.1 
h c , s = 134.0 Btu/hr ft2 R 
bs= 0.5431 
(j)s= 25.24 R 


= 1.007 caFgK 
= 265.2 K 
= 20.35 cal/hr m K 

= 1.222xl0 3 g/m 3 

= 6.5 44x1 0 5 cal/hr m 2 K 
= 14.02 K 


(Recall that the scale value of (f> was set to the reference 
value earlier in this section to find the scale static tem- 
perature.) 


k/k stokes — 0.316 


Pww,s = .08859 psia =610.8 Pa 
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p w s = .02585 psia 
D vS = 2 ,394x 1 0" 4 ft 2 /s 
Sc a , s = 0.6339 
h G , s = 0.1676 lbm/ft 2 s 
<% = 27.72 R 
/lg S = 143.5 Btu/lbm 
«o,s = 0.536 


= 178.2 Pa 
= 2.224xl0" 5 m 2 /s 

= 818.5 g/m 2 s 
= 15.42 K 
= 79.7 cal/g 


This first-guess value of L WCs leads to a scale freezing 
fraction a little higher than the reference value of 0.518, 
so a second estimate will be made by forming a simple 
ratio with the desired n 0 ,s of 0.518: 


LWC S 


t g ( 0-536^1 

m 3 v 0.5 1 8 J 


= 1.035-^- 
m 3 


Both (f) and //as well as the properties needed to calculate 
tin are independent of LWC. Therefore, with the new 
estimate for L WC it is only necessary to reevaluate bp. 


b s = 0.5621 


Wo.s — 0.524 

As this freezing fraction is still too high, a third estimate 
of LWC will be made by extrapolating from the two pre- 
vious results: 


g 


LWC S = 1.000-2-± 1.037 -2- -1.000 -2- 


g | 0.536-0.518 


0.536-0.524 


= 1.052- 


m~ 


This new value of L WCs then gives the following values 
of b s and n 0 ,s- 

b s = 0.5714 


It only remains to determine the scale icing time. This is 
found by setting A cS = A C _ R and solving equation (3.26) 
for the time: 


Ts = 2.26 min 


Table B.l. gives a number of sample reference conditions 
with the corresponding size scaling test conditions found 
using the Ruff Method supplemented with constant We/ 
to find the scale velocity. Reference and scale similarity 
parameters are also given in the table. Case 1 is the size 
scaling example presented in section B.l. Numbers in 
the table were generated with a scaling computer pro- 
gram and were rounded off for presentation. 

B.2. Sample Calculations for LWC Scaling (Olsen 
Method) 

The Olsen method maintains scale and reference accu- 
mulation parameter and freezing fraction the same for 
fixed values of model size, velocity and drop size. Scale 
LWC can be chosen arbitrarily, then scale temperature is 
calculated to provide a match of freezing fraction, and 
scale icing time is calculated by matching accumulation 
parameter. 

B.2.1. Choose Reference Conditions The fol- 

lowing reference conditions will be used, with, again, an 
NACA 0012 airfoil: 


c R = 36 in 
t s t,R~ 19.6 °F 
p tot = 14.5 psia 
V R = 200 mph 
S R = 20 pm 
LWC r = 0.25 g/m 3 
t r = 25 min 


= 91.4 cm 
= -6.9 °C 
= 1.00x10 s Pa 
= 89.4 m/s 


viz is — 0.5 18 

Thus, with LWCs = 1.052 g/m’ the scale and reference 
freezing fractions match. This final value of L WCs is 
only about 5% greater than the original estimate, while 
the experimental uncertainty in LWC is about ±10%. 
Furthermore, for analytical freezing fractions of this 
magnitude, experimental ice shapes do not show discern- 
able differences for freezing fraction variations less than 
10%. Therefore, if the estimated LWCs results in a scale 
freezing fraction within about 10% of the reference 
value, that is usually close enough, and the above se- 
quence of calculations to refine the original estimate of 
LWCs would not have been necessary in a practical ap- 
plication. 


B.2.2, Calculate Values for Reference Properties and 
Parameters The combination of reference velocity, 
drop MVD and L WC given in section B.2.1 is outside the 
current range of calibrated conditions for the IRT. Con- 
sequently, in order to perform a test, it is necessary to 
choose a new LWC that falls within the IRT calibrated 
envelope. Scale model size, velocity and drop size will 
be set to the same values as the reference; thus, K 0 is 
already matched between scale and reference. This pa- 
rameter only needs to be evaluated for the purpose of 
finding /3 0 , which is needed to find n 0 . Freezing fraction 
and accumulation parameter will be calculated for the 
reference case. 

d R = 0.03 16c r 

= 1 . 1 3 8 in =2. 890cm 
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Table B.l 


Examples of Size Scaling 
Ruff Method with Constant We L 

NACA 0012 Airfoils 


Case 

Type 

in 

hb 

°F 

hob 

°F 

v, 

mph 

3 

(am 

LWC, 

g/m 3 

h 

min 

K 0 

% 

A c 

n 0 

b 

°F 

9 

°F 

Re, 

10 4 

We h 

10 6 

M 

1 

Ref 

36.0 

5.8 

9.8 

150.0 

40.0 

1.00 

15.1 

3.35 

0.789 

2.29 

0.518 

0.714 

25.2 

34.6 

15.4 

2.00 

0.208 

n 

Scale 

10.5 

3.5 

17.3 

277.7 

14.4 

1.05 

2.3 

3.35 

0.789 

2.29 

0.518 

0.572 

25.2 

27.8 

7.8 

2.00 

0.386 

2 

Ref 

72.0 

0.0 

4.0 

150.0 

40.0 

1.00 

25.0 

1.72 

0.675 

1.90 

0.559 

0.862 

31.0 

41.9 

31.5 

4.00 

0.209 

II 

Scale 

36.0 

-1.0 

7.1 

212.1 

22.5 

1.09 

8.1 

1.72 

0.675 

1.90 

0.559 

0.804 

31.0 

39.0 

21.7 

4.00 

0.296 

3 

Ref 

72.0 

0.0 

4.0 

150.0 

30.0 

1.00 

25.0 

1.12 

0.583 

1.90 

0.613 

0.745 

31.0 

41.9 

31.5 

4.00 

0.209 

n 

Scale 

36.0 

-1.0 

7.1 

212.1 

17.0 

1.09 

8.1 

1.12 

0.583 

1.90 

0.613 

0.695 

31.0 

39.0 

21.7 

4.00 

0.296 

4 

Ref 

72.0 

0.0 

4.0 

150.0 

20.0 

1.00 

25.0 

0.62 

0.438 

1.90 

0.743 

0.560 

31.0 

41.9 

31.5 

4.00 

0.209 

II 

Scale 

36.0 

-1.0 

7.1 

212.1 

11.4 

1.09 

8.1 

0.62 

0.438 

1.90 

0.743 

0.522 

31.0 

39.0 

21.7 

4.00 

0.296 

5 

Ref 

72.0 

0.0 

7.2 

200.0 

40.0 

1.00 

25.0 

2.03 

0.706 

2.53 

0.471 

1.054 

30.3 

38.8 

41.0 

7.10 

0.279 

II 

Scale 

36.0 

-1.7 

12.6 

282.8 

22.3 

1.01 

8.8 

2.03 

0.706 

2.53 

0.471 

0.917 

30.3 

33.7 

27.7 

7.10 

0.395 

6 

Ref 

72.0 

0.0 

7.2 

200.0 

30.0 

1.00 

25.0 

1.33 

0.620 

2.53 

0.507 

0.926 

30.3 

38.8 

41.0 

7.10 

0.279 

II 

Scale 

36.0 

-1.7 

12.6 

282.8 

16.8 

1.01 

8.8 

1.33 

0.620 

2.53 

0.507 

0.805 

30.3 

33.7 

27.7 

7.10 

0.395 

7 

Ref 

72.0 

0.0 

7.2 

200.0 

20.0 

1.00 

25.0 

0.74 

0.481 

2.53 

0.592 

0.718 

30.3 

38.8 

41.0 

7.10 

0.279 

ii 

Scale 

36.0 

-1.7 

12.6 

282.8 

11.3 

1.01 

8.8 

0.74 

0.481 

2.53 

0.592 

0.624 

30.3 

33.7 

27.7 

7.10 

0.395 

8 

Ref 

72.0 

0.0 

4.0 

150.0 

40.0 

1.00 

25.0 

1.72 

0.675 

1.90 

0.559 

0.862 

31.0 

41.9 

31.5 

4.00 

0.209 

ii 

Scale 

24.0 

-1.9 

10.2 

259.8 

16.1 

1.10 

4.4 

1.72 

0.675 

1.90 

0.559 

0.745 

31.0 

36.2 

17.3 

4.00 

0.363 

9 

Ref 

72.0 

0.0 

4.0 

150.0 

30.0 

1.00 

25.0 

1.12 

0.583 

1.90 

0.613 

0.745 

31.0 

41.9 

31.5 

4.00 

0.209 

ii 

Scale 

24.0 

-1.9 

10.2 

259.8 

12.2 

1.10 

4.4 

1.12 

0.583 

1.90 

0.613 

0.644 

31.0 

36.2 

17.3 

4.00 

0.363 

10 

Ref 

72.0 

0.0 

4.0 

150.0 

20.0 

1.00 

25.0 

0.62 

0.438 

1.90 

0.743 

0.560 

31.0 

41.9 

31.5 

4.00 

0.209 

II 

Scale 

24.0 

-1.9 

10.2 

259.8 

8.2 

1.10 

4.4 

0.62 

0.438 

1.90 

0.743 

0.484 

31.0 

36.2 

17.3 

4.00 

0.363 

11 

Ref 

72.0 

0.0 

7.2 

200.0 

40.0 

1.00 

25.0 

2.03 

0.706 

2.53 

0.471 

1.054 

30.3 

38.8 

41.0 

7.10 

0.279 

II 

Scale 

24.0 

-3.4 

18.1 

346.4 

15.8 

0.92 

5.2 

2.03 

0.706 

2.53 

0.471 

0.779 

30.3 

28.6 

21.6 

7.10 

0.485 

12 

Ref 

72.0 

0.0 

7.2 

200.0 

30.0 

1.00 

25.0 

1.33 

0.620 

2.53 

0.507 

0.926 

30.3 

38.8 

41.0 

7.10 

0.279 

II 

Scale 

24.0 

-3.4 

18.1 

346.4 

11.9 

0.92 

5.2 

1.33 

0.620 

2.53 

0.507 

0.684 

30.3 

28.6 

21.6 

7.10 

0.485 

13 

Ref 

72.0 

0.0 

7.2 

200.0 

20.0 

1.00 

25.0 

0.74 

0.481 

2.53 

0.592 

0.718 

30.3 

38.8 

41.0 

7.10 

0.279 

ii 

Scale 

24.0 

-3.4 

18.1 

346.4 

8.1 

0.92 

5.2 

0.74 

0.481 

2.53 

0.592 

0.530 

30.3 

28.6 

21.6 

7.10 

0.485 

14 

Ref 

72.0 

10.0 

17.2 

200.0 

40.0 

1.00 

25.0 

2.03 

0.706 

2.53 

0.315 

1.055 

20.3 

26.0 

39.6 

7.10 

0.276 

II 

Scale 

24.0 

6.6 

28.1 

346.4 

15.8 

0.77 

6.3 

2.03 

0.706 

2.53 

0.315 

0.647 

20.3 

15.9 

20.8 

7.10 

0.480 

15 

Ref 

72.0 

10.0 

17.2 

200.0 

30.0 

1.00 

25.0 

1.32 

0.620 

2.53 

0.339 

0.926 

20.3 

26.0 

39.5 

7.10 

0.276 

II 

Scale 

24.0 

6.6 

28.1 

346.4 

12.0 

0.77 

6.3 

1.32 

0.620 

2.53 

0.339 

0.569 

20.3 

15.9 

20.8 

7.10 

0.480 

16 

Ref 

72.0 

10.0 

17.2 

200.0 

20.0 

1.00 

25.0 

0.74 

0.480 

2.53 

0.397 

0.717 

20.3 

26.0 

39.5 

7.10 

0.276 

II 

Scale 

24.0 

6.6 

28.1 

346.4 

8.1 

0.77 

6.3 

0.74 

0.480 

2.53 

0.397 

0.440 

20.3 

15.9 

20.8 

7.10 

0.480 

17 

Ref 

72.0 

15.0 

22.2 

200.0 

30.0 

1.00 

25.0 

1.32 

0.620 

2.53 

0.252 

0.926 

15.3 

19.1 

38.8 

7.10 

0.275 

II 

Scale 

24.0 

11.6 

33.1 

346.4 

12.0 

0.60 

8.0 

1.32 

0.620 

2.53 

0.252 

0.443 

15.3 

9.1 

20.5 

7.10 

0.477 
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T s ,,r = 479.3 R 

= 266.3 K 

Substituting V R and T R into equation (A. 38) gives the 
reference Mach number. 

M a , R = 0.2734 


Ttot.R = 486.4 R 

= 270.2 K 

Uot.R = 26.8°F 

= -2.9°C 

Pst.R = 13.77 psi 

= 9.492xl0 4 Pa 

A c , r = 1.265 


p a , R = 0.07754 lbm/ft 3 

= 1.242 g/m 3 

p a .R = 1.130xl0" 5 lbm/ft s 

= 1.681xl0" 4 g/cm s 

Re m = 132.1 


X R = 0.3 167 


p w = 62.43 lbm/ft 3 = 1.000 g/cm 3 

K R = 4.09 


Ko,r = 1.381 


Po.r= 0.629 


Values of the parameters b, $ and <9 used in the freezing 
fraction expression and the freezing fraction itself will 
next be determined. 

Cp.wsR = 1-007 Btu/lbm R 

= 1.007 cal/g K 

Tfiim.R = 485.5 R 

= 269.7 K 

k a ,R = 0.01386 Btu/hr ft R 
Pr a , R = 0.7115 
Re a>R = 1.909xl0 5 
Re 'film , r = 1 . 8 67 x 1 O' 

Nu a ,R = 429.8 

= 20.62 cal/hr m K 

h c , R = 62.84 Btu/hr ft 2 R 

= 3.068xl0 5 cal/hr m 2 K 

b R = 0.1662 


<f> R = 10.69 R 

= 5.94 K 

Pww.r = 0.08859 psia 

= 610.8 Pa 

p WtR = 0.05293 psia 

= 365.0 Pa 

D vR = 2. 3 62x1 O' 4 ft 2 /s 

= 2.194xl0" 5 m 2 /s 

Sc aR = 0.6308 

h aR = 0.07884 lbm/ft 2 s 

= 384.9 g/m 2 s 

d R = 12.38 R 

= 6.88 K 

A fR = 143.5 Btu/lbm 

= 79.70 cal/g 

nog = 0.598 



B.2.3. Choose Scale Conditions All scale conditions 
are the same as the reference except for the scale liquid- 
water content, temperature and time. L WC S was selected 
to fall within the IRT calibrated envelope. 

c s =36in =91.4 cm 

p, ol = 14.5 psia =1.00xl0 5 Pa 

V s = 200 mph =89.4 m/s 

8 S = 20 pm 

LJVCs= 0.60 g/m 3 

B.2.4. Calculate Scale Properties, Parameters and Test 
Conditions The strategy in the Olsen method is to 
adjust temperature so that the scale and reference freez- 
ing fractions match for the given d , V, 5 and L WC. 

Because conditions affecting the modified inertia pa- 
rameter are the same as for the reference situation, the 
leading-edge collection efficiency will have the same 
value as the reference. Therefore, 


As =0.629 

c P ,wsS = 1-007 Btu/lbm R = 1.007 cal/g K 

An initial guess for the value of scale temperature can be 
made by setting it equal to the reference value: 


T st , s = 479.3 R 
M a , s = 0.2734 
T, 0 ,, s = 486.4 R 
p st .s = 13.77 psi 
Tjiim,s = 485.5 R 
k a .s = 0.01386 Btu/hr ft R 
^0,5 = 0.7115 
p a ,s= 0.07754 lbm/ft 3 
Re a , s = 1.909xl0 5 
Rcjiims = 1.867x10' 

Nu a ,s = 429.8 

h Ci s = 62.84 Btu/hr ft2 R 

b s = 0.3988 

</> s = 10.69 R 

Pww,s = 0.08859 psia 

/Vs =0.05293 psia 

D vS = 2.362x 10" 4 ft 2 /s 


= 266.3 K 

= 270.2 K 
= 9.492xl0 4 Pa 
= 269.7 K 
= 20.62 cal/hr m K 

= 1.242xl0 3 g/m 3 


= 3.068xl0 5 cal/hr m 2 K 

= 5.94 K 
= 610.8 Pa 

= 365.0 Pa 
= 2.194xl0" 5 m 2 /s 
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Sc a , s = 0.6308 

h G , s = 0.07884 lbm/ft 2 s = 384.9 g/m 2 s 
0 S = 12.38 R = 6.877 K 

A fs = 143.5 Btu/lbm = 79.7 cal/g 
n 0 , s = 0.2931 


Because this scale freezing fraction is lower than the 
reference value of 0.598, the estimated scale temperature 
was too high. It will be arbitrarily decreased by 5 R and 
the calculations repeated: 

T st ,s= 474.3 R = 263.5 K 

M a , s = 0.2748 

T m ,s= 481.4 R = 267.5 K 

p stiS = 13.76 psi = 9.486xl0 4 Pa 


Tfiim.s = 483.0 R = 268.3 K 

k a , s = 0.01380 Btu/hr ft R = 20.54 cal/hr m K 


Pr a . s = 0.7116 

p a , s = 0.07832 lbm/ft 3 = 1.255xl0 3 g/m 3 
Re a , s = 1.944xl0 5 
Rsfiim, s = 1.883xl0 5 
Nu ai s= 431.8 

h c , s = 62.86 Btu/hr ft2 R = 3.069xl0 5 cal/hr m 2 K 
b s = 0.3987 

0s = 15.69 R = 8.72 K 


Pmv,s = 0.08859 psia =610.8 Pa 

p w> s = 0.04264 psia = 294.0 Pa 

A^ = 2.339xl0" 4 ft 2 /s = 2.173xl0" 5 m 2 /s 


Sc a , s = 0.6313 

h as =0.07883 lbm/ft 2 s =384.9 g/m 2 s 
6 S = 19.64 R = 10.91 K 

A fs = 143.5 Btu/lbm = 79.7 cal/g 


>io,s = 0.4560 


Note that many of the values changed inconsequentially 
when the temperature was changed from the initial esti- 
mate. The calculations can be simplified if desired by 
approximating the parameter bs with 


b S 


, LWC S 
b R 

lwc r 


Using this constant value for bs for all the iterations, only 
0s an d 0s need to be re-evaluated for each new tempera- 
ture estimate. However, for this sample calculation, the 
complete set of evaluations will be shown. 

The next temperature estimate can be found by a linear 
extrapolation of the two previous temperatures with their 
corresponding freezing fractions: 


rp rp , n 0,R n 0,Sl (rp rp 

1 st,S ~ 1 st,Sl H l- z rt,S2 ~ 1 st,Sl 

n 0,S2~ n 0, 


T s = 469.9 R 

M a , s = 0.2761 

7(ot5 = 477.1 R 

t st ,R = 10.24 °F 

t,„,.R= 17.40 °F 

Pst ,s= 13.75 psi 

Tfiim.s = 480.8 R 

k a , s =0.01375 Btu/hr ft R 

Pr a ,s = 0.7116 

p a ,s= 0.07901 lbm/ft 3 

Re a . s = 1.975xl0 5 

Re fi im, s= 1.898x10' 

Nu a , s = 433.5 

h c .s = 62.88 Btu/hr ft 2 R 

bs ~ 0.3989 

0s = 20.06 R 

Pww.s — 0.08859 psia 

Pw.s = 0.03515 psia 

D vS = 2.3 1 9x 1 0 4 ft 2 /s 

*5^,5 = 0.6317 

h Gi s = 0.07882 lbm/ft 2 s 

0 S = 25.65 R 

A fs = 143.5 Btu/lbm 


= 261.1 K 

= 265.1 K 
= -12.09 °C 
= -8.11 °C 
= 9.481xl0 4 Pa 
= 267.1 K 
= 20.46 cal/hr m K 

= 1.266xl0 3 g/m 3 


= 3.070xl0 5 cal/hr m 2 K 

= 11.14 K 
= 610.8 Pa 

= 242.4 Pa 
= 2.154xl0" 5 m 2 /s 

= 384.8 g/m 2 s 
= 14.25 K 
= 79.7 cal/g 


rio,s = 0.5927 

This value of freezing fraction is within 1% of the refer- 
ence value of 0.5982; consequently, T stS = 469.9 R 
(261.1 K) will be accepted as the scale temperature. 
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Finally, spray time is found by matching scale and refer- 
ence accumulation parameter: 

A c ,s = A c r = 1 .265 


T S ~ A c,S 


dsPi 

LWC S V S 


= 10.4 min 


Table B.II gives the reference and scale information for 
several L WC scaling cases using the Olsen method. The 
detailed sample calculations offered in section B.2. were 
for case 1 of the table. 

B.3. Procedures for Scaling Drop Size and Tempera- 
ture 

B.3.1. Drop-Size The calibrated envelope of Appendix- 
C cloud conditions for the 1RT currently covers a range 
of drop sizes from 14 to 50 pm MVD. Figure 4.1.1. 
showed the Mod-1 and Standard nozzle calibrated LWC- 
vs -8 envelopes for an airspeed of 200 mph (89 m/s). In 
section 4.1 some of the facility limitations were dis- 
cussed. In this section, examples will be given for scal- 
ing drop size when desired conditions fall outside 
the capability of the facility to test them. These desired 
conditions may themselves be the result of a size-scaling 
calculation. For example, in Table B.l the scale drop 


sizes calculated for cases 4, 7, 9, 10, 12, 13, 15, 16 and 
17 are all below the 14-pm minimum calibrated drop size 
for the 1RT. However, in section 3.3.3 it was demon- 
strated that the main ice shape was unaffected by a 
change in drop size from 20 to 55 pm. If this insensitiv- 
ity to drop size holds true for smaller drops, it should be 
possible to make at least small changes to the scale drop 
size and still obtain scale ice shapes that adequately rep- 
resent the reference shape. Using case 4 in Table B.l as 
an example, the procedure is: 

(1) Increase 8$ from 11.4 to 15 pm without changing c$. 
Vs or L WCs . The scale cloud conditions now fall within 
the IRT envelope. 

(2) Calculate the new fio.s- Note that this parameter will 
no longer match the reference value, and impingement 
limits will therefore be different for scale and reference 
accretions. For this example Po.s increases from 0.438 to 
0.540. 

(3) Calculate the new scale icing time that will produce a 
match between the scale and reference products fin.sAc.s 
and Po.rA Ci r. 

(4) Calculate a new scale temperature such that nn.s 
matches the reference value. This step is necessary be- 
cause the increase in /3 0 .s would otherwise decrease the 


Table B.II 

Examples ofLWC Scaling 
Olsen Method 


NACA 0012 Airfoils 


Case 

Type 

C, 

in 

*sh 

°F 

Hoh 

°F 

v, 

mph 

3 

jam 

UVC, 

g/m 3 

T, 

min 

K 0 

fio, 

% 

A c 

n 0 

b 

°F 

e, 

°F 

Re, 

10 4 

We h 

10 6 

M 

1 

Ref 

36.0 

19.6 

26.8 

200.0 

20.0 

0.25 

25.0 

1.38 

0.629 

1.26 

0.598 

0.166 

10.7 

12.4 

1.91 

3.55 

0.273 

n 

Scale 

36.0 

10.0 

17.2 

200.0 

20.0 

0.60 

10.4 

1.38 

0.629 

1.26 

0.598 

0.399 

20.3 

25.9 

1.98 

3.55 

0.276 

2 

Ref 

36.0 

23.7 

27.7 

150.0 

20.0 

0.25 

25.0 

1.16 

0.590 

0.95 

0.548 

0.134 

7.4 

9.4 

1.44 

2.00 

0.204 

II 

Scale 

36.0 

17.4 

21.4 

150.0 

20.0 

0.55 

11.4 

1.16 

0.590 

0.95 

0.548 

0.294 

13.7 

18.9 

1.48 

2.00 

0.205 

3 

Ref 

36.0 

21.8 

25.9 

150.0 

40.0 

0.25 

25.0 

3.35 

0.789 

0.95 

0.548 

0.179 

9.2 

12.3 

1.45 

2.00 

0.205 

n 

Scale 

36.0 

7.8 

11.8 

150.0 

40.0 

0.82 

7.6 

3.35 

0.789 

0.95 

0.548 

0.586 

23.3 

32.1 

1.53 

2.00 

0.208 

4 

Ref 

36.0 

20.3 

27.5 

200.0 

20.0 

0.25 

25.0 

1.38 

0.629 

1.26 

0.548 

0.166 

10.0 

11.3 

1.90 

3.55 

0.273 

II 

Scale 

36.0 

15.2 

22.3 

200.0 

20.0 

0.45 

13.9 

1.38 

0.629 

1.26 

0.548 

0.299 

15.1 

18.8 

1.94 

3.55 

0.275 

5 

Ref 

36.0 

18.4 

25.5 

200.0 

40.0 

0.25 

25.0 

3.96 

0.812 

1.26 

0.548 

0.215 

11.9 

14.2 

1.92 

3.55 

0.274 

II 

Scale 

36.0 

5.4 

12.6 

200.0 

40.0 

0.70 

8.9 

3.96 

0.812 

1.26 

0.548 

0.601 

24.9 

32.0 

2.01 

3.55 

0.277 

6 

Ref 

36.0 

16.4 

27.6 

250.0 

20.0 

0.25 

25.0 

1.59 

0.658 

1.58 

0.548 

0.197 

12.9 

12.8 

2.34 

5.55 

0.343 

II 

Scale 

36.0 

12.6 

23.7 

250.0 

20.0 

0.38 

16.4 

1.59 

0.658 

1.58 

0.548 

0.300 

16.8 

18.4 

2.38 

5.55 

0.344 

1 

Ref 

36.0 

14.5 

25.6 

250.0 

40.0 

0.25 

25.0 

4.52 

0.829 

1.58 

0.548 

0.249 

14.9 

15.7 

2.36 

5.55 

0.344 

II 

Scale 

36.0 

4.5 

15.7 

250.0 

40.0 

0.55 

11.4 

4.52 

0.829 

1.58 

0.548 

0.547 

24.8 

29.1 

2.45 

5.55 

0.347 

8 

Ref 

36.0 

12.0 

28.1 

300.0 

20.0 

0.25 

25.0 

1.78 

0.681 

1.90 

0.548 

0.228 

16.1 

14.1 

2.76 

7.99 

0.413 

n 

Scale 

36.0 

8.7 

24.8 

300.0 

20.0 

0.35 

17.9 

1.78 

0.682 

1.90 

0.548 

0.319 

19.5 

18.6 

2.79 

7.99 

0.415 

9 

Ref 

36.0 

10.0 

26.1 

300.0 

40.0 

0.25 

25.0 

5.05 

0.842 

1.90 

0.548 

0.281 

18.2 

16.8 

2.78 

7.99 

0.414 

n 

Scale 

36.0 

1.6 

17.8 

300.0 

40.0 

0.48 

13.0 

5.05 

0.842 

1.90 

0.548 

0.540 

26.5 

27.8 

2.86 

7.99 

0.418 
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scale freezing fraction. For case 4, t s/S has to be reduced 
from -1°F to -5.9°F to give an n os of 0.743, the reference 
value. 

Figure B.l. illustrates this example, showing a map of 
calibrated cloud conditions for the IRT and the specific 
L WC - VI V D combination for this example. Because the 
IRT undergoes improvements from time to time the 
cloud calibration is updated periodically. Consequently, 
the envelopes shown here may not be the same as the 
current calibration for the IRT. Figure B.l. (a) shows the 
cloud envelope for the reference velocity of 150 mph (67 
m/s). The portion of the map obtained with mod-1 noz- 
zles is given with a solid line and for the standard nozzles 
with a dashed line. The solid point on the figure repre- 
sents the reference drop size and liquid water content. 

For the scale velocity of 212 mph (95 m/s) the tunnel 
calibrated cloud contracts to that of figure B.l. (b). The 
solid circle marked ‘1’ indicates the scale drop size and 
LWC as given in Table B.L This condition fell outside 
the nozzle range for the facility, and therefore could not 
be tested without special calibration. The second symbol 
gives the final solution, with S s increased to place the 
condition within the calibrated envelope. 

This example involved a relatively small change in 8 S 
and was therefore simple to adjust to the facility enve- 
lope. This is not always the case. For example, for the 
i/s-scale case 10 in Table B.l 8 S was 8.2 pm. While a 
value within the facility envelope can be found and the 
possibility exists that substituting such a value may be 
valid, tests of the effect of such a large change in MVD 
for this magnitude of drop size have not been made; 



Drop Median Volume Diameter, pm 


therefore, such substitutions carry an unknown level of 
risk. 

Alternately, for situations like case 10, it is possible in 
principle to modify the original set of reference condi- 
tions by using test-condition scaling. We know from 
figures 3.3.6 and 3.3.7 that an encounter with a cloud 
MVD of 40 pm should give the same ice shape as one 
with a drop size of 20 pm. Therefore, the reference con- 
ditions for case 10 can be replaced by a set of conditions 
with the same model size, temperature, velocity and liq- 
uid-water content as shown in Table B.l, but with 8= 40 
pm instead of 20. Because J3 0 will now be different from 
that shown in Table B.l for case 10, the scale time must 
be altered to provide the same /3oA c as that for the case-10 
reference. This new set of reference conditions is now 
scaled to reduce the model size to !4 of the reference 
using size scaling. Scaling by this model size ratio from 
a drop size of 40 pm gives a scale drop size of 16.1 pm, 
which is within the tunnel cloud envelope. This ap- 
proach can only be used to simulate main ice shapes, and 
additional tests with matched //, are needed if impinge- 
ment limit information is desired. Drop-size substitution 
has not yet been fully tested but appears to be promising 
based on the results shown in figures 3.3.6 and 3.3.7. 
However, until further validation is obtained, this method 
should be used cautiously. 

Although not illustrated by any of the cases in Table B.l, 
it is also possible that scale L WC may fall above or be- 
low the calibrated envelope of the facility. However, this 
situation is easily remedied by applying the Olsen 
method to scale L WC to a value within the nozzle array 



Drop Median Volume Diameter, pm 


(a) Envelope for V= 150 mph (67 m/s) (Reference). (b) Envelope for V= 212 mph (95 m/s) (1/2-Scale). 

Figure B.l. IRT Cloud Envelope with Conditions for Scaling Case 4, Table B.l. 
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operating map. An example of this approach was given 
in section B.l. 

B.3.2. Temperature Cases 14 - 17 of Table B.l all 
had scale total temperatures in excess of 27°F. It was 
noted in section 3.5.13 that scale temperatures higher 
than this have been observed to produce ice shapes that 
do not match the reference shapes well. Therefore, for 
each of these cases it would be necessary to find a tem- 
perature below 27°F for testing. The following examples 
illustrate the process. 

For the conditions of case 14 it is only necessary to select 
a new scale static temperature a little lower than the re- 
sult determined by the Ruff scaling method. With c s , V s 
and Ss unchanged, the modifications are: 

(1) Arbitrarily, set the scale static temperature a few 
degrees lower than the value originally found. Thus, t^s 
= 5°F (-15°C). 

(2) Find a new LWC$ for the new static temperature to 
maintain nos = n o,R- Note that for this example, jBo.s for 
the modified scale conditions is the same as the original, 
which also matched j3d,r. The result is LWC S = 0.93 
g/m 3 . 

(3) Calculate a new scale icing time so that A C: s = A cR . 
Then, zy= 5.2 min. 


(4) Determine the new scale total temperature from 
equation (A. 17). The static temperature change from 6.6 
to 5°F has a negligible effect on the Mach number, so 
M aS does not have to be re-calculated. t tolS is then found 
to be 26.5°F (-3°C). 

For case 15 both a lower scale static temperature and a 
higher scale drop size are needed compared with the 
scale values in Table B.L Figure B.2. (a) shows the facil- 
ity cloud envelope with the reference condition as a solid 
circle. The initial scale cloud conditions are given as 
point 1 in figure B.2. (b). The adjustments are made as 
follows with c s and Cs unchanged: 

(1) Set the scale static temperature a few degrees lower 
than the value in Table B.l to reduce the total tempera- 
ture below 27°F. As for case 14 above, t st s is lowered to 
5°F (-15°C). 

(2) Determine the new scale total temperature from 
equation (A. 17). The Mach number is effectively un- 
changed from the value in Table B.l because the static 
temperature undergoes only a small change. Then we 
find t m s = 26.5°F. This is sufficiently below 27°F to be 
acceptable. 

(3) Choose a slightly larger Ss such that it is within the 
facility cloud envelope. In this example, a value of 15 
pm is used. 




Drop Median Volume Diameter, pm Drop Median Volume Diameter, pm 

(a) Envelope for V= 200 mph (89 m/s) (Reference). (b) Envelope for V= 346 mph (155 m/s) (1/3-Scale). 

Figure B.2. IRT Cloud Envelope with Conditions for Scaling Case 15, Table B.II. 
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(4) Calculate the new J3 0 ,s corresponding to S s = 15 pm. 
The result is 0.691. 

(5) Find the new LWC S for the new static temperature 
and § s such that n 0i s = n 0 , R . The result is LWC S = 0.83 
g/nrl 

(6) Reevaluate the scale accumulation parameter so that 

Po,sA c ,s = Pi>,rA c ,r ■ A c ,s = 2.29. 

(6) Calculate a new scale icing time from the final A cS . 
The result is t s = 5.2 min. This completes the calculation 
of all the revised scale conditions. 

The revised scale cloud conditions for case 15 are shown 
as point 2 on figure B.2. (b). 

Adjustments to the scale conditions for case 17 in Table 
B.l will be discussed next. Figure B.3. (a) shows the 
facility cloud envelope with the reference condition as a 
solid circle. The scale cloud conditions from Table B.l 
are indicated as point 1 in figure B.3. (b). The scale total 
temperature for this case is above freezing, and the scale 
drop size is outside the calibrated envelope for the 1RT. 
Thus, again it is necessary to reduce the scale static tem- 
perature and increase the scale drop size from the scale 
values in Table B.l. 

The adjustments are made in the same way as for case 
15, except that larger changes are necessary. Again, c s 
and V s are unchanged: 


(1) Set the scale static temperature sufficiently lower 
than the value in Table B.l. that the total temperature will 
be below 27°F. t slS is set to 5°F (-15°C). Because the 
velocity and static temperature for this example is the 
same as for case 15, the Mach number and total tempera- 
ture will also be the same. The resulting total tempera- 
ture, t, otiS = 26.5°F. As before, this is sufficiently below 
27°F to be acceptable. 

(2) Choose a new S s such that it is within the facility 
cloud envelope. In this example, a value of 20 pm was 
expected to be large enough. It is desirable that this 
value not be any greater than necessary, but because we 
don’t yet know the value of LWC S that corresponds with 
the new temperature and drop size this initial guess for S s 
may need to be adjusted later. 

(3) The p os corresponding to S s = 20 pm is calculated to 
be 0.767. 

(4) The L WC S required to keep n 0 ,s = n QR with the new 
static temperature and 8 S is 1.42 g/m 3 . 

(5) The scale accumulation parameter, A cS , is found to 
be 2.04 by matching p 0 ,sA c , s to P m A cS . 

(6) Calculate the scale icing time from A cS . t s = 2.7 
min. 




Drop Median Volume Diameter, pm 


Drop Median Volume Diameter, pm 


(a) Envelope for V= 200 mph (89 m/s) (Reference). (b) Envelope for V= 346 mph (155 m/s) (1/3-Scale). 

Figure B.3. IRT Cloud Envelope with Conditions for Scaling Case 17, Table B.II. 
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These new scale cloud conditions for case 17 are shown 
as point 2 on figure B.3. (b). As can be seen from the 
figure, the S s - LWC S combination was still just outside 
the calibrated map of the tunnel. Consequently, a 
slightly larger Ss will be selected and the calculations 
repeated. The static temperature does not have to be 
changed. The procedures are the same as steps (2) - (6) 
above, and the resulting values are: 

(1) S s = 22 pm 

(2) Po. s = 0.789 

(3) LWC S = 1.38 g/m 3 

(4) A c , s , = 1.99 

(5) zj=2.7min 

The final scale cloud conditions for case 17 are plotted as 
point 3 on figure B.3. (b). The 8 S - LWC S combination 
now is within the calibrated envelope of the facility. 
This solution requires a fairly large modification to Ss in 
a range of drop sizes for which the effect on ice shape 
has not been tested. Thus, while a solution is possible, 
we do not know how well the final result will simulate 
the desired accretion. 
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